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Symbols 


A 

spectral shape function for 
TBL-TE noise, eqs. (35) 
and (36) 

a, a o 

parameters of shape function 
A , eqs. (37) and (38) 

B 

spectral shape function for 
separation noise, eqs. (41) 
and (42) 

b, bo 

parameters of shape function 
B, eqs. (43) and (44) 

c 

chord length, m 

Co 

medium speed of sound, m/s 

D h 

directivity function for TE 
noise (high-frequency limit), 
eq. (Bl) 

D f 

directivity function for trans- 
lating dipole (low-frequency 
limit), eq. (B2) 

F( St) 

universal spectral shape 
function, eq. (18) 

/ 

frequency, Hz 

Ci 

spectral shape function for 
LBL-VS noise, eq. (57) 

g 2 

R c -dependence for LBL-VS 
noise peak amplitude, eq. (58) 

C 3 

angle dependence for G<i 
function, eq. (60) 

c 4 

peak level function for G 5 , 
eq. (74) 

C 5 

spectral shape function for TE 
bluntness noise, eqs. (75)-(82) 

Gij(f) 

cross-spectrum between 
microphone IvL and M 7 , 
Pa 2 /Hz 

H 

tunnel height, m 

h 

TE thickness (degree of 
bluntness), m 

K,K u AKuK 2 

constants, defined by eqs. (18), 
(47), (48), and (49) 

L 

span, m 

L' 

sectional lift of blade, lift per 
unit span 


c 

turbulence correlation scale, m 

l 

spanwise extent of tip vortex 
at TE, m 

M 

Mach number, Ujc 0 

Mmax 

maximum Mach number in tip 
vortex formation region 

R C 

Reynolds number based on 
chord length, cU/v 

Rij(r) 

cross-correlation between 
microphones Mi and Mj, Pa 2 

r 

distance from source to 
observer, m 

S(f ) 

spectrum of self-noise, Pa 2 /Hz 

St, St 1 , St,, St2 

Strouhal numbers defined for 
TBL-TE and separation noise 
scaling, section 5.1. 

st', st; 

Strouhal numbers defined 
for LBL-VS noise scaling, 
section 5.2. 

St" 

Strouhal number defined for 
tip vortex formation noise, 
section 5.3. 

St'" 

Strouhal number defined for 
TE-bluntness- vortex-shedding 
noise, section 5.4. 

t 

time, s 

u 

free-stream velocity, m/s 

Uc 

convection velocity, m/s 

X 

streainwise axis, see fig. B3, m 

y 

lateral axis, m 

z 

vertical axis, m 

C*T1P 

angle of attack of airfoil tip to 
oncoming flow, deg 

a TIP 

corrected angle of attack of 
airfoil tip, eq. (66), deg 

at 

airfoil angle of attack refer- 
enced to tunnel streamwise 
axis, deg 

a* 

effective aerodynamic angle 
of attack, corrected for open 
wind tunnel effects, deg 

6 

boundary-layer thickness, m 

6* 

boundary- layer displacement 
thickness, m 


V 



r tip vortex strength, m 2 /s 0 

0 angle from source streamwise 

axis x to observer, see fig. B3, 



deg 

1/3 

9 

boundary- layer momentum 



thickness, m 

Abbreviations: 

V 

kinematic viscosity of medium, 

BL 


m 2 /s 

LBL 

T 

time delay, s 

LE 

$ 

angle from source lateral axis y 

LHS 


to observer, see fig. B3, deg 

Mi 

0 

cross-spectral phase angle, deg 



angle parameter related to sur- 
face slope at TE, section 5.4, 

OASPL 


deg 

RHS 

Subscripts: 


SPL 

avg 

average 

TBL 

e 

retarded coordinate 

P 

pressure side of airfoil 

TE 

UTRC 

3 

suction side of airfoil 


TIP 

tip of blade 

VS 

TOT 

total 

2D 

a 

angle dependent 

3D 


vi 


for <5 q, <5q, and 9q is for airfoil 
at zero angle of attack, refer- 
ence value 

1 /3-octave spectral 
presentation 

boundary layer 
laminar boundary layer 
leading edge of airfoil blade 
left-hand side 

microphone number i for i = 1 
through 9, see fig. 4 

overall sound pressure level, 
dB 

right-hand side 

sound pressure level, spectrum, 
dB (re 2 x 10 -5 Pa) 

turbulent boundary layer 

trailing edge of airfoil blade 

United Technologies Research 
Center 

vortex shedding 
two-dimensional 
three-dimensional 
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Summary 

An overall prediction method has been developed for the self-generated noise 
of an airfoil blade encountering smooth flow. Prediction methods for individual 
self-noise mechanisms are semiempirical and are based on previous theoretical 
studies and the most comprehensive self-noise data set available. The specially 
processed data set, most of which is newly presented in this report, is from a 
series of aerodynamic and acoustic tests of two- and three-dimensional airfoil 
blade sections conducted in an anechoic wind tunnel. Five self-noise mecha- 
nisms due to specific boundary- layer phenomena have been identified and mod- 
eled: boundary-layer turbulence passing the trailing edge, separated-boundary- 
layer and stalled-airfoil flow, vortex shedding due to laminar- boundary- layer 
instabilities, vortex shedding from blunt trailing edges, and the turbulent vor- 
tex flow existing near the tips of lifting blades. The data base, with which the 
predictions are matched, is from seven NACA 0012 airfoil blade sections of dif- 
ferent sizes (chord lengths from 2.5 to 61 cm) tested at wind tunnel speeds up 
to Mach 0.21 (Reynolds number based on chord up to 3 x 10 6 ) and at angles of 
attack from 0° to 25.2°. The predictions are compared successfully with pub- 
lished data from three self-noise studies of different airfoil shapes, which were 
tested up to Mach and Reynolds numbers of 0.5 and 4.6 x 10 6 , respectively. 
An application of the prediction method is reported for a large-scale-model he- 
licopter rotor and the predictions compared well with data from a broadband 
noise test of the rotor, conducted in a large anechoic wind tunnel. A computer 
code of the methods is given for the predictions of 1/3-octave formatted spectra. 



1. Introduction 

Airfoil self-noise is due to the interaction be- 
tween an airfoil blade and the turbulence produced 
in its own boundary layer and near wake. It is 
the total noise produced when an airfoil encounters 
smooth nonturbulent inflow. Over the last decade, 
research has been conducted at and supported by 
NASA Langley Research Center to develop funda- 
mental understanding, as well as prediction capabil- 
ity, of the various self-noise mechanisms. The interest 
has been motivated by its importance to broadband 
helicopter rotor, wind turbine, and airframe noises. 
The present paper is the cumulative result of a se- 
ries of aerodynamic and acoustic wind tunnel tests 
of airfoil sections, which has produced a comprehen- 
sive data base. A correspondingly extensive semi- 
empirical scaling effort has produced predictive 
capability for five self-noise mechanisms. 

1.1. Noise Sources and Background 

Previous research efforts (prior to 1983) for the 
broadband noise mechanisms are reviewed in some 
detail by Brooks and Schlinker (ref. 1). In fig- 
ure 1, the subsonic flow conditions for five self-noise 
mechanisms of concern here are illustrated. At high 
Reynolds number R c (based on chord length), turbu- 
lent boundary layers (TBL) develop over most of the 
airfoil. Noise is produced as this turbulence passes 
over the trailing edge (TE). At low R c , largely lam- 
inar boundary layers (LBL) develop, whose instabil- 
ities result in vortex shedding (VS) and associated 
noise from the TE. For nonzero angles of attack, the 
flow can separate near the TE on the suction side of 
the airfoil to produce TE noise due to the shed tur- 
bulent vorticity. At very high angles of attack, the 
separated flow near the TE gives way to large-scale 
separation (deep stall) causing the airfoil to radiate 
low-frequency noise similar to that of a bluff body in 
flow. Another noise source is vortex shedding occur- 
ring in the small separated flow region aft of a blunt 
TE. The remaining source considered here is due to 
the formation of the tip vortex, containing highly tur- 
bulent flow, occurring near the tips of lifting blades 
or wings. 

1A.1. Turbulent- Boundary- Layer- Trailing-Edge 

(TBL-TE) Noise 

Using measured surface pressures, Brooks and 
Hodgson (ref. 2) demonstrated that if sufficient infor- 
mation is known about the TBL convecting surface 
pressure field passing the TE, then TBL-TE noise 
can be accurately predicted. Schlinker and Amiet 
(ref. 3) employed a generalized empirical description 
of surface pressure to predict measured noise. How- 
ever, the lack of agreement for many cases indicated 
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Turbulent-boundary-layer— trailing-edge 
noise 
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Laminar-boundary-layer — vortex-shedding 
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Trailing-edge-bluntness— vortex-shedding 
noise 



Figure 1. Flow conditions producing airfoil blade self- noise. 
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a need for a more accurate pressure description than 
was available. Langley supported a research effort 
(ref. 4) to model the turbulence within boundary lay- 
ers as a sum of discrete “hairpin” vortex elements. In 
a parallel and follow-up effort, the present authors 
matched measured and calculated mean boundary- 
layer characteristics to prescribed distributions of the 
discrete vortex elements so that associated surface 
pressure could be determined. The use of the model 
to predict TBL-TE noise proved disappointing be- 
cause of its inability to show correct trends with an- 
gle of attack or velocity. The results showed that to 
successfully describe the surface pressure, the history 
of the turbulence must be accounted for in addition 
to the mean TBL characteristics. This level of tur- 
bulence modeling has not been attempted to date. 

A simpler approach to the TBL TE noise prob- 
lem is based on the Ffowcs Williams and Hall (ref. 5) 
edge-scatter formulation. In reference 3, the noise 
data were normalized by employing the edge-scatter 
model with the mean TBL thickness 6 used as a 
required length scale. When 6 was unknown, sim- 
ple flat plate theory was used to estimate 6. Spec- 
tral data initially differing by 40 dB collapsed to 
within 7 dB, consistent with the results of the ap- 
proach discussed above using surface pressure mod- 
els. The extent of agreement between data sets was 
largely due to the correct scaling of the velocity de- 
pendence, which is the most sensitive parameter in 
the scaling approach. The dependence of the overall 
sound pressure level on velocity to the fifth power 
had been verified in a number of studies. The ex- 
tent to which the normalized data deviation was due 
to uncertainty in <5 was addressed by Brooks and 
Marcolini (ref. 6) in a forerunner to the present re- 
port. For large R c and small angles of attack, which 
matched the conditions of reference 3, the use of mea- 
sured TBL thicknesses <5, displacement thicknesses 
<5*, or momentum thicknesses 9 in the normalization 
produced the same degree of deviation within the 
TBL-TE noise data. Subsequently, normalizations 
based on boundary- layer maximum shear stress mea- 
surements and, alternately, profile shape factors were 
also examined. Of particular concern in reference 6 
was that when an array of model sizes, rather than 
just large models, was tested at various angles of at- 
tack, the normalized spectrum deviations increased 
to 10 or even 20 dB. These large deviations indicate 
a lack of fidelity of the spectrum normalization and 
any subsequent prediction methods based on curve 
fits. They also reinforce the conclusion from the 
aforementioned surface pressure modeling effort that 
knowledge of the mean TBL characteristics alone is 
insufficient to define the turbulence structure. The 
conditions under which the turbulence evolves were 


found to be important. The normalized data ap- 
peared to be directly influenced by factors such as 
Reynolds number and angle of attack, which in pre- 
vious analyses were assumed to be of pertinence only 
through their effect on TBL thickness S (refs. 3 and 

7 ). 

Several prediction schemes for TBL TE noise 
have been used previously for helicopter rotor noise 
(refs. 3 and 8) and for wind turbines (refs. 9 and 
10). These schemes have all evolved from scaling law 
equations which were fitted to the normalized data 
of reference 3 and, thus, are limited by the same con- 
cerns of generality discussed above. 

1.1.2. Separation- Stall Noise 

Assessments of the separated flow noise mecha- 
nism for airfoils at moderate to high angles of at- 
tack have been very limited (ref. 1). The relative 
importance of airfoil stall noise was illustrated in 
the data of Fink and Bailey (ref. 11) in an airframe 
noise study. At stall, noise increased by more than 
10 dB relative to TBL-TE noise, emitted at low an- 
gles of attack. Paterson et al. (ref. 12) found evidence 
through surface to far field cross-correlations that for 
mildly separated flow the dominant noise is emitted 
from the TE, whereas for deep stall the noise radiated 
from the chord as a whole. This finding is consistent 
with the conclusions of reference 11. 

No predictive methods are known to have been 
developed. A successful method would have to ac- 
count for the gradual introduction of separated flow 
noise as airfoil angle of attack is increased. Beyond 
limiting angles, deep stall noise would be the only 
major contributing source. 

1.1.3. Laminar- Boundary- Layer- Vortex- 

Shedding (LBL- VS) Noise 

When a LBL exists over most of at least one 
side of an airfoil, vortex shedding noise can oc- 
cur. The vortex shedding is apparently coupled 
to acoustically excited aerodynamic feedback loops 
(refs. 13, 14, and 15). In references 14 and 15, the 
feedback loop is taken between the airfoil TE and 
an upstream “source” point on the surface, where 
Tollmien-Schlichting instability waves originate in 
the LBL. The resulting noise spectrum is composed 
of quasi-tones related to the shedding rates at the 
TE. The gross trend of the frequency dependence 
was found by Paterson et al. (ref. 16) by scaling on a 
Strouhal number basis with the LBL thickness at the 
TE being the relevant length scale. Simple flat plate 
LBL theory was used to determine the boundary- 
layer thicknesses 6 in the frequency comparisons. 
The use of measured values of 6 in reference 6 veri- 
fied the general Strouhal dependence. Additionally, 
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for zero angle of attack, Brooks and Marcolini (ref. 6) 
found that overall levels of LBL VS noise could be 
normalized so that the transition from LBL VS noise 
to TBL TE noise is a unique function of R c . 

There have been no LBL VS noise prediction 
methods proposed, because most studies have em- 
phasized the examination of the rather erratic fre- 
quency dependence of the individual quasi-tones in 
attempts to explain the basic mechanism. However, 
the scaling successes described above in references 6 
and 16 can offer initial scaling guidance for the de- 
velopment of predictions in spite of the general com- 
plexity of the mechanism. 

LI. 4- Tip Vortex Formation Noise 

The tip noise source has been identified with the 
turbulence in the local separated flow associated with 
formation of the tip vortex (ref. 17). The flow over 
the blade tip consists of a vortex with a thick viscous 
turbulent core. The mechanism for noise production 
is taken to be TE noise due to the passage of the 
turbulence over the TE of the tip region. George and 
Chou (ref. 8) proposed a prediction model based on 
spectral data from delta wing studies (assumed to 
approximate the tip vortex flow of interest), mean 
flow studies of several tip shapes, and TE noise 
analysis. 

Brooks and Marcolini (ref. 18) conducted an ex- 
perimental study to isolate tip noise in a quantitative 
manner. The data were obtained by comparing sets 
of two- and three-dimensional test results for differ- 
ent model sizes, angles of attack, and tunnel flow ve- 
locities. From data scaling, a quantitative prediction 
method was proposed which had basic consistency 
with the method of reference 8. 

1.1.5 . Trailing- Edge- Bluntness- Vortex-Shedding 

Noise 

Noise due to vortex shedding from blunt trailing 
edges was established by Brooks and Hodgson (ref. 2) 
to be an important airfoil self-noise source. Other 
studies of bluntness effects, as reviewed by Blake 
(ref. 19) and Brooks and Schlinker (ref. 1), were only 
aerodynamic in scope and dealt with TE thicknesses 
that were large compared with the boundary- layer 


displacement thicknesses. For rotor blade and wing 
designs, the bluntness is likely to be small compared 
with boundary- layer thicknesses. 

Grosveld (ref. 9) used the data of reference 2 to 
obtain a scaling law for TE bluntness noise. He found 
that the scaling model could explain the spectral 
behavior of high-frequency broadband noise of wind 
turbines. Chou and George (ref. 20) followed suit 
with an alternative scaling of the data of reference 2 
to model the noise. For both modeling techniques 
neither the functional dependence of the noise on 
boundary-layer thickness (as compared with the TE 
bluntness) nor the specifics of the blunted TE shape 
were incorporated. A more general model is needed. 

1.2. Overview of Report 

The purpose of this report is to document the de- 
velopment of a self-noise prediction method and to 
verify its accuracy for a range of applications. The 
tests producing the data base for the scaling effort 
are described in section 2, In section 3, the mea- 
sured boundary- layer thickness and integral parame- 
ter data, used to normalize airfoil noise data, are doc- 
umented. The acoustic measurements are reported in 
section 4, where a special correlation editing proce- 
dure is used to extract clean self-noise spectra from 
data containing extraneous test rig noise. In sec- 
tion 5, the scaling laws are developed for the five self- 
noise mechanisms. For each, the data are first nor- 
malized by fundamental techniques and then exam- 
ined for dependences on parameters such as Reynolds 
number, Mach number, and geometry. The resulting 
prediction methods are delineated with specific calcu- 
lation procedures and results are compared with the 
original data base. The predictions are compared in 
section 6 with self-noise data from three studies re- 
ported in the literature. In appendix A, the data 
processing technique is detailed; in appendix B, the 
noise directivity functions are defined; and in appen- 
dix C, an application of the prediction methods is re- 
ported for a helicopter rotor broadband noise study. 
In appendix D, a computer code of the prediction 
method is given. 
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2. Description of Experiments 

The details of the measurements and test facil- 
ity have been reported in reference 6 for the sharp 
TE two-dimensional (2D) airfoil model tests, in ref- 
erence 18 for corresponding three-dimensional (3D) 
tests, and in reference 2 for the blunt TE 2D airfoil 
model test. Specific information applicable to this 
report is presented here. 

2.1. Models and Facility 

The models were tested in the low-turbulence po- 
tential core of a free jet located in an anechoic cham- 
ber. The jet was provided by a vertically mounted 
nozzle with a rectangular exit with dimensions of 
30.48 x 45.72 cm. The 2D sharp TE models are 
shown in figure 2. The models, all of 45.72-cm span, 
were NACA 0012 airfoils with chord lengths of 2.54, 
5.08, 10.16, 15.24, 22.86, and 30.48 cm. The models 
were made with very sharp TE, less than 0.05 mm 
thick, without beveling the edge. The slope of the 
surface near the uncusped TE corresponded to the 
required 7° off the chord line. The sharp TE 3D mod- 
els, shown in figure 3, all had spans of 30.48 cm and 
chord lengths that were the same as the five largest 
2D models. The 3D models had rounded tips, defined 
by rotating the NACA 0012 shape about the chord 
line at 30.48-cm span. An NACA 0012 model of per- 
tinence to the present paper, which is not shown here, 
is the blunt-TE airfoil of reference 2, with a chord 
length of 60.96 cm. Details of the blunt TE of this 
large model are given in section 5. 

The cylindrical hubs, shown attached to the mod- 
els, provided support and flush-mounting on the side 
plates of the test rig. At a geometric tunnel angle 
of attack at of 0°, the TE of all models was located 
61.0 cm above the nozzle exit. The tunnel angle at is 
referenced to the undisturbed tunnel streamline di- 
rection. In figure 4, an acoustic test configuration for 
a 3D model is shown. A 3D setup is shown so that 
the model can be seen fitted to the side plate. The 
side plates (152.4 x 30.0 x 1 cm) were reinforced 
and flush mounted on the nozzle lips. For the 2D 
configurations, an additional side plate was used. 

2.2. Instrumentation 

For all of the acoustic testing, eight 1.27-cm- 
diameter (1/2-in.) free-field- response microphones 
were mounted in the plane perpendicular to the 2D 
model midspan. One microphone was offset from 
this midspan plane. In figure 4, seven of these 
are shown with the identification numbers indicated. 
Microphones Ml and M2 were perpendicular to the 
chord line at the TE for at = 0°. The other 
microphones shown were at radii of 122 cm from 
the TE, as with Ml and M2, but were positioned 


30° forward (M4 and M7) and 30° aft (M5 and M8). 
The data acquisition and processing approaches are 
described in appendix A. 

For the aerodynamic tests the microphones to 
the right in figure 4 were removed and replaced by 
a large three-axis computer-controlled traverse rig 
used to position hot-wire probes. The miniature 
probes included both cross-wire and single-wire con- 
figurations. In figure 5, a cross-wire probe is shown 
mounted on the variable- angle arm of the traverse rig. 
Again, for clarity, a 3D airfoil model is shown. The 
probes were used to survey the flow fields about the 
models, especially in the boundary-layer and near- 
wake region just downstream of the trailing edge. 

2.3. Test Conditions 

The models were tested at free-stream velocities 
U up to 71.3 m/s, corresponding to Mach num- 
bers up to 0.208 and Reynolds numbers, based on a 
30.48-cm-chord model, up to 1.5 x 10 6 . The tunnel 
angles of attack at were 0°, 5.4°, 10.8°, 14.4°, 19.8°, 
and 25.2°. The larger angles were not attempted 
for the larger models to avoid large uncorrectable 
tunnel flow deflections. For the 22.86-cm- and 
30.48-cm-chord models, at was limited to 19.8° and 
14.4°, respectively. 

For the untripped BL cases (natural BL develop- 
ment), the surfaces were smooth and clean. For the 
tripped BL cases, BL transition was achieved by a 
random distribution of grit in strips from the lead- 
ing edge (LE) to 20 percent chord. This tripping 
is considered heavy because the chordwise extent of 
the strip produced thicker than normal BL thick- 
nesses. It was used to establish a well-developed 
TBL even for the smaller models and at the same 
time retain geometric similarity. The commercial 
grit number was No. 60 (nominal particle diam- 
eter of 0.29 mm) with an application density of 
about 380 particles/cm 2 . An exception was the 
2.54-cm-chord airfoil which had a strip at 30 per- 
cent chord of No. 100 grit with a density of about 
690 particles/cm 2 . 

2.4. Wind Tunnel Corrections 

The testing of airfoil models in a finite-size open 
wind tunnel causes flow curvature and downwash 
deflection of the incident flow that do not occur in 
free air. This effectively reduces the angle of attack, 
more so for the larger models. Brooks, Marcolini, and 
Pope (ref. 21) used lifting surface theory to develop 
the 2D open wind tunnel corrections to angle of 
attack and camber. Of interest here is a corrected 
angle of attack a* representing the angle in free air 
required to give the same lift as at would give in the 
tunnel. One has from reference 21, upon ignoring 
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Figure 2. Two-dimensional NACA 0012 airfoil blade models. 


Figure 3. Three-dimensional NACA 0012 airfoil blade models. 
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Figure 4. Test setup for acoustic tests of a 3D model airfoil. 
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Hot-wire probe 


Trailing edge 


3D airfoil model 


Figure 5. Tip survey using hot-wire probe. 
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small camber effects 


«* = «</ C (l) 

where 

C = (1 + 2a) 2 + y/l2a 

and 

a = (tt 2 /48 )(c/H) 2 


The term c is the airfoil chord and H is the tunnel 
height or vertical open jet dimension for a horizon- 
tally aligned airfoil. For the present 2D configura- 
tions, a,/a< equals 0.88, 0.78, 0.62, 0.50, 0.37, and 
0.28 for the models with chord lengths of 2.54, 5.08, 
10.16, 15.24, 22.86, and 30.48 cm, respectively. 
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3. Boundary-Layer Parameters at the 
Trailing Edge 

The purpose of this section is to present measured 
boundary- layer thicknesses from reference 21 and to 
document corresponding curve fit scaling equations 
to be employed in the normalization of the airfoil 
self-noise data. 

The data presented are the result of hot-wire 
probe measurements made in the boundary- layer/ 
near-wake region of the sharp TE of the 2D air- 
foil models. The probes were traversed perpendic- 
ular to the model chord lines downstream of the TE. 
These measurements were made at 0.64 mm from 
the TE for the 2.54-cm-chord airfoil and at 1.3 mm 
for the other airfoils. The integral BL parameters — 
displacement thickness 6* and momentum thickness 
6 — were calculated from mean velocity profiles with 
the BL/near-wake thickness 6 specified. The thick- 
ness 6 is that distance from the airfoil surface where 
the mean velocity reaches 99 percent of the potential 
flow stream velocity. The values of 6 were chosen by 
carefully examining the respective turbulent veloc- 
ity and Reynolds stress distributions as well as the 
mean profiles. For all cases, the estimated accuracy 
of 6 is within ±5 percent for the turbulent-boundary- 
layer (TBL) flow and ±10 percent for the laminar 
and transitional flows, whereas the error range for 
the integral thicknesses 6* and 9 is less (ref. 21). 

3.1. Scaled Data 

The thicknesses 6 and integral properties 6* and 9 
at the TE of the sharp TE 2D airfoil models at at = 
0° are given in figure 6 for both the artificially tripped 
and the untripped boundary- layer conditions. The 
subscript 0 for the thicknesses indicates that the 
airfoil is at zero angle of attack. The parameters 
are normalized by the chord length c and are given 

3.2. Calculation Procedures 


as a function of Reynolds number based on the 
chord R c . As R c increases, the thicknesses decrease 
for both the tripped and the untripped boundary 
layers. The tripped boundary layers are almost 
uniformly thicker than the corresponding untripped 
boundary layers. One should refer to reference 21 for 
details of the boundary- layer character for the cases 
of figure 6. In general, however, one can say that the 
tripped boundary layers are fully turbulent for even 
the lowest R c . The untripped boundary layers are 
laminar or transitional at low R c and become fully 
turbulent for high R c . In figure 6, the boundary-layer 
thickness data are approximated by curve fits whose 
equations are specified in the following section. 

Angle-of-attack effects on the thickness parame- 
ters are given at free-stream velocities of 71.3 and 
39.6 m/s for the untripped and tripped BL airfoils in 
figures 7 and 8. The parameters are normalized by 
those measured for the corresponding cases at zero 
angle of attack, given in figure 6. The data are plot- 
ted against the corrected angle a * of equation (1). 
The collapse of the data is much improved over that 
when at is used (ref. 21). In general, the data show 
that for increasing a* (or at) the thicknesses increase 
on the suction side because of the increasing sever- 
ity of the adverse pressure gradient. The converse is 
true for the pressure side, where the pressure gradi- 
ent becomes increasingly favorable. Also included in' 
figures 7 and 8 are curve fits to the data. For the 
pressure side of the airfoils, the curves are the same 
for the tripped and untripped cases. The suction side 
curves differ, reflecting differences in the angle depen- 
dence of where the TE boundary layer separates and 
finally stalls the airfoil. 

In reference 21, the data are discussed and com- 
pared with flat plate experimental results and results 
from boundary-layer prediction codes. 


The boundary-layer thickness parameters at the TE of a symmetric NACA 0012 airfoil at zero angle of 
attack are approximated by the curve fits to the data of figure 6. The expressions for the curve fits for 
boundary- layer thickness <5, displacement thickness <5*, and momentum thickness 9 are, for the heavily tripped 
boundary layer, 

6q/c= lot 1 892-0 9045 lo S *r+0.0596(log /^) 2 ] ( 2 ) 


0.0601 i?r°- 114 


6Z/ C = 


(R c < 0.3 x 10 6 ) 


10 [3.411— 1.5397 log i* c +0.1059(log R c ) 2 ] ^ > q 3 x 1 q 6 ) 


9q/c = 


0.0723 fi" 0 - 1765 


(R c < 0.3 x 10 6 ) 


10 [0. 5578-0. 7079 log -I- 0.0404 (log #D 2 ] > 0 3 X 10^) 


(3) 

(4) 
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where the zero subscripts indicate zero angle of attack, zero lift on these symmetric airfoils. For the untripped 
(natural transition) boundary layers , 

<5 0 /c = lot 1 6569-0 9045 lo 8 flc+0-0596(log iic) 2 ] ( 5 ) 


S*/c= lot 30187-1 5397 lo 8 -Rc+0.1059(log Re) 2 ] 


( 6 ) 


Qq/c = lot 0 - 2021- °- 7079 lo S +0.0404 (log Re) 2 ] 


( 7 ) 


The boundary- layer thicknesses for the airfoils at nonzero angle of attack, in terms of the zero-angle-of- 
attack thicknesses and the corrected angles a*, are given in figures 7 and 8. The expressions for the curve fits 
for the pressure side, for both the tripped and the untripped boundary layers , are 


(jP _ Jq[ — 0. 04175a, +0.00106aJ] 

, n [— 0.0432a, +0.001 13a?] 

Jj* 1U 
°0 

_ 10 [-0.04508a,+0.000873a2] 

#0 


( 8 ) 

( 9 ) 

( 10 ) 


For the suction side, the parametric behavior of the thicknesses depends on whether the boundary layers are 
attached, separated near the trailing edge, or separated a sufficient distance upstream to produce stall. For 
the suction side for the tripped boundary layers (fig. 7), 


1q 0. 0311a, 

(0° < a* < 5°) 


0.3468(10 0 1231q *) 

(5° < a* < 12.5°) 

(11) 

5.7 1 8(10 0 0258a *) 

(12.5° < a, < 25°) 


IqO. 0679a* 

(0° < a. < 5°) 


0.381(10° 1516a ») 

(5° < a, < 12.5°) 

(12) 

14.296(10 0 0258a *) 

(12.5° < a* < 25°) 


lQ0.0559a+ 

(0° < a, < 5°) 


0.6984(10°' O869q *) 

(5° <«,< 12.5°) 

(13) 

4.0846(10°° 258a *) 

(12.5° < a* < 25°) 
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For the suction side for the untripped boundary layers (fig. 8), 


6 * 

*0 


r 10 0.03114 q. ( 0 o < am < 7 5 0 ) 

i 0.0303(10°- 2336q *) (7.5° < a* < 12.5°) 

k 12(10°° 258a *) (12.5° < a, < 25°) 


6 

6 


* 

S_ 

* 

0 


6 jl 

Oq 


{ 10 0.0679a, (fio < a* < 7.5°) 

0.0162(10°- 3066q -) (7.5° < a* < 12.5°) 

52.42(10 00258q *) (12.5° < a* < 25°) 

' 10 0 0559a * ( 0 ° < a* < 7 . 5 °) 

< 0.0633(10°- 2157a *) (7.5° <a ( < 12.5°) 

k 14.977(10 0 0258q *) (12.5° < a* < 25°) 


(14) 


(15) 


( 16 ) 



4. Acoustic Measurements 

The aim of the acoustic measurements was to de- 
termine spectra for self- noise from airfoils encoun- 
tering smooth airflow. This task is complicated by 
the unavoidable presence of extraneous tunnel test 
rig noise. In this section, cross-correlations between 
microphones are examined to identify the self-noise 
emitted from the TE in the presence of other sources. 
Then, the spectra of self-noise are determined by per- 
forming Fourier transforms of cross-correlation data 
which have been processed and edited to eliminate 
the extraneous contributions. The results are pre- 
sented as 1 /3-octave spectra, which then form the 
data base from which the self-noise scaling predic- 
tion equations are developed. 

4.1. Source Identification 

The upper curves in figure 9 are the cross- 
correlations, R\ 2 {t) = (p\{t)p2(t + t)), between the 
sound pressure signals p\ and p 2 of microphones Ml 
and M2 identified in figure 4. Presented are cross- 
correlations both with and without the tripped 30.48- 
cm-chord airfoil mounted in the test rig. Because the 
microphones were on opposite sides of, and at equal 
distance from, the airfoil, a negative correlation peak 
occurs at a signal delay time r of 0. This correlation 
is consistent with a broadband noise source of dipole 
character, whose phase is reversed on opposing sides. 
When the airfoil is removed, the strong negative peak 
disappears leaving the contribution from the test rig 
alone. The most coherent parts of this noise are from 
the lips of the nozzle and are, as with the airfoil noise, 
of a dipole character. The microphone time delays 
predicted for these sources are indicated by arrows. 
The predictions account for the effect of refraction of 
sound by the free-jet shear layer (refs. 22 and 23), as 
well as the geometric relationship between the micro- 
phones and the hardware and the speed of sound. 

The lower curves in figure 9 are the cross- 
correlations, #45 (r), between microphones M4 and 
M5. The predicted delay times again appear to cor- 
rectly identify the correlation peaks associated with 
the noise emission locations. The peaks are positive 
for i?45(r) because both microphones are on the same 
side of the dipoles' directional lobes. The noise field 
is dominated by TE noise. Any contribution to the 
noise field from the LE would appear where indicated 
in the figure. As is subsequently shown, there are 
contributions in many cases. For such cases the neg- 
ative correlation peak for R\ 2 (t) would be the sum 
of the TE and LE correlation peaks brought together 
at r = 0 and inverted in sign. 

In figure 10, the cross-correlations #45 (r) are 
shown for tripped BL airfoils of various sizes. The 
TE noise correlation peaks are at txe = -0.11 ms 


for all cases because at at = 0°, the TE location 
of all models is the same. The LE location changes 
with chord size, as is indicated by the change in the 
predicted LE noise correlation peak delay times. 

For the larger airfoils in figure 10, the TE con- 
tribution dominates the noise field. As the chord 
length decreases, the LE noise peaks increase to be- 
come readily identifiable in the correlation. For the 
smallest chord the LE contribution is even some- 
what more than that of the TE. Note the extraneous, 
but inconsequential, source of discrete low-frequency 
noise contributing to the 22.86-cm-chord correlation, 
which can be readily edited in a spectral format. 

It is shown in reference 6 that the LE and TE 
sources are uncorrelated. The origin of LE noise 
appears to be inflow turbulence to the LE from the 
TBL of the test rig side plates. This should be the 
case even though the spanwise extent of this TBL 
is small compared with the portion of the models 
that encounter uniform low- turbulence flow from the 
nozzle. Inflow turbulence can be a very efficient 
noise mechanism (ref. 24); however its full efficiency 
can be obtained only when the LE of the model 
is relatively sharp compared with the scale of the 
turbulence. The LE noise contributions diminish for 
the large chord because of the proportional increase 
in LE radius with chord. When this radius increases 
to a size that is large compared with the turbulent 
scale in the side plate TBL, then the sectional lift 
fluctuations associated with inflow turbulence noise 
are not developed. 

4.2. Correlation Editing and Spectral 

Determination 

The cross-spectrum between microphones Ml and 
M2, denoted G\ 2 {f), is the Fourier transform of 
J?12 (r). If the contributions from the LE, nozzle lips, 
and any other coherent extraneous source locations 
were removed, G 12 (f) would equal the autospectrum 
of the airfoil TE self-noise, S(f). Actually the rela- 
tionship would be G\ 2 (f) = S(f) exp[i(27r fv $ e ±7r)], 
where i = y/—l and I s the delay time of the TE 
correlation peak. This approach is formalized in ref- 
erence 2. 

In reference 6, the spectra were found from Gi2(/) 
determined with the models of the test rig after a 
point-by-point vectorial subtraction of G 12 (f) deter- 
mined with the airfoil removed. This was equivalent 
to subtracting corresponding i?i2( T ) results, such as 
those of figure 9, and then taking the Fourier trans- 
form. This resulted in “corrected" spectra which 
were devoid of at least a portion of the background 
test rig noise, primarily emitted from the nozzle lips. 
The spectra still were contaminated by the LE noise 
due to the inflow turbulence. 
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Figure 9. Cross-correlations for two microphone pairs with and without airfoil mounted in test rig. 
c = 30.48 cm; BL tripped; a t = 0°; V = 71.3 m/s. Arrows indicate predicted values of r. (From ref. 6.) 



Figure 10. Cross-correlations between microphones M4 and M5 for tripped BL airfoils of various chord sizes. U = 71.3 m/s. 
Arrows indicate predicted values of r. (From ref. 6.) 
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In the present paper, most spectra presented 
were obtained by taking the Fourier transform of 
microphone-pair cross-correlations which had been 
edited to eliminate LE noise (see details in appen- 
dix A). The microphone pairs used included M4 and 
M5, M4 and M8, and M4 and M2. These pairs pro- 
duced correlations where the TE and LE noise peaks 
were generally separated and readily identifiable. Re- 
ferring to figure 10 for i? 45 (r), the approach was to 
employ only the left-hand side (LHS) of the TE noise 
peak. The LHS was “folded” about r at the peak 
(txe) to produce a nearly symmetrical correlation. 
Care was taken in the processing to maintain the ac- 
tual shapes near the very peak, to avoid to the extent 
possible the artificial introduction of high-frequency 
noise in the resulting spectra. Cross-spectra were 
then determined which were equated to the spectra 
of TE self-noise. 

The data processing was straightforward for the 
larger chord airfoils because the LE and TE peaks 
were sufficiently separated from one another that the 
influence of the LE did not significantly impact the 
TE noise correlation shapes. For many of the smaller 
airfoils, such as those with chord lengths of 2.54, 5.08, 
and 10.16 cm shown in figure 10, the closeness of the 
LE contribution distorted the TE noise correlation. 
A processing procedure was developed to effectively 
“separate” the TE and LE peaks to a sufficient dis- 
tance from one another, within the correlation pre- 
sentation, so that the correlation folding of the LHS 
about txe produced a more accurate presentation of 
the TE noise correlation shape. The separation pro- 
cessing employed symmetry assumptions for the TE 
and LE noise correlations to allow manipulation of 
the correlation records. This processing represented 
a contamination removal method used for about one- 


quarter of the spectra presented for the three small- 
est airfoil chord lengths. Each case was treated in- 
dividually to determine whether correlation folding 
alone, folding after the separation processing, or not 
folding at all produced spectra containing the least 
apparent error. In appendix A, details of the edit- 
ing and Fourier transform procedures, as well as the 
separation processing, are given. 

4.3. Self-Noise Spectra 

The self-noise spectra for the 2D NACA 0012 
airfoil models with sharp TE are presented in a 
1 /3-octave format in figures 11 to 74. Figures 11 
to 43 are for airfoils where the boundary layers have 
been tripped and figures 44 to 74 are for smooth sur- 
face airfoils where the boundary layers are untripped 
(natural transition). Each figure contains spectra for 
a model at a specific angle of attack for various tun- 
nel speeds. Note that the spectra are truncated at 
upper and lower frequencies. This editing of the spec- 
tra was done because, as described in appendix A, a 
review of the narrow-band amplitude and phase for 
all cases revealed regions where extraneous noise af- 
fected the spectra in a significant way (2 dB or more). 
These regions were removed from the 1 /3-octave 
presentations. 

The spectra levels have been corrected for shear 
layer diffraction and TE noise directivity effects, as- 
detailed in appendix B. The noise should be that for 
an observer positioned perpendicular to, and 1.22 m 
from, the TE and the model midspan. In terms of the 
directivity definitions of appendix B, r € — 1.22 m, 
0 e = 90°, and 4> e — 90°. In section 5 (beginning 
on p. 51), the character and parametric behavior of 
the self-noise, as well as the predictions which are 
compared with the data, are discussed. 
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Figure 13. Self-noise spectra for 30.48-cm-chord airfoil with tripped BL at a* = 10.8° (a, = 
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Figure 14. Self-noise spectra for 30.48-cm-chord airfoil with tripped BL at at = 14.4° (a* = 




















































































1 

10 

20 .2 

1 

Frequency, kHz 
3 ) U = 71.3 m/s 



Frequency, kHz 
(b) U = 55.5 m/s 




Frequency, kHz 
(c) U = 39.6 m/s 


Frequency, kHz 
(d) U = 31.7 m/s 


Figure 27. Self-noise spectra for 10.16-cm-chord airfoil with tripped BL at at = 5.4° (a* = 3 
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Figure 28. Self-noise spectra for 10.16-cm-chord airfoil with tripped BL at a* = 10.8° (a* = 6.7°). 



Figure 29. Self-noise spectra for 10.16-cm-chord airfoil with tripped BL at a t = 14.4° (a* = 8.9°). 
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Figure 36. Self-noise spectra for 5.08-cm-chord airfoil with tripped BL at at — 19.8° (a* = 15.4°). 



Figure 37. Self-noise spectra for 5.08-cm-chord airfoil with tripped BL at at = 25.2° (a* = 19.7°). 
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Figure 44. Self- noise spectra for 30.48-cm-chord airfoil with untripped BL (natural transition) at 
(a, = 0°). 
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Figure 45. Self-noise spectra for 30.48-cm-chord airfoil with untripped BL at a t = 5.4° (a* = 
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Figure 55. Self-noise spectra for 15.24-cm-chord airfoil with untripped BL at a t = 10.8° (a* = 5.4°). 



(a) U = 71.3 m/s (b) U = 39.6 m/s 

Figure 56. Self- noise spectra for 15.24-cm-chord airfoil with untripped BL at a t = 14.4° (a* = 7.2°). 
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Figure 59. Self-noise spectra for 10.16-cm-chord airfoil with untripped BL at a* = 0° (a* = 
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Figure 60. Self-noise spectra for 10.16-cm-chord airfoil with untripped BL at a t = 5.4° (a* = 
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Figure 63. Self-noise spectra for 10.16-cm-chord airfoil with untripped BL at a t = 19.8° (a* = 12.3°). 



(a) U = 71.3 m/s (b) U = 39.6 m/s 

Figure 64. Self-noise spectra for 10.16-cm-chord airfoil with untripped BL at at = 25.2° (a* = 15.6°). 
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5. Spectral Scaling 

In this section, the scaling laws are developed 
for the five self-noise mechanisms. The spectra of 
figures 11 to 74 form the basis of the scaling for 
three of the mechanisms: turbulent-boundary-layer- 
trailing-edge (TBL-TE) noise and separation noise 
were scaled from the tripped boundary- layer cases, 
and laminar-boundary-layer-vortex-shedding (LBL- 
VS) noise was scaled from the untripped cases. For 
the tip vortex formation noise mechanism, both the 
data and the scaling approach are obtained from ref- 
erence 18. Finally, for TE-bluntness vortex-shedding 
noise, spectral data from the study of reference 2, as 
well as previously unpublished data from that study, 
form the basis of scaling analysis. 

5.1 . Turbulent-Boundary-Layer-Trailing-Edge 
Noise and Separated Flow Noise 

What has become traditional TE noise scaling is 
based on the analysis of Ffowcs Williams and Hall 
(ref. 5). For the problem of turbulence convecting at 
low subsonic velocity U c above a large plate and past 
the trailing edge into the wake, the primary result is 


As mentioned in section 1, some of the airfoil self- 
noise spectral data of the present report were pre- 
sented, in uncorrected form, in reference 6, and nor- 
malized in the manner of equation (18) using mea- 
sured values of 5. It was found that, contrary to what 
was previously assumed (e.g., refs. 25 and 3), the nor- 
malized levels, spectral shape, and Strouhal number 
were not independent of airfoil size, airfoil angle of 
attack, and free-stream velocity. However, the lim- 
ited scope of the paper, as well as the uncertainty 
caused by the aforementioned extraneous noise con- 
tamination of the uncorrected spectra, prevented a 
clear definition of the functional dependences. The 
corrected spectra of the present report are used to de- 
termine the parametric dependences and to account 
for these in the spectral scaling. 

5.1.1. Scaled Data 

Zero angle of attack. In figure 75, 1/3-octave 
spectra for four airfoil sizes, each at four tunnel 
speeds, are scaled. The spectra are obtained from 
figures 11, 20, 26, and 32. The angle of attack is zero 
and the boundary layers are tripped. The form of 
the normalization is 


c 0 



(17) 


Scaled SPLj^ = SPL^ — 



(19) 


where (p 2 ) is the mean-square sound pressure at the 
observer located a distance r from the edge. The 
medium density is pq , t/ 2 is the mean-square tur- 
bulence velocity, cq is the speed of sound, L is the 
spanwise extent wetted by the flow, and £ is a charac- 
teristic turbulence correlation scale. The directivity 
factor D equals 1 for observers normal to the surface 
from the TE. The usual assumptions for boundary- 
layer flow are that v f a U c oc U and £ oc 6 or 6*, 
where 6 and 6* are, respectively, the boundary-layer 
thickness and displacement thickness. Fink (ref. 25), 
when normalizing airframe noise data where TBL- 
TE noise was believed to be dominant, assumed 
a universal spectrum shape F( St) for the noise, 
where St is the Strouhal number f6/U. The shape 
F(St) depended only on the ratio of St to its peak 
value Stp ea k. This gave the following normalization 
form for the 1 /3-octave sound pressure level spectral 
presentation: 


SPL 1/3 — 10 log 



= F{St) + K 


(18) 


with SPL ^3 = OASPL + F(St) and where K is an 
empirical constant which was determined when the 
velocity U is given in units of knots. 


where Mach number replaces the velocity in knots, 6q 
replaces <5, and r e replaces r. The retarded observer 
distance r € equals here the measured value, 122 cm 
(see appendix B). For the right side of equation (19) 
to be accurately expressible by the form F(St) + K of 
equation (18), the scaled spectra of figure 75 should 
be identical to one another for all cases. However, 
the peak Strouhal number, spectral shape, and scaled 
level vary significantly. 

For each spectrum in figure 75, a symbol indi- 
cates the approximate spectral peak location. The 
peak locations were based on gross spectral shapes 
and trends rather than specific peak maximums. The 
peak Strouhal number, St pea k = (/6*/[/) pea k, an d 
scaled levels corresponding to these peak locations 
are shown in figures 76 and 77, respectively, as a 
function of Reynolds number R c . These data are 
also presented in table 1 (at the back of this re- 
port). Included in the figures are the other cases for 
tripped BL airfoils of different chord lengths. Also 
included are data at nonzero angle of attack for sub- 
sequent discussion. The displacement thicknesses for 
the suction side, 6*, are used for these normaliza- 
tions. In figure 76, St pea k for zero angle of attack 
(solid symbols) shows no clear i? c -dependence, but a 
Mach number dependence is apparent. The horizon- 
tal lines through the data correspond to the function 
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Figure 75. Scaled 1/3-octave spectra of tripped BL airfoils at at — 0° (a* = 0°). Symbols indicate approximate spectral peak locations 
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Figure 76. Peak Strouhal number for TBL-TE noise versus Reynolds number. Numbers aligned with data are chord sizes in 
inches (for brevity of notation). 
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Figure 77. Peak scaled level for TBL TE noise versus Reynolds number. Numbers aligned with data are chord sizes in inches, 






Sti = 0.02M~ 0,6 for the presented values of Mach 
number and is taken to approximate the behavior of 
St peak- For the scaled levels in figure 77, a continu- 
ous function, designated as K\ } that is comprised of 
i2 c -dependent segmented lines is drawn to approx- 
imate the zero-angle-of-attack data. Other choices 
for a function to approximate these data are possible 
but the one shown, which is chosen to be constant 
for high i? c , was found to be compatible with higher 
Reynolds number data obtained from other studies, 
as is shown subsequently. Note that the behavior of 
K i at very low R c is at most academic because of the 
lack of importance of this TBL-TE noise mechanism 
in this range. 

In figure 78(a), a shape function denoted by A is 
proposed as representative of the 1/3-octave spectral 
shape of the TBL TE noise mechanism. (Fig. 78(b) 
presents a corresponding shape function for sepa- 
rated flow noise.) The spectrum A is a function of the 
ratio St/St pea k that is symmetric about St/St pea k = 
1.0. The spectral width or broadness depends on Rc. 
Two extremes in A are shown corresponding to so- 
called maximum and minimum Reynolds numbers. 
Intermediate values of R c require interpolation. As 
seen in figure 75, the larger chords have the broad- 
est TBL-TE spectra. The spectrum A was matched 
to these and the other chord lengths. The specific 
details of A and the other functions are given in the 
calculation procedures section (5.1.2.). 

One of the key results of reference 2 is that 
each side of an airfoil with well-developed boundary 
layers produces TBL-TE noise independently of the 
other side. This is not in conflict with our scaling 
approach for the symmetric airfoil at zero angle of 
attack. Consistency of this with equation (19) merely 
requires a level adjustment (—3 dB) of the scaling 
equations to account for the equal contributions of 
the two sides to the total spectrum. For the pressure 
and suction sides, i = p or s, 

rf ) 

= A (If) + <*> - 3 > < 20 > 

where St* = (f6*/U). The total TBL-TE noise for 
zero angle of attack then is 

SPLtbl-TE = 10 log (10 SPL ‘/1° + io spl p/ 10 ) (21) 

where a 1 /3-octave presentation for spectra is 
understood. 

Nonzero angle of attack. In figure 79, scaled noise 
spectra are presented for the same tripped BL airfoil 


models as in figure 75, but here the angle of attack 
is varied while holding tunnel velocity constant at 
U = 71.3 m/s. The tunnel angles of attack at are 
given along with the effective angles a*. The level 
normalization approach and Strouhal scaling are the 
same as in figure 75 except that here the displacement 
thickness of the suction side of the airfoil 8* is used. 
For increasing a * the peak Strouhal number and 
level increase and the spectra become sharper at 
the peaks. Beyond limiting values of a *, roughly 
corresponding to stall, substantial changes occur to 
the scaled spectra. 

If equations (20) and (21) were used to predict 
the spectra in figure 79 and the predictions scaled 
accordingly, one would find for increasing angle of at- 
tack that peak Strouhal number would remain con- 
stant, peak level would decrease, and the spectral 
shape would become broader at the peak. This is 
because the suction side contribution would remain 
dominant and that of the pressure side would shift to 
higher frequencies at reduced levels. These trends, 
of course, are virtually opposite to those observed. 
The approach that is now taken is to postulate at 
nonzero angles of attack an additional contribution 
to the spectrum that controls the spectral peak. To 
justify this, one could hypothesize that the spectrum 
is the total from attached TBL contributions, as for- 
mulated in equations (20) and (21), and a contri- 
bution from a separated portion of the TBL on the 
suction side. The modeling approach, however, is not 
without conflict at the low Reynolds numbers, as is 
discussed subsequently. Model details are developed 
below, after establishing the Strouhal and level scal- 
ing behavior for the angle cases. 

In figure 79, for each spectrum, symbols indicate 
the approximate peak Strouhal locations. As in fig- 
ure 75, the locations of the peaks were based on 
gross trends and shapes of the spectra rather than 
precise peaks. These values of St peak are included 
in figure 76 for the various chords, speeds, and an- 
gles of attack, along with the zero angle values pre- 
viously discussed. Again little direct i? c -dependence 
is noted for Stp^. The basic trends observed can be 
explained by velocity and angle dependence. The val- 
ues of St pea k are plotted versus corrected angle of at- 
tack a* in figure 80. For reference, the chord lengths 
(in units of inches for presentation convenience) are 
given. Through the data are drawn data-fit lines des- 
ignated as St 2 , corresponding to two velocity values. 
At a* =0°, St 2 becomes the function St i of figure 76. 
In the hand-fitting procedure to determine St 2 , some 
preference was given to the higher speed cases. This 
preference is discussed subsequently with regard to 
Strouhal peak level scaling. As for the substantial 


Scaled SPL* = SPL^ — 10 log ( M l 


54 



Function A 
level, dB 


Function B 
level, dB 



(a) Function A for TBL-TE noise, equations (35) to (40). 



(b) Function B for separated flow noise, equations (41) to (46). 

Figure 78. One-third-octave spectral shapes as functions of Strouhal and Reynolds numbers. 
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data scatter of figure 80, some comments are war- 
ranted. It was found that if one used the actual mea- 
sured values of <5, (where available) in the Strouhal 
scaling, one would have a similar degree of scatter 
to that shown in figure 80, where scaled values of 
6* (eq. (12)) were used. Also if untripped BL air- 
foil results were plotted, for those limited number of 
cases where the LBL-VS shedding source is not ap- 
parent in the spectra, the scatter and trend would be 
about the same as those shown in figure 80. Other 
deviations of the data from the St2 lines occur at 
mid to high angles of attack, where the low-frequency 
parts of the spectra were limited by the experimental 
high-pass filtering and thus values of St p( . a ^ were in- 
accurately large. The behavior of St2 seen in figure 80 
at the higher angles of attack (where the horizontal 
lines are placed lower than the data) was chosen to 
approximately correct this bias. 

The scaled levels corresponding to spectral peaks 
chosen in figure 79 are shown in figure 77 with 
the other cases. The previously indicated conflict 
within the data base for the proposed modeling ap- 
proach, which hypothesizes contributions from two 
attached TBL’s and an angle-dependent separation- 
related portion, is seen in figure 77. Peak levels for 
the two smallest chord lengths, except at the highest 
speeds, significantly decrease as the angle of attack 
increases from zero. This is incompatible with the 
modeling approach. A choice is made to ignore the 
conflicting low Reynolds number data in the model 
development. While admitting that the inclusion of 
the low Reynolds number behavior would conceptu- 
ally be desirable for completeness of the modeling, 
the exclusion is believed justifiable because of the 
greater interest in higher Reynolds number condi- 
tions. The TBL TE noise mechanism is not con- 
sidered important for low Reynolds numbers. Even 
if this were not the case, it is not certain that the 
present test flow conditions with heavy leading-edge 
tripping for airfoils at nonzero angles of attack prop- 
erly represent the mechanism, especially for higher 
angles where relaminarization of the pressure-side 
boundary layer is possible. Regardless, the results 
of the scaling are compared subsequently with the 
spectra of all the data to allow a direct assessment of 
the effect of modeling choices. 

The scaled levels of figure 77 for chord lengths 
of 10.16, 15.24, 22.86, and 30.48 cm are plotted 
in figure 81 versus a*. If the portion of these 
levels that cannot be accounted for by the modeling 
of equations (20) and (21) can be extracted, this 
portion would be designated as the separated flow 
noise contribution. Calculations were performed by 
taking into account that the Strouhal dependence 
of A in equation (20) would follow Sti of figure 76 


rather than St 2 of figure 80, which applies to that 
portion extracted. The extracted levels are given 
in figure 82. These extracted levels are normalized 
by subtracting the zero-angle-of-attack function of 
figure 77 (K\) for the particular chord lengths and 
speeds. Although substantial scatter is present, a 
basic trend of increasing importance for increasing 
angle and speed is seen. Drawn through the data is 
a function designated as K 2 — A’l which represents 
a partially observed, partially postulated dependence 
on velocity and angle of attack. The assigned spectral 
shape for this additive source is function B, which 
is given in figure 78(b) and is defined in a manner 
similar to function A of figure 78(a) to have a width 
which is dependent on chord Reynolds number. 

The resulting scaling model for the angle- 
dependent noise SPL a is 



= B (fe ) +K2 <22) 

where this represents the separated-boundary-layer 
noise contribution to the total noise. The total TBL 
TE and separation noise is then 

SPLtot = 10 log (l0 SPLa/1 ° + 10 SPL '/ 10 

4. 10 SPI V1°) (23) 

During development of the scaling procedures, 
equations (20), (22), and (23) were compared with 
spectra for all tripped BL airfoils and with spectra 
for the untripped BL airfoils for which TBL TE noise 
appeared to significantly contribute. Analyses of 
comparisons resulted in optimization of curves A and 
B, as well as development of the specific calculation 
procedures. The analysis found that better results 
are obtained when the Strouhal dependency of the 
suction-side spectrum SPL s is (St 1 + St2 )/2 rather 
than Sti. It was found that for better SPL agree- 
ment, one should make an adjustment in pressure- 
side level SPLp (defined as AKi in the following sec- 
tion) as a function of angle of attack and Reynolds 
number based on the displacement thickness 6*. This 
adjustment diminishes the pressure-side contribution 
for increasing angle and decreasing velocity. Also it 
was found that the drastic spectral shape changes 
that occur at sufficiently high angles of attack, near 
stall, are roughly simulated by a calculation proce- 
dure change. At the value of a* corresponding to 
the peak of the appropriate K 2 curve, the spectral 
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contributions SPL 5 and SPLp in equation (23) are 
eliminated and the B curve of equation (22) is re- 
placed by an A curve corresponding to a value of R c 
which is three times the actual value. 


The calculation procedures are specified in the 
next section followed by comparison with the spectral 
data base. 



Figure 82. Angle-dependent scaled noise levels as referenced to zero angle of attack, TBL-TE noise model. Numbers aligned 
with data are chord sizes in inches. 


5. 1.2. Calculation Procedures 

The total TBL-TE and separation noise spectrum in a 1 /3-octave presentation is predicted by 

SPLjot = 10 log (io SPL «/ 10 + io sp W 10 + io spl p/ 10 ) 

where 


SPLp = 10 log 
SPL* = 10 log 


' 6*M rj LD h ' 


' 6*M 5 LD h 


+ ^(l7) +(K '- 3) + AK » 

+ A (fir) + {Kl ~ 3) 


and 


SPL a = 10 log 


' 6*M b LD 


‘)**(» 


+ k 2 


for angles of attack up to (a*)o, an angle to be defined later in this section. At angles above (a*)o, 


SPLp = — oc 


SPL s = — oo 


(24) 


(25) 

(26) 


(27) 


(28) 

(29) 
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and 


SPL a = 10 log 


6* g M 5 LD e \ 

rl ) 



(30) 


where A! is the curve A but for a value of R c which is three times the actual value. The directivity functions 
Dfr and D £ are given in appendix B by equations (Bl) and (B2), respectively. 

The Strouhal definitions are (see figs. 76 and 80) 


f6* 

St p = IjB- St 3 

* W3|. 

II 

(31) 

Sti = 0.02M -0 ' 6 


(32) 

oT Sti + St2 

St l= 2 


(33) 

1 

(a* < 1.33°) 


10 0.0054(q,-1.33) 2 

(1.33° < a* < 12.5°) 

(34) 

4.72 

(12.5° < a*) 



For the spectral shape function definitions, we first consider the function A of figure 78(a). As discussed, 
the function A for a particular Reynolds number R c is obtained from an interpolation of the curves A nrdx and 
A m in, corresponding to chosen values, (R c ) max and {R c )mm- The two curves are defined as 


^min( a ) 


r v/67.552 - 886.788a 2 - 8.219 (a < 0.204) 

< -32.665a + 3.981 (0.204 < a < 0.244) 

k -142.795a 3 + 103.656a 2 - 57.757a 4- 6.006 (0.244 < a) 


(35) 


and 


( ^67.552 - 886.788a 2 - 8.219 


(a < 0.13) 


^4max( a ) — \ 


-15.901a + 1.098 

-4.669a 3 + 3.491a 2 - 16.699a + 1.149 


(0.13 < a < 0.321) 
(0.321 < a) 


(36) 


where a is the absolute value of the logarithm of the ratio of Strouhal number, St = St p or St 5 , to the peak 
Strouhal number, St pea ]< = Sti, Sti, or St 2 : 


a = I log( St /Stpeak ) I (37) 

The absolute value is used because the spectral shape is modeled to be symmetric about a = 0. 

The interpolative procedure includes defining a value, ao(Rc), at which the spectrum has a value of —20 dB. 
This -20 dB corresponds to a horizontal axis intercept in figure 78(a) for an interpolated curve. The function 
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ao(R c ) is given by 


0.57 (Rc < 9.52 x 10 4 ) 

a 0 (fl c ) = { (-9.57 x 10- 13 )(fl c - 8.57 x 10 5 ) 2 + 1.13 (9.52 x 10 4 < R c < 8.57 x 10 5 ) 

1.13 (8.57 x 10 5 < R c ) 

An interpolation factor A^(ao) is determined from 

— 20 — A m i n (ao) 


(38) 


^fl( a o) = 


A max ( a o) ^4min( a o) 


(39) 


where A m j n (ao) and A m ax(«o) are the A max and A m \ n spectra evaluated at a$. The spectrum shape A can now 
be evaluated for any frequency by computing the Strouhal number St and the corresponding a and using the 
interpolation factor. The result for use in equations (25), (26), and (30) is 


A(a) — ^min( a ) Aft(ao) [A max (a) ^min^)] 


(40) 


The function B in equation (27) and shown plotted in figure 78(b) is calculated in a manner similar to 
function A above. The two curves i? max and S m j n , through which B is obtained from interpolation, are 


B m m(b) = 


[ n/16.888 - 886.788 6 2 - 4.109 ( 6 < 0.13) 

-83.6076 + 8.138 (0.13 < 6 < 0.145) 

-817.8106 3 + 355.2106 2 - 135.0246 + 10.619 (0.145 < 6) 


(41) 


and 


where 


f \/l6.888 - 886.7886 2 - 4.109 ( 6 < 0.10) 

B m ax(6) = { -31.3306 + 1.854 (0.10 < 6 < 0.187) 

— 80.5416 3 + 44.1746 2 - 39.3816 + 2.344 (0.187 < 6) 

6 = | log(St s /St 2 )| 


(42) 


(43) 


The spectral shape B for intermediate values of R c have horizontal axis intercepts at —20 dB in figure 78(b) 
for values of 6 of 


M*c) = < 


0.30 (R c < 9.52 x 10 4 ) 

(-4.48 x 10 _13 )(fi c - 8.57 x 10 5 ) 2 + 0.56 (9.52 x 10 4 < R c < 8.57 x 10 5 ) 


(44) 


0.56 


(8.57 x 10 5 < Rc) 


The interpolation factor Bft(bo) is defined as 




20 -Bmin(^o) 


^max(6o) ^min(^o) 

and thus the result for use in equation (27) is 

B(b) — B m j n (b) + Bn(bo) [Bmax(6) — ^min(^)] 


(45) 


(46) 


61 



The amplitude function K\ in equations (25) and (26) is plotted in figure 77 and is given by 


( -4.31 log(flc) + 156.3 (Rc < 2.47 x 10 5 ) 


K i 


-9.0 log(ZZ c ) + 181.6 (2.47 x 10 5 < R c < 8.0 x 10 5 ) 

128.5 (8.0 x 10 5 < R c ) 


( 47 ) 


The level adjustment previously mentioned for the pressure-side contribution for nonzero angles of attack 
appears as AK\ in equation (25). This is given by 


AK[ = 


a* 1.43 log (/?,,*) — 5-29j 


0 


(Rs; < 5000) 
(5000 <R S .) 


(48) 


where R$* is the Reynolds number based on pressure-side displacement thickness. 

The amplitude function K 2 of equations (27) and (30) is plotted for some values of M in figure 82 and is 
given as 

f -1000 (a* < 70 “ 7) 


K 2 = K 1 + < 


v //?2-( / 3/ 7 )2 (a ,- 7o )2 +A) 


(70 - 7 < a* < 70 + 7) 


(49) 


-12 


(70 + 7 < a*) 


where 


7 = 27.094M + 3.31 
/ 3 = 72.65 M + 10.74 


70 = 23.43M + 4.651 
00 = -34.19 M - 13.82 


(50) 


The angle definitions above are in units of degrees and are taken as positive in sign. The K 2 definition 
above is valid for all values of a*, even when the calculation of the total noise in equation (24) switches from 
the use of equations (25), (26), and (27) for assumed attached TBL flow to equations (28), (29), and (30) for a 
supposedly stalled flow condition. The angle where the switch occurs, specified previously as (a*)o, is taken to 
be equal to the peak of the K 2 function defined by 70 in equation (50) or whenever a* exceeds 12.5°, whichever 
is first. 


5.1.3, Comparison With Data 

The scaling predictions of TBL-TE and separa- 
tion noise are compared with the noise data in fig- 
ures 11 to 43 for the tripped BL airfoils. The cal- 
culations used the appropriate values of 6 * from sec- 
tion 3 and the directivity functions from appendix B 
(where r e = 1.22 m, © e = 90°, and $ e = 90°). The 
total self-noise is given as well as the individual noise 
components of TBL TE noise from the suction and 
pressure sides and separation noise. The predictions 
follow the shapes and levels of the data, especially 
for the larger airfoils and the lower angles of attack 
where the scaling accuracy was most emphasized. 
Predictions of TBL TE and separation noise are also 
shown for the untripped BL airfoils in figures 44 


to 74. For the many untripped cases where these 
sources are predicted to be dominant, the agreement 
is generally good. Even where the LBL VS noise 
dominates, the TBL-TE and separation contribu- 
tions help with the overall spectral agreement. 

5.2, Laminar-Boundary-Layer-Vortex- 

Shedding Noise 

As previously described in section 1, laminar- 
boundary-layer instabilities couple with acoustic 
feedback to produce quasi-tonal noise. In contrast to 
TBL-TE noise, there are no LBL VS noise scaling 
methods established in the literature because of the 
erratic behavior of the multiple tones in the narrow- 
band spectra and the general complexity of the 
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mechanism. Two key results from the literature 
which provide initial scaling guidance are (1) that 
the gross trend of the frequency dependence was 
found to scale on a Strouhal basis, with the relevant 
length scale being the laminar-boundary-layer thick- 
ness at the airfoil trailing edge (ref. 16), and (2) that 
on the basis of the limited data from the data base of 
the present paper as reported in reference (6), overall 
levels tended to coalesce to a unique function of R c 
when normalized in the fashion of TBL-TE noise. 

The scaling approach taken herein is similar to 
that taken for TBL TE noise in the last section in 
that a universal spectral shape and Strouhal depen- 
dency is modeled in terms of boundary- layer param- 
eters, Mach number, angle of attack, and Reynolds 
number. The use of 1/3-octave spectra, rather than 
narrow band, permits such an approach because the 
broad spectral bands overlap the tonal frequency 
spacing to give smoother and generally single-peaked 
spectra. 

5.2.1. Scaled Data 

Scaled 1 /3-octave sound pressure level spectra for 
four airfoil sizes, each at four tunnel speeds, are 
presented in figure 83 from figures 44, 53, 59, and 65. 
The angle of attack for all is zero and the boundary 
layers are untripped. The normalization employs 

Scaled SPL 1/3 = SPL 1/3 - 10 log (m 5 ^) (51) 

for level scaling and 

St' = Ip- (52) 

for Strouhal frequency scaling. For the symmetric 
airfoils at zero angle of attack, 6 p = 6 S = <5o- 
The scaling approach differs from the TBL TE noise 
scaling because of the use of 6 p , the boundary- layer 
thickness on the pressure side of the airfoil, rather 
than <5*, the boundary- layer displacement thickness 
on the suction side. The use of 6 p as the pertinent 
length scale follows from reference 16 and was found 
to give seemingly better results in initial scaling of 
the present data base than 6* and by far better than 
c, <5 5 , or 6* for angles of attack other than zero. 

In figure 83(a) for the large 30.48-cm-chord air- 
foil, the spectra appear to be of smooth broad hump 
shapes. There is no apparent contribution to the 
spectra from LBL VS noise which is peaked in char- 
acter. The boundary layers are fully turbulent in 


the vicinity of the trailing edge at all four tunnel 
speeds (ref. 21), so no laminar vortex shedding is 
established. The noise produced is assumed to be 
all TBL-TE noise. In figure 83(b) for the 15.24-cm- 
chord airfoil, the broad spectral shapes are changed 
by the addition of a peak when the flow velocity is 
diminished. The peak levels increase with decreasing 
velocity. Although the boundary layer is turbulent at 
the trailing edge at all velocities shown, laminar flow 
exists over larger portions of the airfoil at the lower 
velocities. As mentioned for the LBL-VS noise mech- 
anism, any spectral peaks containing a number of 
tonal contributions should scale with Strouhal num- 
bers based on boundary- layer thickness. This is the 
case in figure 83(b) with St' » 0.27. For the shorter 
10.16-cm-chord airfoil, in figure 83(c), the LBL VS 
noise peaks become even more dominant for decreas- 
ing velocity. Note also the changing Strouhal depen- 
dence, not noted in previous studies. The shorter 
5.08-cm-chord airfoil, in figure 83(d), has even more 
pronounced level and Strouhal dependence with ve- 
locity variations. 


Whereas figure 83 shows the dependence of LBL 
VS noise on velocity for the various airfoil sizes at 
zero angle of attack, figure 84 shows the effect of 
angle of attack a* of the airfoils at a velocity of 
71.3 m/s. The spectra for the 30.48-cm-chord airfoil, 
shown in figure 84(a), change from being dominated 
by TBL-TE noise, for a* =0°, to being dominated 
by LBL-VS noise, for a* = 4.0°. So even with a 
large Reynolds number ( R c = 1.52 x 10 6 ), LBL-VS 
noise occurs. With increasing a*, the boundary layer 
on the pressure side becomes more laminar over a 
sufficiently large portion of the chord to result in in- 
creased shedding and corresponding noise. For the 
15.24-cm-chord airfoil (Rc = 7.58 x 10 5 ), shown in 
figure 84(b), the LBL-VS noise increases with a * un- 
til a certain value is reached where it diminishes. At 
a* = 7.2°, no apparent shedding noise is shown. At 
a* = 9.9°, the noise changes appreciably to that for 
stalled flow as discussed in the last section. The use 
of 6 p as the characteristic length scale apparently re- 
sults in a proper Strouhal scaling for the shedding 
noise peaks; but, as expected, the spectra for a* = 
0°, 7.2°, and 9.9°, which are dominated by TBL- 
TE and separated flow noise, diverge in this nor- 
malized format. A similar angle-dependent behavior 
where spectra do not coalesce is seen for the 
10.16-cm-chord airfoil, in figure 84(c), where LBL 
VS noise is apparent at a* = 0° and 3.3° but not 
at the higher angles. For the 5.08-cm-chord model, 
figure 84(d) shows large-amplitude LBL VS noise at 
a* = 0° and 4.2°. 
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The LBL-VS noise portions of the spectra 
(figs. 83 and 84) are rather invariant with respect 
to spectral shape. Based on this observation, a func- 
tion Gi (shown in fig. 85) was chosen as a shape to 
represent the LBL-VS contribution to the self-noise 
1/3-octave spectra for all cases. The level of G\ at 
St' = Stp eak has a value of —3.5 dB. The reference 
level of 0 dB is the integrated total of G\. To permit 
an orderly study of the Reynolds number and angle 
dependences of the spectral data, the shape G\ was 
matched to the individual spectra to obtain reference 
overall peak levels and Strouhal numbers. Emphasis 
was placed on matching the global spectral shape of 
G\ to the data rather than matching 1 /3-octave band 
peak or overall levels. Reference peak locations are 
indicated by the symbols in figures 83 and 84. 

In figure 86, the chosen values of Stp eak are plot- 
ted versus the Reynolds number R c for the 42 cases 
where LBL-VS noise is prominent. The values are 
also given in table 2 (at the back of this report) 
along with the effective angles of attack a* corre- 
sponding to at- For a* = 0, Stp eak is approximately 
constant at low R c and increases with R c in the mid- 
range of Rc shown. The values of Stp eak are lower 

for nonzero angles of attack. A function Stj is drawn 
to approximate the data of zero angle of attack. A 
constant value for Stj is chosen for high R r . where 
no zero-angle-of-attack data are present, because the 
value permits a simple angle dependence definition 
for Stp eak . In figure 87, Stp eak is normalized by St j 
and plotted versus a*. For each of the six airfoils, 
the line described by io _0 °4a* approximates the an- 
gle dependence. 

5.2.2. Calculation Procedures 


The reference peak scaled levels which correspond 
t° Stp eak in figure 86 are plotted versus R c in fig- 
ure 88. To show general trends more clearly, the 
symbols are replaced by the value of a*, rounded off 
to the nearest whole degree (see table 2 for more ex- 
act values). In this format it is seen that for each 
a* the scaled levels tend to increase, peak, and de- 
crease as Rc increases. For the larger angles of at- 
tack, the peak levels are lower and the corresponding 
values of R c are larger. Superimposed on the data 
are curves of identical shape, called here “level shape 
curves,” which are positioned in a monotonically de- 
creasing fashion to approximately correspond to the 
data trends with angle variation. The angles indi- 
cated for each curve position should not necessarily 
match the angle values listed for the data because 
the data values are rounded off in the figure, as men- 
tioned. The intent is to use the curves, with their 
functional relationship to a* and R c shown in fig- 
ure 88, to represent the amplitude definition of LBL 
VS noise. In the following calculation procedures sec- 
tion, a function G 2 specifies the curve shape, G3 is 
the angle dependence for the level of the G 2 curve, 
and a reference ( Ac )o value is defined as a function 
of angle to specify the Reynolds number dependence. 
The success of the functions in normalizing the data 
is shown in figure 89 where peak scaled 1/3-octave 
level minus G3 is compared with the function G 2 . 
In this format the individual angle numbers should 
ideally match the G 2 curve. Although the agreement 
shown is certainly not complete, it is regarded here as 
acceptable. Note that much better curve fits to the 
data would be possible if a requirement for mono- 
tonic functional behavior had not been imposed on 
G 3 and (R c )o- 


The LBL VS noise spectrum in a 1/3-octave presentation is predicted by 


SPLlbl-VS = 10 log 



_Rc_ 

(Rc) 0 . 


+ G' 3 (<**) 


The Strouhal definitions are (see figs. 86 and 87) 



0.18 ( Rc < 1.3 x 10 5 ) 

St'j = | 0.001756RO- 3931 (1.3 x 10 5 < Rc < 4.0 x 10 5 ) 

l 0.28 (4.0 x 10 5 < Rc) 


(53) 


(54) 


(55) 
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and 


St 


/ 

peak 


= St', X io-°- 04a ‘ 


(56) 


The directivity function D f, is given by equation (Bl) in appendix B. The function Gi defines the spectral 
shape, in terms of the ratio of Strouhal number to its peak, as (see fig. 85) 


r 39.8 log(e) - 11.12 (e < 0.5974) 

98.409 log(e) + 2.0 (0.5974 < e < 0.8545) 


G i(e) = l -5.076 


^2.484 - 506.25[log(e)] 2 


(0.8545 < e < 1.17) 


-98.409 log(e) + 2.0 (1.17 < e < 1.674) 

„ -39.8 log(e) - 11.12 (1.674 < e) 


(57) 


where e = St'/Stp eak . The peak scaled level shape curve G 2 depends on Reynolds number and angle and is 
(see figs. 88 and 89) 


G 2 (d) = { 


77.852 log(d) + 15.328 
65.188 log(d) + 9.125 
— 1 14.052[log(rf)] 2 
— 65. 188 log(d) + 9.125 
-77.852 log(d) + 15.328 


(d < 0.3237) 
(0.3237 < d < 0.5689) 
(0.5689 < d < 1.7579) 
(1.7579 < d < 3.0889) 
(3.0889 < d) 


where d = R c /(R c )q and the reference Reynolds number is 


(Rc) o = 


|Q0.215a*+4.978 

^QO.^Oa.+s^ea 


(a* < 3.0) 
(3.0 < a*) 


The angle-dependent level for the shape curve is 


G 3 (a*) = 171.04 - 3.03a* 


(58) 


(59) 


(60) 
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5.2.3. Comparison With Data 

The spectral predictions from the above equations 
are compared with the untripped BL airfoil noise 
data in figures 44 to 74. The great sensitivity of 
this mechanism to angle and velocity change can 
be clearly seen. In many respects the prediction 
agreement in shape, level, and actual occurrence of 
LBL-VS noise is good. Also as indicated in the 
last section, the combined contributions of LBL VS, 
TBL-TE, and separation noise are important to the 
total predictions for this untripped BL airfoil data. 

5.3. Tip Vortex Formation Noise 

The prediction method proposed in this section 
for tip vortex formation noise is that developed by 
Brooks and Marcolini (ref. 18). The study isolated 
this high-frequency broadband self-noise by compar- 
ing aerodynamic and acoustic test results of both 
two-dimensional (2D) and three-dimensional (3D) 
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airfoil models shown in figures 2 and 3, respectively. 
The premise of the tip noise determination method 
was that 3D models produce both tip noise and TBL- 
TE noise, while the 2D models produce only the lat- 
ter. The study produced a prediction method in gen- 
eral agreement with the physical model of the mech- 
anism first proposed by George, Najjar, and Kim 
(ref. 17). The noise is associated with the turbu- 
lence in the locally separated flow region at the tip 
of a lifting blade, where the tip vortex is formed. The 
flow field is illustrated in figure 90 for an airfoil blade 
tip at an angle of attack ajIP to the flow of velocity 
U . The flow over the blade tip consists of a vortex 
of strength T with a thick viscous core whose span- 
wise extent at the TE is l. The recirculating flow 
within the core is highly turbulent . The mechanism 
of noise production is taken to be TE noise due to 
the passage of this turbulence over the edge and into 
the wake. 


The tip vortex formation noise spectrum in a 1 /3-octave presentation is predicted by 


SPL tip = 10 log 


{Ml M^D j^ 


- 30.5(logSt" + 0.3) 2 + 126 


(61) 


The Strouhal number is 

St" = JL- (62) 

L'max 

The directivity function D h is given by equation (Bl) in appendix B. The second term on the right side of 
equation (61), which gives the frequency dependence, is a parabolic fit about a peak Strouhal number of 0.5. 
The spanwise extent at the TE of the separation due to the tip vortex is, for the tested rounded tip, 


^/c ~ 0.008axip 


(63) 


where c is the chord length and a X ip ( see discussion below) is the angle of attack of the tip region to the 
oncoming flow. The maximum Mach number M ma x of the flow within or about the separated flow region at 
the trailing edge is 

Mmax/M » (1 + 0.036a XI p) (64) 

where M is the Mach number of the oncoming flow to the airfoil tip region. The velocity corresponding to 
A/max is 

Ui nax = .Mmax (65) 


Note that in the use of equations (63) and (64) to determine l and M max ,a:TIP is correctly regarded as 
the actual angle of attack of the tip to the oncoming flow when the blade under consideration has a large 
aspect ratio (large span), is untwisted, and encounters uniform flow over its span. This is the reference case 
in reference 18. When the tip loading characteristics differ from those for the reference case, such as for some 
rotor and propeller blades, a T jp must be redefined according to computed sectional loading. The redefined 


QfXIP IS 


q TIP - 


( djjjdy \ 
\{dL'/ dy) re f ) y-.TIP 


(66) 
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where c*xip is the geometric angle and Li is the lift per unit span at the spanwise position y. The sectional 
lift slope dL f /dy near the tip is taken to be proportional to the tip vortex strength F (of fig. 90). When 
tip loading is found to be high, the predicted tip noise levels increase. The use of o^ip rather than axip 
in equations (63) and (64) generalizes the solution for arbitrary aspect ratios, blade twist, and spanwise flow 
variations. Reference 18 contains examples which provide guidance in the evaluation of equation (66) for aspect 
ratio, as well as tunnel testing effects. 

The tip noise prediction equations are based on data from airfoils with rounded blade tips. Of interest is a 
flat (square-off or cut-off) tip geometry which reference 8 considered, along with rounded tips, in calculations 
employing limited tip flow measurements reported in the literature. The different tip geometries required a 
different definition of the separated flow region size t. In applying the prediction equations of the present paper 
for flat tips, it does not appear appropriate to use the definition of reference 8. The constants in equation (61) 
reflect the definition of £ in equation (63). The measurements did not confirm the definition of £ for rounded 
tips proposed by reference 8. For consistency, the following definition for £ is proposed for flat tips for the 
present prediction equations: 


f 0.0230 + 0.0169qj IP (0° < a f Tlp < 2°) 

e/c = { (67) 

{ 0.0378 + 0.0095a' rip (2° < a' Tlp ) 


This definition of l approximately accounts for differences between the definition of reference 8 and that of 
equation (63) for rounded tips. There is at present no experimental confirmation of equation (67). 


5.3.2. Comparison With Data 

Noise data from reference 18 (fig. 7) are presented 
in figure 91 along with predictions of tip noise and 
the combined contributions of TBL TE and separa- 
tion noise. The rounded tip 3D model has a chord of 
15.24 cm and a span of 30.5 cm. The corresponding 
2D model has a span of 45.7 cm so its noise spec- 
trum levels in the figure were adjusted downward by 
1.8 dB (based on a 101og(L) dependency) to obtain 
that expected for a 30.5-cm span. The difference be- 
tween the 2D and 3D spectra should be that due 
to tip noise. The predictions in figure 91 for TBL- 
TE and separation noise, which employed the an- 
gle a* = 0.5(10.8°) to account for the wind tunnel 
correction, should ideally match the 2D model spec- 
trum. The tip noise prediction adds to the prediction 
to obtain a total which should match the 3D model 
spectrum. The tip noise prediction involved the use 
of equation (66) because of the finite extent of the 
span as well as open wind tunnel influences. Based 
on the lift distributions presented in reference 18, 
the tip angle becomes a^j P = 0.71(10.8°). While a 
slight overprediction at higher frequencies is seen in 
figure 91 for this particular example, the differences 
between levels with and without tip noise are the 
same for both data and prediction. The comparison 
shows consistency and compatibility not only with 
the data but also between the self-noise prediction 
methods. 


5.4. Trailing-Edge-Bluntness-Vortex- 
Shedding Noise 

In this section, the experiment of reference 2 is 
briefly described, published and previously unpub- 
lished TE bluntness noise data from the study are 
presented, and a prediction method is developed. 

5. 1^.1. Experiment 

The Brooks-Hodgson experiment (ref. 2) em- 
ployed an experimental arrangement similar to that 
reported in section 2 of the present paper with re- 
spect to hardware and acoustic measurement. How- 
ever, in reference 2, the model airfoil tested was large 
with a 60.96-cm chord length. When BL tripping was 
used, 2.0-cm-wide strips of No. 40 grit were applied 
at 15 percent of the chord. Rather than the TE be- 
ing sharp, the model TE thickness, or bluntness, was 
h = 2.5 mm. Figure 92 shows the TE region of the 
airfoil. The TE geometry was rounded at the two 
edges and flat between the rounded edge portions, 
which each comprised about one-third of the 2.5-mm 
thickness. The thickness h was varied, with edges 
of similar geometry, by alternately attaching exten- 
sions on the edge, as illustrated in figure 92(a). Also 
tested were sharp-edge ( h = 0) plate extensions 15.24 
and 30.48 cm long, as shown in figure 92(b). Another 
sharp-edge extension (not shown) was a 2.54-cm-long 
“flap” extension placed at 17.5° off the chord mean 
axis at the trailing edge. In addition, blunt plate ex- 
tensions were tested which were 15.24 cm long with 
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h = 2.5 and 4.8 mm and 30.48 cm long with h = 
4.8 mm. These extensions with rounded TE corners 
are shown in figure 92(c). Tape, 0.08 mm thick, was 
used to provide a smooth surface transition from the 
airfoil to the extensions. 

Presented in figure 93, from reference 2, are power 
spectral noise data of the airfoil at four flow veloc- 
ities. The airfoil is at zero angle of attack and the 
boundary layers are tripped. The microphone ob- 
server position is r e = 1.22 m and 0 e = 90° with re- 
spect to the model trailing edge. For two speeds, the 
spectra are given for the four TE thicknesses of fig- 
ure 92(a). The spectral results for the sharp, h = 0, 
TE cases should be all due to TBL-TE noise. The 
bluntness contributes additively at high frequencies 
to the spectrum levels. The values given for h/6* 
in figure 93 differ slightly from those specified in 
reference 2 because 6* here is calculated from the 
BL thickness scaling equations of the present paper. 
Data are presented in reference 2 for the sharp ge- 
ometries of figure 92(b), as well as the mentioned 
17.5° sharp flap extension. These geometries give es- 
sentially the same spectra as the sharp extension of 
figure 92(a). This demonstrates that TBL-TE noise 
is rather invariant with regard to geometry changes 
in the edge region, as long as the TE is sharp and 
the boundary layers are substantially the same. 

Trailing-edge bluntness noise spectra in a 
smoothed 1 /3-octave format are presented in fig- 
ure 94 for the edge geometries of figures 92(a) and 
92(c). These spectra are the result of a spectral 
subtraction process between the total spectra and 
the corresponding sharp TE spectra and should thus 
represent the bluntness contribution only. With the 
exception of the eight spectra also represented in 
figure 93, the data have not been previously pub- 
lished. The indicated values of h/6* for the ex- 
tensions are based on calculations of 6* for the TE 
of the airfoil without the extensions. This is justi- 
fied by indications that the boundary layers did not 
substantially change over the zero pressure gradient 
extension plates due to the influence of the upstream 
adverse pressure gradient (ref. 2). The spectrum for 
the airfoil with h = 2.5 mm and h/6 * = 1.15 in fig- 
ure 94 is for naturally transitional boundary layers; 
all others are for tripped boundary layers. 

5.4-2. Scaled Data 

The spectra of figure 94, as well as limited fre- 
quency data of Blake (ref. 19), form the foundation 
of the scaling approach. As with the scaling ap- 
proach for TBL-TE and LBL-VS noise, the level, 
frequency, and spectral shape are modeled as func- 
tions of flow and geometric parameters. For the level 
and frequency definition, we chose the peak of the 


spectral humps as the reference. The peak value of 
Strouhal number, defined as 

cjn _ /peak^ /* 0 \ 

^ l peak - JJ — 

is plotted versus the thickness ratio h/6* in fig- 
ure 95. The Strouhal numbers increase with increases 
in thickness ratio. The Strouhal numbers for the 
plate extensions of figure 92(c) are uniformly higher, 
for the same thickness ratios, than for the edge ex- 
tensions of figure 92(a). Also shown are two results 
obtained from Blake (ref. 19). Blake presents sur- 
face pressure data for a large array of plate edge 
geometries — all for very large values of h/6* (with 
the exception of the ref. 2 data reported and the one 
case shown in fig. 95 at h/6 * = 5.19). Blake, for most 
data, employed Strouhal relationships which depend 
on special wake stream thicknesses, and convection 
velocities not available without measurements. From 
Blake, however, it is obvious that different TE ge- 
ometries have different frequency dependences, con- 
sistent with the result of figure 95 that Strouhal num- 
bers for the flat plate extension and the airfoil TE 
geometries differ. The primary difference between 
the geometries is that the NACA 0012 airfoil has a 
beveled or sloping surface upstream of the trailing 
edge with a solid angle $ of 14° and the flat plate 
has F = 0°. The result shown from Blake in fig- 
ure 95 at h/6* = 5.19 is for a plate with F = 12.5° 
and nonrounded TE corners. In figure 95, parallel 
curves are fitted to the data. The curves, designated 
with values of F, are defined on the basis of a match 
point at h/6* = 20 for F = 0°. From Blake’s scaling 
for a thick flat plate (h/6* large) with nonrounded 
TE corners, one can determine that fh/U = 0.21 
at h/6* — 20. The curve for F = 14° intercepts 
Blake’s F = 12.5° result, but this is deemed an ac- 
ceptable deviation from the curve fit. For scaling 
purposes, values of Stp" ak for F values other than 0° 
and 14° could be determined by linear interpolation 
as described in the calculation procedure section to 
follow. 

For amplitude scaling, the peak values of the 1/3- 
octave spectra of figure 94 were normalized as 

M bb hL\ 

r e ) 

(69) 

The 5.5 power for Mach number dependence was 
determined to give better overall scaling success than 
either a 5 or 6 power. Figure 96 shows the scaled 
levels plotted versus the thickness ratio h/6*. As 
in figure 95 for the Strouhal dependency, the scaled 
levels are uniformly higher for the plates than for the 


Scaled peak SPL^ = Peak SPLj^ - 10 log ^ 
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h = 2.5 mm 



1.9 mm 

l h = 1 .1 mm h = 0 (sharp trailing 
" ~~ x ' edge) 



*- T railing-edge extensions 


(a) TE Extensions, Which are Alternately Attached 
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Figure 92. Illustration of trailing-edge extensions and plates. Smooth surface transition is provided for all geometries. 
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Blunt TE (h = 2.5 mm) 

Blunt TE (h = 1 .9 mm) 

Blunt TE (h = 1.1 mm) 



Figure 93. Spectral density for TE noise for 60.96-cm-chord airfoil with various degrees of TE bluntness. Tripped BL; at = 0° 
0 e = 90°. Level referenced to 1-Hz bandwidth. Data from reference 2. 
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Plate with h = 4.76 mm 
Length = 30.48 cm, M = .206 
h/5* = 1.16 


Plate with h = 4.76 mm 
Length = 15.24 cm, M = .206 
h/6* = 1.16 


Airfoil with h = 2.54 mm: 


Plate with 
h = 2.54 mm 
Length = 15.24 cm 
M = .206 
h/6* = .62 


Frequency, kHz 


\ 


I Airfoil: 

\ — h = 1 .9 mm 
| M = .203 
. h/6* = .46 

' h = 1 .1 mm 

\ M = .203 

' h/ 5* = .27 


h = 1 .9 mm, 

M = .113, h/6* = .40 


— h = 1.1 mm, 

M = .113, h/6* = .23 


Figure 94. TE-bluntness-vortex-shedding noise, extracted from data of figure 93, data for untripped BL, and data with plate 
extensions (fig. 92(c)) attached. 













edge extensions for the same thickness ratios. The 
levels increase with increasing thickness ratios. The 
edge extension data for the two smaller thicknesses 
of h = 1.1 and 1.9 mm at M = 0.113 deviate most 
from a straight line trend. Because of signal-to- 
noise concerns in the specification of these points, 
these data have the least confidence in the figure 
and are thus ignored in the specification of a curve 
fit. However, the accuracy of the resultant scaling 
equations in predicting these data is subsequently 
examined. The curve fits, designated as G±(h/8* , #), 
shown for the data are straight lines which are chosen 
to level off at h/6* = 5. The curve fit behavior at 
high h/6* is admittedly rather arbitrary, but there 
are no noise data available for guidance, unlike in the 
above Strouhal scaling where some frequency data 
from Blake are used. Fortunately, in practice, the 
likely values of h/6* to be found for rotor blades and 
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wings should be in the range where data are present 
and scaling confidence is greatest. 

Given the specification of the functions St^ ak and 
G 4, a definition of the spectral shape completes the 
scaling. Spectral curve fits for the data of figure 94 
are shown for the airfoil TE extensions, 4/ = 14°, and 
for the plate extensions, V = 0°, in figures 97(a) and 
97(b), respectively. The shapes reflect the observa- 
tions that the spectra are sharper for the plates for 
the same h/6*, and the spectra widen in the lower 
frequencies for decreased h/6* values for both the 
plates and the edge extensions. The spectral curve fit 
is specified as the function G*,(h/6*, $) whose peak 
level is 0 dB and whose shape is defined in terms of 
St'7St"' eak . The specification of G5 for in-between 
values of ^ would be an interpolation between the 
limiting cases shown in figures 97(a) and 97(b). 


The TE bluntness noise spectrum in a 1 /3-octave presentation is predicted by 


SPLblunt = 10 log 



(70) 


The directivity function D h is given by equation (Bl) in appendix B. The Strouhal definitions are (see fig. 95) 


/// _ f_h 

' “ U 


(71) 


and 


c*'" - 

peak 


( 0.212 - 0.0045fr 

I 1 + 0.235 (/i/6| vg ) _1 - 0.0132 (ft/^lvg) 2 


( 0.2 < fc/£ vg ) 


(72) 


l 0.1(h/<5* vg ) + 0.095 - 0.00243'!' 


(Wvg < 0.2) 


The V^avg term is the ratio of TE thickness 
thickness £* vg , where 


(degree of bluntness) h to the average boundary- layer displacement 


6 


* 

avg “ 


6* P + 6*s 
2 


(73) 


The angle 'I' is the solid angle, in degrees, between the sloping surfaces upstream of the trailing edge. For an 
edge on a flat plate * = 0°, whereas * = 14° for an NACA 0012 airfoil. The determination for this parameter 
for other TE geometries is discussed in section 6 and appendix C. 

The peak level of the spectrum is determined from the function G 4 (see fig. 96) where 


G 4 (Wv g,«) 


' 17.5 log (h/<5* vg ) + 157.5 - 1.114* (h/6* vg < 5) 

, 169.7-1.114* (5<V«; vg ) 


(74) 


The shape of the spectrum is defined by the function G 5 (see figs. 97(a) and 97(b)) where the calculation 
procedure involves an interpolation between the spectra for * = 0° and 14° as follows: 


G 5 



St' 1 


St"' 

3 t peak ; 


— (<^5)^=0° + 0.0714* [(G5 )^, =14 o — (G 5)\f/ = o°] 
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where 


{G$)y =U o = \ 


f TUT] + k 

2 . 50 ^ 7 ^ 2 - 2.5 


(»7 < m) 

(m<v< o) 


v/1.5625 - 1194.997/2 - 1.25 (0 < r/ < 0.03616) 

t -155.5437? + 4.375 (0.03616 < r?) 


v = iog(st"7st"; ak ) 


m = < 


0.1221 

(h/6* vs < 0.25) 

— 0.2175(/i/<5* vg ) + 0.1755 

(0.25 < h/6* vg < 0.62) 

-0.0308(h/<5* vg ) + 0.0596 

(0.62 < h/S* vg <1.15) 

0.0242 

(1.15 < h/S* vg ) 

0 

(h/S* vg < 0.02) 

68.724 (h/S* vg ) - 1.35 

(0.02 < h/S * vg < 0.5) 

308.475(/i/^av g ) - 121.23 

(0.5 < h/6* avg < 0.62) 

224. 811(/i/<5* vg ) -69.35 

(0.62 < h/6* vg < 1.15) 

1583.28(/i/<5* vg ) - 1631.59 

(1.15 < h/6£ vg < 1-2) 

268.344 

(1.2 < h/6t vg ) 


and 


VQ = 


m 2 /x 4 


6.25 + m 2 /z 2 


fc = 2.5 \ / 1 


(?) 


mr/o 


The spectrum (G'5)ar=o° is obtained by computing equations (76) through (81), as one would for (Gs)^ : 
but replacing (h/6* vg ) by (h/t 5 avg )' where 


h 


J avg 


6.724 


- 4.019 


av g , 


h 

,^av g> 


+ 1.107 


(76) 


(77) 


(78) 


(79) 


(80) 


(81) 
= 14 °) 

(82) 
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5.4-4- Comparison With Data 

Noise spectra for the airfoil with different TE 
thicknesses (geometry of fig. 92(a)) are presented for 
the flow Mach numbers of M = 0.21 and 0.12 in 
figures 98 and 99. respectively. The data were ob- 
tained by digitizing the spectra of figure 93 and con- 
verting these to 1/3-octave levels. The prediction 
curves shown are those of TBL TE and bluntness 
noise sources. For the sharp TE of figures 98(a) 


and 99(a), there is no bluntness contribution. Over- 
prediction is seen for the TBL TE noise at the low- 
est frequencies and some underprediction is appar- 
ent in the higher frequencies for the highest flow 
speed. For the nonzero TE thicknesses the blunt- 
ness noise contributes to the total spectra at high 
frequencies and renders good comparisons with the 
data. Good agreement is found even for the afore- 
mentioned smaller thickness cases at low Mach num- 
ber (figs. 99(c) and 99(d)). 
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Figure 98. One- third- octave presentation of spectra of figure 93 at U = 69.5 m/s with predictions for various degrees of bluntness. 
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Figure 99. One-third-octave presentation of spectra of figure 93 at U = 38.6 m/s with predictions for various degrees of bluntness. 
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6* Comparison of Predictions With 
Published Results 

The scaling law predictions are compared in this 
section with data from self-noise studies of airfoil sec- 
tions performed at the United Technologies Research 
Center (UTRC). 

6.1. Study of Schlinker and Amiet 

Schlinker and Amiet (ref. 3) conducted tests in 
the UTRC Acoustic Research Tunnel to study TBL- 
TE noise from a cambered helicopter blade section. 
The cross section of the 40.6-cm-chord and 53.3-cm- 
span model is shown in figure 100. As in the present 
NASA Langley studies, the model was mounted on 
sidewalls and spanned the width of the open tun- 
nel jet, so that the flow across the model was two- 
dimensional. The nozzle providing the flow had a 
rectangular exit of dimensions of 29 cm x 53.3 cm. 
To isolate the TBL-TE noise from facility back- 
ground noise, a directional microphone system was 
used. The experimental configuration, illustrating 
the shear layer refraction effect on the TE noise re- 
ceived by the directional microphone, is shown in fig- 
ure 101. The Mach numbers tested ranged from 0.1 
to 0.5 and the tunnel angle of attack at varied from 
—0.4° (zero lift for this cambered airfoil) to 12°. 

6.1.1. Boundary ‘Layer Definition 

Because only TBL TE noise measurements were 
desired, the boundary layers were tripped by apply- 
ing thin serrated aluminum tape at the blade loca- 
tions indicated in figure 100. The tape thickness 
was on the order of the BL displacement thickness 
at the points of application, providing minimum sur- 
face protrusion to avoid unnaturally large TBL thick- 
nesses downstream. This “light” trip is in contrast 
to the present study where the trips were “heavy” for 
reasons discussed. 

Hot-wire measurements were made in the 
boundary- layer/ near- wake region at the TE of the 
model. In figure 102, measured BL thicknesses are 
plotted versus Mach number for various tunnel angles 
of attack a*. These data are from figure 17 of refer- 
ence. 3. At zero lift, at = -0.4°, in figure 102(a), the 
BL thicknesses 6o on the pressure and suction sides 
are approximately the same. This should be expected 
since they developed under approximately the same 
adverse pressure gradient. Included in figure 102(a) 
are corresponding values of BL displacement thick- 
nesses, which were calculated by the present authors 
from velocity profiles presented in reference 3 (6* was 
not a quantity of interest in ref. 3). In figures 102(b) 
and 102(c), 6/c values are shown for at = 7.6° and 
12°, respectively. Comparing figures 102(a), 102(b), 


and 102(c), one can see that as angle of attack in- 
creases, increases and 6 p decreases. 

These measurements are compared with the 
thickness scaling equations of the present paper. 
First, equations (2) and (3) are used to calculate the 
BL thickness ratio 6q/c and displacement thickness 
ratio <5q/c. To make the calculations agree with the 
data of figure 102(a), all calculated values of 6q/c 
and £q jc were multiplied by a factor 0.6. This factor 
is taken to be the adjustment in equations (2) and 
(3) needed to make them appropriate for the “light” 
trip of reference 3. Next, the corrected angles of at- 
tack are determined by (1) adding 0.4° to at so that 
the tunnel angle is referenced to the zero-lift case 
and (2) using equation (1), with c = 40.6 cm and 
H = 79 cm, to obtain a* = 0°, 3.9°, and 6.1° for 
at — -0.4°, 7.6°, and 12°, respectively. These val- 
ues of a* are now used in equations (8) and (11) to 
obtain 6 p /6q and 6 s /6q, respectively. The resultant 
values of fj 9 /c and 6 p /c are compared with the data 
in figures 102(b) and 102(c). 

6.1.2. Trailing- Edge Noise Measurements and 

Predictions 

Trailing-edge noise spectra in a 1/3-octave pre- 
sentation are given in figure 103 for the airfoil at 
at = —0.4° with Mach number ranging from M = 
0.1 to 0.5. The data were obtained by the direc- 
tional microphone system at differing orientations to 
the airfoil. Shear layer corrections and directional 
microphone gain adjustments were made so that the 
data shown represent the noise radiated from a unit 
length of L = 0.3048 m of the TE span, at an ob- 
server distance of r e = 3 m, and an observer angle 
0 e which is specified in the figure. Figures 104 and 
105 contain spectra for the airfoil at at = 7.6° and 
12°, respectively. 

The TBL-TE and separation noise spectra were 
predicted using the calculation procedures of the 
present paper. The values of a*, <5*, and 6 p used 
were calculated as described in the previous section. 
Because of the BL trips and the 2D flow, no LBL-VS 
or tip noise calculations were made. In performing 
the calculations for TE blunt ness noise, one has to 
assign values of the TE thickness h and the TE flow 
angle parameter 'F. The thickness was indicated 
in reference 3 to b e h = 0.38 mm but the shape 
of this small TE region was not given. A value of 
$ = 17° has been used in the prediction because it 
gives reasonable prediction-data comparisons. 

In figures 103 to 105, the predictions are com- 
pared with the measurements. As in the presentation 
of figures 11 to 74, the individual noise contributions 
are shown, along with the total summed spectra. 
The prediction-data comparisons are good, especially 
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Suction surface 
trip location 



Figure 100. Cross section of Sikorsky rotor blade (ref. 3). Span is 53.3 cm and chord length is 40.6 cm. 



microphone 

reflector 

Figure 101. UTRC experimental configuration of reference 3, showing the effect of tunnel flow and shear layer refraction on the 
directional microphone alignment. 
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Frequency, kHz 
(c) U = 102.2 m/s 



Frequency, kHz 
(b) U = 68.1 m/s 



Frequency, kHz 
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Frequency, kHz 
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Figure 104. Noise spectra for Sikorsky airfoil at a* = 7.6° (a* = 3.9°) from reference 3 compared with prediction of present 
paper. 
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(a) U = 68.1 m/s (b) U = 102.2 m/s 


Figure 105. Noise spectra for Sikorsky airfoil at a* = 12° (a* = 6.1°) from reference 3 compared with prediction of present paper. 
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considering that the predictions are based empirically 
on a different airfoil section and that the noise mea- 
surement methods were quite different. There does 
appear to be a mild overprediction of the TBL TE 
noise, although not consistently so. The extent of 
agreement in the spectra where the TE bluntness 
noise contributes is substantially due to the afore- 
mentioned choice of = 17° (the previously used 
'F = 14° would result in a contribution about 3 dB 
higher than that shown). 

6,2. Study of Schlinker 

The tests of Schlinker (ref. 26) were similar in 
design to that of reference 3, whose measurement 
configuration is shown in figure 101. The 2D airfoil 
model, however, was an NACA 0012 section (as in the 
present study) with a chord length of c = 22.9 cm. 
Again, the aim of the tests was to measure TBL TE 
and not LBL-VS noise. However, no BL trip was 
used at zero angle of attack because no LBL VS noise 
was identified (except for the lowest speed tested and 
those data were not presented). At at = 6°, the LBL- 
VS noise was pervasive so a trip was placed on the 
pressure side at 30 percent of the chord to eliminate 
the LBL-VS noise. 

The TE noise spectra at various tunnel velocities 
are shown in figures 106 and 107 for the airfoil at 
at = 0° and 6°, respectively. The data were pro- 
cessed so that the levels shown are for the full airfoil 
span of L = 53.3 cm and an omnidirectional observer 
positioned at r e = 2.81 m and © e = 90°. For this air- 
foil, the corrected angles of attack, using equation (1) 
with c = 22.9 cm and H = 79 cm, are a* = 0° and 
3.9° for at = 0° and 6°, respectively. The predic- 
tions shown in figure 106 for zero angle of attack are 
for TBL-TE, LBL-VS, and TE bluntness noise. The 
values of and 6q used in the predictions were ob- 
tained from equations (5) and (6), for an untripped 
BL airfoil. The predictions shown in figure 107 for 
at = 6° are for only TBL-TE, separation, and TE 
bluntness noise, since the LBL-VS noise was elim- 
inated by the pressure side tripping. The required 
values of <5* were calculated from equation (14), for 
an untripped BL. However, the values of 6* were de- 
termined from equations (3) and (9), for a tripped 
BL and then multiplying the result by 0.6 (to re- 
flect the “light” trip condition as discussed for the 
Schlinker and Amiet study). For the calculations for 
TE bluntness noise, there was no guidance from the 
paper for the specification of h and 'F. A reasonable 
TE thickness of h = 0.63 mm was assumed and the 
TE flow angle parameter was set at 4* = 23°, because 
it gave good agreement for the high frequencies in fig- 
ures 106 and 107. The overall agreement between the 
total predictions and the data appears good. 
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6.3. Study of Fink, Schlinker, and Amiet 

Fink, Schlinker, and Amiet (ref. 27) conducted 
tests in the UTRC tunnel to study LBL-VS noise 
from three airfoil geometries. The untripped BL air- 
foil models had an NACA 0012 planform and their 
geometries are shown in figure 108. The first had 
a constant-chord length of 11.4 cm across the span 
while the other two were spanwise tapered, having 
linearly varying chord lengths along the span. Of 
the tapered airfoils, the first had a taper ratio of 2 
to 1 with chord length varying from 15.2 cm down 
to 7.6 cm. The other airfoil had a taper ratio of 4 to 
1 with chord length varying from 18.3 cm to 4.6 cm. 
The span was L = 79 cm for all cases. Because the 
levels of the LBL-VS noise were sufficiently intense 
compared with the tunnel background noise, a direc- 
tional microphone system was not used to measure 
the noise. Instead, far-field spectra were obtained 
with individual microphones placed on an arc of 
2.25-m radius about the midspan of the models. The 
noise data from reference 27 which are presented 
in the present report are all from a microphone for 
which 0 e « 90°. 

Reference 27 presented most noise data in narrow- 
band form at various bandwidths to allow exami- 
nation of the tonal character of the LBL-VS noise. 
To compare these data with the predictions of the 
present paper, the narrow-band data were digitized 
and converted to 1/3-octave presentations. As a 
check on this procedure, as well as a check on the con- 
sistency of the data presented in reference 27, overall 
sound pressure levels (OASPL) were computed from 
the digitized data and compared with overall levels 
reported from direct measurement. The values gen- 
erally agreed to within 1.0 dB. 

For the constant-chord airfoil at at = 4°, 1/3- 
octave spectra are shown in figure 109 for various 
tunnel velocities between U = 37 m/s and 116 m/s. 
The number of spectral bands, as well as the fre- 
quency range, presented for the spectra varies for the 
different speeds. This variation is due to the differ- 
ent narrow-band analysis ranges used in reference 27, 
as all available data were used to generate the 1/3- 
octave band spectra. For U — 37 m/s, figure 109(a), 
the spectrum is flat at the lower frequencies but is 
peaked between 1 and 3 kHz. From the narrow-band 
presentation of reference 27 (fig. 22), one finds that 
the flat portion is dominated by broadband noise, 
which is characteristic of tunnel background contam- 
ination. It is noted again that these spectra are single 
microphone results from which the background noise 
has not been subtracted. The spectral peak region 
is due to the presence of three quasi-tones, repre- 
senting the LBL-VS noise portion. At U = 52, 64, 
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Figure 108. Airfoil models of reference 27. All dimensions are in centimeters. 


and 79 m/s, in figures 109(b), 109(c), and 109(d), 
the spectra are very peaked because of the dominat- 
ing contributions from large numbers (10 to 15) of 
LBL-VS quasi-tones. At U = 98 and 116 m/s, in 
figures 109(e) and 109(f), the spectra are less peaked 
because of a somewhat decreased number of quasi- 
tones which become submerged within broadband 
background noise (which itself increases with speed). 

The strong velocity dependence of the noise is 
seen clearly in figure 110 (from fig. 25 of ref. 27) 
where the OASPL is plotted as a function of veloc- 
ity. The overall levels were directly measured, for the 
noise between 200 Hz and 20 kHz, rather than deter- 
mined by integrating measured spectra. The levels 
rise and then stabilize with increases in velocity. The 
resumed increase in levels at the highest speeds (ap- 
proximately 100 m/s) is where the background noise 
appears to become dominant. 

Compared with the data in figures 109 and 110 
are noise predictions of LBL-VS, TBL-TE, and sep- 
aration noise. No consideration was given to blunt- 
ness noise because of the lack of information about 
the TE geometry as well as the fact that LBL VS 
noise dominates the predictions where comparative 
data are available. For the BL thickness determina- 


tions, the equations of section 3 for untripped bound- 
ary layers were used. The corrected angles of attack 
were calculated from equation (1), with c = 11.4 cm 
and H = 53 cm, which rendered a* = 0° and 1.9° 
for at = 0° and 4°, respectively. These were em- 
ployed with the prediction equations for an observer 
at r e — 2.25 m, 0 e = 90°, and 4> e = 90°. The 
predictions in figure 109 give good comparisons, ex- 
cept that the peak frequencies are lower than pre- 
dicted. The previously described background noise 
contributions explain the differences for the lowest 
and highest speeds. For the predictions of OASPL in 
figure 110, the spectra for LBL-VS, TBL-TE, and 
separation noise were summed. Predictions are pre- 
sented for not only at = 4° but also at = 0°, 2°, and 
6°. This is done to show the great sensitivity of the 
predictions to airfoil angle of attack. It is seen that 
the data would most agree with predictions for about 
a t ~ 5° rather than at = 4°. This could be inter- 
preted to mean that the agreement is on the order of 
possible experimental bias error in angle definition. 

The tapered-chord airfoils were used in refer- 
ence 27 to provide a continuous variation in expected 
vortex tone frequency to compare with an analogous 
rotating constant-chord blade. The tone variation 
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Figure 110. Overall sound pressure level versus velocity for constant-chord airfoil from reference 27 compared with predictions 
of present paper. 




was found not to be continuous; however the ta- 
pered models did produce spectra containing a large 
number of peaks spread over a somewhat wider fre- 
quency range than those for the constant-chord air- 
foil at about the same velocities. In figure 111, 1 73- 
octave spectra are shown for the 2-to-l taper air- 
foil at a t = 4° for tunnel velocities between U = 27 
and 107 m/s. The data are similar to those for the 
constant-chord model, except that the peaks are gen- 
erally less well defined. In figure 112, corresponding 
OASPL variations with tunnel velocity are shown for 
a* = 4°. Also in this figure, OASPL is shown for a 
range of velocities where at = 0°. 

The predictions shown in figures 111 and 112 were 
obtained by dividing the models into 10 segments of 
constant chord (where actual chord length for each 
segment varied according to the blade taper), then 
making predictions for each segment, and summing 
on a pressure-squared basis the contributions of each. 
Angle-of-attack corrections for each segment were 
made by calculating the correction based on the mean 
chord (11.4 cm) across the span. This correction 
was then applied to the angle of attack for each of 
the blade segments. The corrected angles, therefore, 


were the same as for the constant- chord model, that 
is, a* = 0° and 1.9° for at — 0° and 4°, respectively. 
The comparisons between predictions and data for 
the 2-to-l taper airfoil appear about as good as those 
for constant-chord comparisons. It appears that the 
predictions for OASPL at at = 4° (fig. 112) would 
best agree if a t « 3.5° had been used rather than 4°. 
This again indicates that agreement is on the order 
of possible experimental angle definition error. The 
OASPL comparisons for zero angle of attack show 
the predicted trends to be quite good but the levels 
to be overpredicted by 5 to 7 dB. 

In figures 113 and 114 are the data and prediction 
comparisons for the 4-to-l taper model at at = 0°. 
The predictions are not as good as for the constant- 
chord and the less tapered model, although the data 
still fall within a predictive range of at = 2° to 3°. 
One should bear in mind that the flow behavior in 
the vicinity of the tapered models would be expected 
to deviate from the idealized 2D behavior assumed 
to be occurring over the small spanwise segments 
employed for the predictions. This makes it difficult 
to assess the meaning of the comparison deviations 
for the tapered models. 
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Figure 112. Overall sound pressure level versus velocity for 2-to-l tapered-chord airfoil from reference 27 compared with 
predictions of present paper. 
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Figure 114. Overall sound pressure level versus velocity for 4-to-l tapered-chord airfoil from reference 27 compared with 
predictions of present paper. 
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7. Conclusions 

This paper documents the development of an 
overall prediction method for airfoil self-noise. The 
approach is semiempirical and is based on previous 
theoretical studies and data from a series of aero- 
dynamic and acoustic tests of isolated airfoil sections. 
The acoustic data processing employed a correlation 
editing procedure to obtain self-noise spectra uncon- 
taminated by extraneous noise. Five self-noise mech- 
anisms, each related to specific boundary- layer phe- 
nomena, are identified in the data and modeled. For 
each mechanism, the data are first normalized by fun- 
damental techniques using scaled aerodynamic pa- 
rameters. The spectral shape, level, and frequencies 
are then examined and modeled for dependences on 
parameters such as Reynolds number, Mach number, 
and geometry. 

The modeling accuracy of the resulting self-noise 
prediction methods is established by comparing pre- 
dictions with the complete data base. The methods 
are shown to have general applicability by compar- 
ing predictions with airfoil self-noise data reported in 
the literature from three studies. A successful appli- 
cation of the methods is reported for a large-scale- 
model helicopter rotor broadband noise test. 

Conclusions can be drawn regarding the spe- 
cific self-noise mechanisms. For the turbulent- 
boundary-layer-trailing-edge noise and separation 
noise sources, an accurate and generally applica- 
ble predictive capability is demonstrated, especially 
for the important conditions of high Reynolds num- 
ber and low to moderate angle of attack. The 
mechanism which can dominate the spectra for low 
Reynolds number, laminar-boundary- layer- vortex- 
shedding noise, is also demonstrated to have good 
predictive capability. For this quasi-tonal noise 


mechanism, there are some issues, not fully ad- 
dressed herein, about how to apply the formulations 
in the most appropriate way to different airfoil ge- 
ometries. The tip vortex formation noise source ap- 
pears to be well predicted, although its relative lack 
of importance compared with the other self-noise 
sources prevents a full assessment of accuracy. The 
trailing-edge-bluntness-vortex-shedding noise source 
is shown to be very important and predictable by the 
method developed. For this source, there is an asso- 
ciated “flow angle” parameter which is found to be 
constant for any given trailing-edge geometry, but 
is difficult to determine a priori. However, for ap- 
plication of the bluntness noise prediction method, 
reasonable estimates for this parameter can be made 
based on the examples in this report. 

The unique prediction capability presented should 
prove useful for the determination of broadband noise 
for helicopter rotors, wind turbines, airframe noise, 
and other cases where airfoil shapes encounter low- 
to moderate-speed flow. For modern propeller de- 
signs, the present equations should be applied with 
some caution because the high-speed, high-loading, 
and skewed-flow conditions existing about propeller 
blades do not match the low- to moderate-speed and 
generally 2D flow conditions of the present data base. 
The computer codes given herein can be readily in- 
corporated into existing or future noise prediction 
codes. The documentation provided in this report 
should provide the means to evaluate where and how 
any needed future refinements can be made in the 
prediction codes for particular applications. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
April 19, 1989 
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Appendix A 

Data Processing and Spectral 
Determination 

In section 4, the special processing approach used 
to determine the self-noise spectra for the 2D airfoil 
models was summarized. Details are given here. 

A.l. Data Acquisition and Initial Processing 

Signals from the microphones shown in figure 4 
were recorded during the test on a 14-channel FM 
analog tape recorder, operated to provide a flat fre- 
quency response up to 40 kHz. Individual amplifiers 
were used to optimize signal-to-noise ratio for each 
microphone channel, and pure-tone and white-noise 
insertions were used to calibrate amplitude and phase 
response, respectively. These calibrations and signal- 
conditioning techniques were the same as in refer- 
ence 2, where additional details are given. The data 
were reduced from tape on a spectrum analyzer in- 
terfaced with a minicomputer. Pairs of microphones 
were used to obtain 1024-point cross-correlations at 
an analysis range of ±4.167 milliseconds. 

A. 2. Correlation Editing 

The correlation records are modified to eliminate 
contributions from extraneous noise sources prior to 
taking the Fourier transforms to obtain the spectra. 
The first step is to remove, to the extent possible, 
the noise from the test hardware by subtracting the 
correlation #45(7-) without the airfoil in place (the 
background noise) from i?45(r) with the airfoil in 
place. (See fig. 9.) The resulting record should 
then be comprised of correlation peaks from the 
desired TE noise, LE noise, and other extraneous 
noise related to interaction between the model and 
test rig not accounted for in the subtraction. The 
TE and LE noise peaks in the cross-correlation are 
assumed to represent the autocorrelation of the TE 
and LE noise, respectively. 

To eliminate the LE contribution, the correlation 
record on the right-hand-side (RHS) of t xe is dis- 
carded and replaced by the mirror image of the left- 
hand-side (LHS). However, for this folding process, 
it was found that it is important to preserve the ba- 
sic shape of the TE peak to more accurately repre- 
sent the spectra at higher frequencies. Because this 
is a digital correlation, made up of discrete points 
which are At apart, it is likely that the true TE 
noise peak falls somewhere between two discrete val- 
ues of r. Folding about a discrete point instead of the 
actual effective peak center would introduce error by 
distorting the peak shape. In figure Al, the discrete 
points of the TE correlation peak are illustrated to 


show how the folding was accomplished. The dis- 
crete center is at t X e, whereas the effective center is 
to the left. The correlation values at t xe + At and 
r TE + 2 At must not be changed to avoid modifying 
the shape near the very peak. The correlation value 
at tje + 2 At is projected to the LHS to intercept a 
line connecting t xe - 3 At and t xe - 2 At. This de- 
fines the constants a and b which are shown. These 
constants then are used to interpolate between the 
points on the LHS to determine values at the points 
on the RHS, that is 

R (tte + A At) = — — R (tje — N Ar) 

At 

+ ^i?(r TE - (A^ + l)Ar) (Al) 

for N > 2. The entire LHS of the correlation is 
folded about the effective peak center using this 
interpolation scheme. 


Effective peak 
center 



Figure Al. Sample correlation peak. 

A.3. Separation of TE and LE Peaks 

As indicated in section 4, for some of the correla- 
tion data for the three smaller chord lengths, the LE 
and TE peaks are so close that the LE contribution 
overlaps and distorts the TE peak shape. For many 
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such cases a procedure was found to successfully re- 
move the distortion prior to implementing the TE 
peak folding process described above. This procedure 
is explained by way of example for the 5.08-cm-chord 
airfoil shown in the bottom trace of figure A2. The 
predicted locations of the TE and LE noise peaks in 
the correlations are indicated and agree well with the 
peaks in the actual trace. Note the proximity of the 
two peaks. 

The procedure to separate these two peaks in- 
volves combining the original #45 (r) at the bottom 
of figure A2 with time-shifted versions of itself, so 
that the peaks are separated by larger time delays. 
The procedure depends on the implied symmetry of 
the LE and TE peaks, inherent in the assumption 
that they represent the autocorrelations for the LE 
and TE noise, respectively. The first step is to invert 
i?45 (r) in sign, and reverse it in time, by “flipping” 
the correlation about tle- The result of combining 
these two curves is seen in the second trace from the 
bottom of figure A2, denoted i? 45 (r). The two peaks 
seen here are the original TE noise peak, and an in- 
verted TE noise peak at 2 tl£ — rj £. There is some 
increase in level and some distortion in the correla- 
tion record away from the peaks, as should be ex- 
pected. The LE noise peak has been removed, but 
the inverted TE noise peak still affects the original 
peak at tje* To remove the inverted peak, the initial 
R^{t) must be shifted by 2 (tee - 7te) and summed 
with the previous result. This produces the third 
curve from the bottom in figure A2, denoted /2 45 (t). 
The TE noise peak has remained at tte? while the 
LE noise peak is now at 3 tle — 2tte- The peaks are 
now separated in time so that details of each peak 
can be seen. Note that as the peaks no longer affect 
one another’s shape, their basic symmetry is evident. 
This helps to validate the initial assumption that the 
peaks represent the autocorrelations of TE and LE 
noise. If the peaks must be further separated, this 
procedure can be successively repeated, with the re- 
sults of the next two iterations seen in the top two 
traces of figure A2, i? 45 (r) and i? 45 (r). It should 
be noted that only the inner portion of the corre- 
lation is shown (the correlation was performed for 
±4.067 ms). Because of the data record manipula- 
tions, much of the outer portions of the correlations 
did not overlap and were thus zeroed out. 

A.4. Determination of Spectra 

Once the correlation records, or their modified 
forms after the separation processing, are folded 
about the effective peak center, the resulting TE 
noise correlations are transformed into spectra of 
the noise. Because the correlation record lengths 
are reduced by varying amounts (typically 20 per- 


cent) because of the editing described above, the use 
of fast Fourier transform techniques is not conve- 
nient. Instead, regular Fourier transform techniques 
are used in an approach based on chapter 9 of ref- 
erence 28. In summary, a data window is applied 
to the correlation (eq. 9.116, ref. 28) and is used to 
provide the real and imaginary portions of the spec- 
trum (eqs. 9.167-9.168, ref. 28). The resulting cross- 
spectra (eqs. 9.172-9.174, ref. 28) are presented in 
terms of magnitude and phase. 

With the cross-spectra produced, amplitude cor- 
rections are applied to account for shear layer effects, 
using the technique of Amiet (ref. 22), as well as self- 
noise directivity effects, which are described in ap- 
pendix B. The spectrum for each microphone pair 
was corrected to an effective position of 90° with re- 
spect to the airfoil chord line. The combined effect of 
both of these corrections tended to be small, with the 
corrections for many test conditions being less than 
1 dB. Since cross-spectra were obtained, the correc- 
tions for each of the two microphones involved were 
averaged to correct the cross-spectral magnitude. 

The results obtained from this method are given 
in figure A3 for the example correlation records of 
figure A2. Figure A3(a) shows the cross-spectrum 
obtained from the correlation of the original R^ir) 
record, which is the bottom curve of figure A2, while 
figure A3(b) shows the cross-spectrum obtained after 
folding the /? 4 5(r) record about the TE noise peak. 
Note that the cross-spectral phase 0 is a partial indi- 
cator of how well the cross-spectrum represents the 
total TE self-noise. Ideally the phase should vary 
linearly with frequency, 4> = 360° / r TE- The breaks 
seen in this phase line and the corresponding spectral 
peculiarities indicate regions adversely affected by 
contamination, which was not removed by the back- 
ground subtraction and, in the case of figure A3(b), 
the folding process. The contamination from the LE 
is seen to primarily affect the cross-spectrum of fig- 
ure A3(a) below around 4 kHz. Folding the correla- 
tion removes most of this, leaving a dip in the spec- 
trum of figure A3(b) at about 1.5 kHz. Figure A3(c) 
shows the spectrum for the third curve from the bot- 
tom in figure A2, i? 45 (r), which is the modified corre- 
lation after two manipulations have separated the TE 
and LE noise peaks. The phase difficulty and spectral 
dip at about 1.5 kHz in figure A3(b) are eliminated 
in figure A3(c). Figure A3(d) shows the spectrum 
for the top curve of figure A2, i? 45 (r), which is for 
four manipulations. This spectrum is similar to that 
of figure A3(c) except for some apparent increase in 
contamination at the low- and high-frequency ends. 
For the airfoil presented here, a choice was made 
to use the spectrum of figure A3(c), based on two 
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Figure A3. Cross-spectra obtained by transforming correlations of figure A2. 








separation manipulations, to represent the self-noise. 
The lower and upper limits to which the spectrum is 
believed to be accurate are from about 0.8 to 13 kHz. 
For the other airfoils of this paper, similar evaluations 
of the limits were made and the spectra are cut off 
beyond these limits in their presentation, as indicated 
in section 4. 

To increase confidence, all the 2D airfoil spectra 
presented in figures 11 to 74 were found by averaging 
independently determined spectra from two micro- 
phone pairs. After the shear layer and directivity 
corrections were applied, the spectra from the two 
microphone pairs generally agreed to within 1 dB. 
In tables 1 and 2, the data processing and manip- 
ulations, and whether the correlations were folded 
or not prior to taking the spectra, are specified for 
each test case. It is seen that for the three larger 
airfoils, no correlation manipulations were needed to 
separate the LE and TE correlation peaks. For the 
three smaller airfoils, performing two separation ma- 


nipulations was advantageous for about a quarter of 
the cases. The table shows that a substantial num- 
ber of correlations were not folded. For airfoils at 
sufficiently high angles of attack, low frequencies can 
dominate the noise. This results in large correlation 
humps, rather than the relatively sharp peaks which 
are needed in the folding process. For these cases, the 
raw cross-correlations are transformed, with only the 
background subtraction being performed. Also the 
correlations were not folded in the presence of strong 
LBL-VS noise. This noise can dominate all other 
self-noise sources, as well as the LE noise contami- 
nation, negating the need for the correlation editing. 
This correlation editing would have proved difficult, 
in any case, since vortex shedding produces noise at 
small bands of frequencies, appearing as damped si- 
nusoidals in the correlation, which tended to mask 
other peaks. The effect of folding the correlation in 
such cases was examined, however, and found to have 
little effect on the spectra. 
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Appendix B 
Noise Directivity 

The purpose of this appendix is to define the 
directivity functions D ^ and D £, which are employed 
in the tunnel noise data processing and proposed 
for use in the prediction equations for the self-noise 
sources. 

B.l. Retarded Coordinates 

The retarded coordinate system is explained by 
first referring to figure B1 where the airfoil is at zero 
angle of attack to the tunnel flow. If the velocity 
were zero everywhere, sound from the model which 
reaches the microphones (M2 is shown) would follow 
the ray path defined by the measured distance r m 
and angle 0 m . But with the velocity in the free jet 
equal to U , the ray which reaches the microphones 
follows first the radiation angle 0 C until it encounters 
the shear layer where it is refracted. It emerges at 
angle 0^ with an amplitude change and travels to the 
microphone. The theoretical treatment employed in 
this study for the angle and amplitude corrections 
is that due to Amiet (refs. 22 and 23). A convenient 
reference for the corrected microphone measurements 
is a retarded coordinate system where the source 
and the observer are at corrected positions. The 
angle 0 e is referenced to a retarded source position 
and a corrected observer position where the distance 
between the positions is r e = r m . As defined, if 
there were no shear layer present with flow extending 
to infinity, the center of the wave front emitted 
from the source would be at the retarded source 
position when the wave front reaches the corrected 
observer position. The retarded coordinates are 
equivalent to the emission time coordinates employed 
in the literature, for example, see reference 29, for 
moving sources and stationary observers. Figure B2 
shows a source flyover geometry corresponding to 
the open jet wind tunnel geometry of figure Bl. 
Physical equivalence between the cases is attained by 
accounting for the Doppler-related frequency shifts 
due to the relative motion between the source and 
observer in one instance and no relative motion in 
the other. There are no Doppler-related amplitude 
corrections required between the flyover and wind 
tunnel cases as the effect of the flow on the source 
definition is already included in the wind tunnel 
environment. 

B.2. Directivity Functions 

In figure B3, a 3D retarded coordinate system is 
defined where the origin is located at the trailing edge 
of a thin flat plate, representing an airfoil. The flat 


plate is in rectilinear motion of velocity U in direction 
of the negative x e axis. The observer is stationary. 
Tr ailing-edge noise is produced when boundary- layer 
turbulence and its associated surface pressure pat- 
tern convect downstream (with respect to the plate) 
at a velocity U c (Mach number M c ) past the trail- 
ing edge. If the noise-producing turbulence eddies 
are sufficiently small and the convection velocities 
are sufficiently large to produce acoustic wavelengths 
much shorter than the chord length, the directivity 
can be shown to be (based on analysis of Schlinker 
and Amiet, ref. 3) 


D h (9 e ,*e) 


2sin 2 (0 e /2) sin 2 

(1 + A/cos© p )[l 4- (M — A/ r )cos© r ] 2 


(Bl) 


where the h subscript indicates the high-frequency 
(or large-chord) limit for D . The overbar on D ^ 
indicates that it is normalized by the TE noise ra- 
diated in the 0 e = 90° and = 90° direction, so 
£>^(90°, 90°) = 1. For the flyover plane ($*> = 90°), 
equation (Bl) is the same as equation (32a) of ref- 
erence 3. In reference 3, the equation was compared 
favorably with measured airfoil TE noise results, for 
limited M and 0 e ranges, as well as with previous 
theoretical results. The directivity expression used 
in reference 2 was found to give virtually identical 
results for low Mach numbers. 

Although developed for when the velocity U is 
parallel to the plate along the x e axis, equation (Bl) 
can be applied when the plate or airfoil is at an 
angle of attack a to the flow. In application (refer 
to fig. B3), one should define the angles with respect 
to a coordinate system that is fixed with respect to 
the airfoil with the x e axis fixed along the chord line, 
rather than one where the x e axis is fixed along the 
direction of motion. Note, however, that any analysis 
of Doppler frequency shifts (not treated in this paper) 
should reference angles with respect to the direction 
of motion. Applications of equation (Bl) at angles 
of attack should result in little additional error to 
that already built into the relation. Because it was 
derived with the plate assumed to be semi-infinite, 
Dfr becomes inaccurate at shallow upstream angles 
(0 e — ► 180°), when applied to finite airfoils even 
for high frequencies. As frequency is lowered, the 
wavelengths become larger with respect to the chord 
and the directivity becomes increasingly in error. 
However, D ^ should be of sufficient accuracy to 
define the directivity of all the self-noise sources 
discussed because of their high-frequency character. 
The one exception to this is the stalled-airfoil noise. 

When the angle of attack of the airfoil is increased 
sufficiently, the attached or mildly separated TBL 
flow on the suction side gives way to large-scale 
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Figure B2, Emission from a source moving with a constant velocity. 
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Figure B3. Flat plate in rectilinear motion. 


separation. The turbulence eddies are then compa- 
rable in size to the airfoil chord length and the eddie 
convection speeds are low. The directivity for this 
low-frequency noise is more properly defined as that 
of a translating dipole, which is 


Dt(Qe,*e) 


sin 2 © e sin 2 $ e 

(1 + M COS ©e ) 4 


(B2) 


where the l subscript indicates a low-frequency limit. 
The coordinate system and comments about an- 
gle definitions in equation (Bl) apply also in equa- 
tion (B2). Equation (B2) is employed for the di- 
rectivity in the expression for stalled flow noise 
(eq. (30)). 


For the noise data reduction in the present study, 
equation (Bl) was used in the determination of the 
self-noise spectral levels for the reference observer po- 
sition, at r e = 122 cm and 0 e = 90°. First, shear 
layer refraction corrections were calculated to deter- 
mine the spectral level adjustments, to add to mea- 
sured values, and a resultant source-observer location 
at r e and © e . This was done while keeping track of 
the actual physical coordinates of the trailing edge 
which varied with airfoil angle of attack. Finally, 
equation (Bl), with $ = 90° and an assumed con- 
vection Mach number of M c « 0.8 A/, was used to 
determine final level adjustments required to match 
results to the © e = 90° location. 
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Appendix C 

Application of Predictions to a Rotor 
Broadband Noise Test 

An acoustics test of a 40-percent scale model 
BO-105 helicopter main rotor was conducted in the 
German-Dutch Wind Tunnel (DNW). Figure Cl 
shows an overview of the test setup in the large open 
anechoic test section. The 4-meter-diameter rotor is 
shown positioned in the flow between the nozzle on 
the right and the flow collector on the left. A key aim 
of the test was to produce a large benchmark aero- 
acoustic data base to aid and verify rotor broadband 
noise prediction. In reference 30, the present authors 
compared data with predictions of rotor broadband 
self-noise for a number of rotor operating conditions. 
The predictions employed the self-noise prediction 
methods, which are documented in section 5 of the 
present paper, and the NASA ROTONET program 
(ref. 31) to define rotor performance and to sum con- 
tributions of noise from individual blade segments. 

In this appendix, the experiment is not reviewed 
in detail nor are data-prediction comparisons pre- 
sented, as reference 30 is complete in this regard. 
Rather, reference 30 is complemented by specifying 
how the self-noise prediction methods of the present 
paper were applied. Given below is a summary of 
the rotor prediction method, a definition of the rotor 
blade geometry and test modifications and a specifi- 
cation of input parameters for the individual source 
predictions. The degrees of success of data-prediction 
comparisons in reference 30 are discussed along with 
recommended refinements to the prediction methods. 

To produce a rotor prediction, the rotor geometry 
definition and flight conditions, specified as thrust, 
rotor angle, rotor speed, flight velocity, and trim con- 
dition, are provided as inputs to the ROTONET ro- 
tor performance module. The particular module used 
assumes a fully articulated rotor with rigid blades 
and a simple uniform inflow model. The module de- 
termines local blade segment velocities and angles of 
attack for a number of radial and azimuthal posi- 
tions. Ten radial segments were considered at 16 az- 
imuthal positions. The BL thicknesses and other pa- 
rameters needed are calculated. The noise due to 
each source is predicted for each blade segment, and 
the ROTONET noise radiation module is used to 
sum contributions from all blade segments to obtain, 
after accounting for Doppler shifts and the actual 
number of blades, the noise spectrum at the observer. 

As indicated in reference 30, the accuracy of 
predictions depends on a number of factors including 
the accuracy of the performance module used. One 
may question the quasi-steady assumptions used in 
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defining the local BL characteristics, which ignore 
unsteadiness and resultant hysteresis effects. Likely 
more important is how well the aeroacoustic scaling 
determined from low-speed data extends to higher 
speed. The Mach number at the tip of the blades is 
0.64 for rotor hover, whereas the 2D airfoil model 
tests were limited to Mach 0.21. Also there are 
questions on how to apply scaling obtained from 
symmetrical NACA 0012 sections with particular TE 
geometries to the cambered NACA 23012 rotor blade 
with different TE geometries. 

The model rotor is a 40-percent-scale, four- 
bladed, hingeless BO- 105 rotor, with a diameter of 
4.0 m and a chord of 0.121 m. A blade and its details 
are shown in figure C2. The blades have -8° linear 
twist and a 20-percent cutout from the hub center. 
The effects of several blade modifications were exam- 
ined, including (1) application of Carborundum grit 
from the blade leading edge to 20 percent chord to 
match the BL trip condition for the 2D blade sec- 
tions described in section 2 of the present paper, 
(2) taping of the TE with 0.064-mm-thick plastic 
tape to modify the “step tab 17 geometry, and (3) at- 
tachment of a rounded tip to each blade (the stan- 
dard blades have a squared-off tip). 

C.l. Boundary-Layer Definition 

With the local blade segment mean velocities and 
angles of attack determined by the rotor performance 
module, the equations of section 3 were directly ap- 
plied to determine the BL thicknesses required in the 
noise predictions. Most noise comparisons in refer- 
ence 30 are for the blades with untripped BL. For 
the tripped BL, the fact that the BL trip conditions 
for the rotor blades matched the 2D test models as- 
sured the appropriateness of using the equations for 
a heavy trip rather than modifying the equations as 
required for the UTRC comparisons reported in sec- 
tion 6. For all BL thickness calculations, the aero- 
dynamic angles of attack were used in the equations. 
The aerodynamic angle is referenced to the zero lift 
angle, which is —1.4° from the geometric angle for 
the NACA 23012 airfoil. 

C.2. TBL-TE and Separation Noise 

Prediction 

Given the definitions of segment chord length, 
span width, velocity, aerodynamic angle of attack, 
and BL thicknesses, the calculation of TBL-TE and 
separation noise is straightforward as specified in 
section 5. From the data-prediction comparisons of 
reference 30, it is concluded that the TBL-TE and 
separation noise calculations demonstrated a good 
predictive capability for these mechanisms. The 
rotor was tested from hover to moderately high flight 
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Figure Cl. Test setup in DNW for helicopter main rotor broadband noise study reported in reference 30. 
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speeds for various climb and descent rates at differ- 
ent thrust settings. Diagnostics included 1/2 rotor 
speed tests and the BL tripping tests. It is noted that 
the TBL-TE and separation noise predictions for a 
number of rotor conditions fell below contributions 
of LBL-VS, especially at the 1/2 rotor speed, and of 
TE bluntness noise. This represents a limitation of 
the comparisons which prevents sweeping statements 
regarding predictive accuracy of TBL-TE and sepa- 
ration noise sources. Still the agreements were quite 
good except when the rotor operated at full speed 
(tip speed of M = 0.64) and the boundary layers 
were tripped. Then the noise was underpredicted by 
about 6 dB. It is believed that for this high speed the 
heavy trip disturbed the flow substantially, made it 
dissimilar to the 2D model cases, where the speed 
was limited to M = 0.21, and perhaps changed the 
controlling noise mechanisms. Comparisons for the 
tripped BL rotor at 1 /2 speed and the untripped BL 
blades at full and 1/2 speed produced good results. 

C.3. LBL-VS Noise Prediction 

The comparisons for LBL-VS noise in refer- 
ence 30 showed, for a broad range of rotor condi- 
tions, very good predictions. As with the TBL-TE 
and separation noise predictions, the calculation of 
LBL-VS noise is straightforward given the specifica- 
tion of local flow conditions at the blade segments. 
A special note should be made for one key param- 
eter involved in the calculations. The angle of at- 
tack a* employed in the LBL-VS noise prediction 
(eqs. (53) to (60)) was the geometric angle rather 
than the aerodynamic angle for the NACA 23012 air- 
foil section. The BL thickness calculations, however, 
used the aerodynamic angle, as previously stated. 
The use of the geometric angle for the noise calcula- 
tion is justified by (1) the better rotor data-prediction 
comparisons found using the geometric rather than 
the aerodynamic angle and (2) the lack of guidance 
one has in applying the acoustic scaling laws which 
were based on symmetrical airfoil results, to airfoils 
that are cambered. Remember that the controlling 
mechanism of LBL-VS noise is the presence of aero- 
acoustic feedback loops between the trailing edge and 
an upstream location on the airfoil surface where 
laminar instabilities occur. This geometric connec- 
tion indicates that a purely aerodynamic angle def- 
inition for the LBL VS mechanism would not likely 
be correct. An alternate viewpoint of the angle defi- 
nition problem would be that the aerodynamic angle 
should be used, which would increase the noise pre- 
dicted over that measured, but that allowance should 
be made for the fact that the inflow to the rotor blade 
segments is not the assumed smooth quasi-steady 
flow. The presence of sufficiently unsteady flow con- 


ditions over portions of the rotor would prevent the 
establishment of the LBL-VS mechanism and related 
noise. Limiting LBL-VS noise production to some 
measure of inflow turbulence offers promise as a re- 
finement to the self-noise prediction method. 

C.4. Tip Vortex Formation Noise 

The tip noise predictions were made for both the 
rounded and the squared-off blade tips tested. The 
performance module was used to determine the local 
flow velocities and angles for the tips at different 
azimuth locations. The ajip used was the NACA 
23012 aerodynamic angle. Because the tip loading 
characteristics for the rotor blades differed from the 
reference case of the tip noise model, which was an 
untwisted large-aspect-ratio blade with uniform flow 
over the span, the sectional lift term of equation (66) 
was evaluated. The sectional lift slopes for the rotor 
blades were analyzed by employing a lifting-surface 
model adapted from reference 18. The velocity and 
angle of attack were linearly varied over the span near 
the tip of the lifting surface blade. It was found that 
the tip loading is increased over the reference case 
by a small amount. For equation (66), the redefined 
ajjp angle was then given by ttjjp = l.laxip- 

The predictions for tip noise in reference 30 were 
in all cases significantly below predictions for TBL 
TE noise. This makes it impossible to truly assess 
the accuracy of the tip noise modeling for the ro- 
tor. However, since the data comparisons with the 
total levels predicted were good for both low and nor- 
mal rotor speeds, the tip noise is apparently well pre- 
dicted. It is noted that a review of data for a number 
of rotor cases, not all given in reference 30, indicated 
no significant effect due to the blade tip modification. 
This is in line with prediction for this rotor. 

C.5. TE-Bluntness-Vortex-Shedding Noise 

Given the flow definition for the blade segments 
from the performance module, the bluntness predic- 
tions require the specification of thickness h and flow 
angle parameter ty. As with the UTRC test compar- 
isons of section 6, it is not clear how to apply scaling 
laws obtained from an airfoil with a particular TE 
geometry to a rotor blade with a different TE. For 
the step tab TE geometry, shown in figure C2, h was 
specified as the actual 0.9 mm and was taken as 
14°, which is actual solid angle of the surface at the 
TE of the NACA 23012 airfoil (same as the NACA 
0012 airfoil). However, because of the 0.5-mm step 
5 mm upstream of the TE, 0.5 mm was added to 
the calculated value of <5* vg to approximately account 
for the anticipated step-caused BL flow deficit. For 
the TE tape modification case, <5* vg was taken as 
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that calculated, because the step was removed, but 
h was increased by four tape thicknesses. Had the 
tape remained fully attached to the TE surface (see 
fig. C2) during the test, two thicknesses would have 
been added. The flow angle was taken as 18°. The 
choice of this specific number was rather arbitrary, 
but is in line with that used for the UTRC compar- 
isons (section 6) for rounded trailing edges. The tape 
rounded the TE bluntness which should reduce the 
persistence of and noise due to the separated flow in 
the near wake. The larger angle value (18° com- 
pared with 14°) results in less noise predicted. 

The comparisons of reference 30 obtained using 
the above “reasonable” choices for the TE parame- 
ters give good results for all 1/2 rotor speed cases. 


For the full rotor speed cases the levels were con- 
sistently overpredicted. This is believed to be due 
to a speed dependence for the blunt ness mechanism 
that could not have been anticipated from the low 
speed airfoil data, from which the scaling laws were 
developed. Subsequent analysis indicates that much 
better agreement with data could have been obtained 
if the bluntness noise contribution was eliminated for 
blade segments exceeding Mach numbers of 0.45 or 
0.5. This is in some conflict with comparisons in sec- 
tion 6 for the blade section noise of Schlinker and 
Amiet (ref. 3), which shows apparently strong blunt- 
ness noise at M = 0.43 and 0.5. However, based 
on the rotor results, an upper limit of 0.45 for the 
bluntness noise contribution is recommended as a re- 
finement to the prediction method. 
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Appendix D 
Prediction Code 

The airfoil self-noise prediction method is avail- 
able as a computer code written in standard 
FORTRAN 5 specifically for the Digital Equipment 
Corp. VAX-1 1/780 series machine running under the 
VMS operating system. To the extent possible, the 
code has been made machine independent. There is 
one input file to the code and one output file. In- 
put consists of user supplied NAMELIST parameters 
while output is a table of 1 /3-octave centered fre- 
quencies with corresponding sound pressure levels for 


each noise mechanism followed by their total. The 
user selects which of the mechanisms to calculate. 

The airfoil section for which a prediction is desired 
is assumed to be composed of a number of segments, 
each having its own chord, span, angle of attack, free- 
stream velocity, trailing-edge bluntness, and angle 
parameter, as well as observer directivity angles and 
distance. This permits a variety of configurations 
such as taper, twist, spanwise-varying free-stream 
velocity (for rotor blades), etc. The user may specify 
as many or few segments as desired depending on the 
complexity of the geometry. Characteristics for each 
segment are specified in the input file, which contains 
the FORTRAN variables given in table Dl. 


Table Dl. Segment Characteristics Specified in Input File 


FORTRAN name 

Symbol 

Description 

NSEG 


Number of segments 

C 

c 

Chord length, m 

l 

L 

Span, m 

R 

r e 

Observer distance, m 

THETA 


Observer angle from x-axis, deg 

PHI 


Observer angle from y-axis, deg 

ALPSTAR 

a* 

Aerodynamic angle of attack, deg 

ALPHTIP 

a TIP 

Tip flow angle, deg 

H 

h 

Trailing-edge blunt ness, m 

PSI 

* 

Trailing-edge angle, deg 

U 

U 

FYee-stream velocity, m/sec 

ITRIP 


0— Use untripped BL condition 

1 — Use tripped BL condition 

2 —Use lightly tripped BL condition 

IL AM 


1 — Compute LBL-VS noise 
0 — Do not compute LBL-VS noise 

ITURB 


1 — Compute turbulent TBL-TE noise 
0 — Do not compute TBL-TE noise 

IBLUNT 


1— Compute TE bluntness noise 
0 — Do not compute TE blunt ness noise 

ITIP 


1 — Compute tip noise 
0 — Do not compute tip noise 

IROUND 


.TRUE. -Use rounded tip in tip calculation 
.false. —Use square tip in tip calculation 

vise 

V 

Kinematic viscosity, m 2 /sec 

CO 

Co 

Speed of sound, m/sec 
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The prediction shown in figure 45(a) was obtained 
using the following input: 

$ INDATA 


NSEG 

= 

1, 

C 

= 

. 3048 , 

L 

= 

.4572, 

R 

= 

1.22, 

THETA 

= 

90. , 

PHI 

= 

90. , 

ALPSTAR 

= 

1.516, 

U 

= 

71.3, 

ITRIP 

= 

o, 

ILAM 

= 

1. 

ITURB 

= 

1, 


Send 

Note that all parameters need not be included in the 
input if their default values are desired (see program 
listing for default values). In this example, only the 
laminar and turbulent mechanisms are computed and 
the untripped boundary layer condition is used in 
both mechanisms. The airfoil consists of one segment 
of constant chord and the observer is 122 cm directly 
beneath the trailing edge at the midspan. The free- 
stream velocity has a constant value of 71.3 m/sec 
across the span. For this example, the output file is 
given in table D2. 

Similarly, the prediction shown in figure 91, was 
obtained using the following input: 

SINDATA 

NSEG = 10 . , 

C = 10*. 1524 

L = 10*. 0305 


R 

= 

10* 1.22, 

THETA 

= 

10* 90. , 

PHI 

— 

10* 90. , 

ALPSTAR 

= 

10* 5.4, 

ALPHTIP 

= 

7.7, 

U 

- 

10* 71.3, 

ITRIP 

= 

1, 

ILAM 

= 

0 , 

ITURB 

= 

1. 

ITIP 

— 

1. 

ROUND 

= 

.TRUE. , 


Send 


This is an example of a multisegmented case where 
each segment has the same geometry and inflow 
conditions. Turbulent-boundary-layer noise and tip 
noise are calculated where the tip is rounded and at 
an effective angle of attack of 7.7°. All segments 
are summed to yield a total prediction for each 
mechanism as shown in table D3. 

For the VAX-11/780 machine running under 
VMS, the following commands will compile, link, and 
execute the code (assumed to reside on 
PREDICT. FOR), read input from a file EXAMPLE. IN, 
and write results to a file EXAMPLE. OUT: 

$ FOR PREDICT 
$ LINK PREDICT 
$ ASSIGN EXAMPLE. IN F0R004 
$ ASSIGN EXAMPLE. OUT F0R005 
$ RUN PREDICT 

The details of execution for other machines or operat- 
ing systems may vary. A listing of the code follows. 
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Table D2. Output File From Prediction Code for Test Case of Figure 45(a) 


ONE-THIRD OCTAVE 
SOUND PRESSURE LEVELS 



PRESSURE 

SUCTION 

SEPARATION 





FREQUENCY { HZ ) 

SIDE TBL 

SIDE TBL 

SIDE TBL 

LAMINAR 

BLUNTNESS 

TIP 

TOTAL 

100.000 

20.654 

28 . 704 

-100.000 

-17.142 

0.000 

0.000 

29.336 

125.000 

24.461 

31.965 

-100 .000 

-13 . 285 

0 .000 

0.000 

32 . 676 

160 .000 

28 . 291 

35.244 

-75.254 

-9.018 

0 .000 

0 .000 

36 .042 

200.000 

31.437 

37.937 

-49.243 

-5.161 

0 .000 

0 .000 

38.815 

250.000 

34.309 

40 .400 

-27 . 506 

-1.304 

0 .000 

0 .000 

41.356 

315.000 

37.023 

42.736 

-9 .030 

2 .690 

0.000 

0.000 

43.768 

400 .000 

39 .577 

44.949 

6 .266 

6.820 

0 .000 

0 .000 

46 . 057 

500 .000 

41.761 

46.859 

17.532 

10.677 

0 .000 

0 .000 

48.034 

630 .000 

43.845 

48 . 706 

26.603 

14.671 

0 .000 

0.000 

49 . 954 

800.000 

45.839 

50.503 

33.718 

18 . 801 

0.000 

0.000 

51.849 

1000 .000 

47.581 

52 . 106 

38.756 

22.658 

0.000 

0 .000 

53.568 

1250 .000 

4 9 . 2 3 3 

53.664 

42.692 

26.515 

0.000 

0.000 

55.255 

1600 .000 

50.987 

55. 368 

46.294 

30.782 

0 .000 

0.000 

57.106 

2000 .000 

52.533 

56.907 

49.334 

37.725 

0 .000 

0.000 

58.817 

2500 .000 

54 .074 

57.750 

51.298 

47.262 

0.000 

0.000 

60.167 

3150.000 

55 . 570 

57.500 

50 . 766 

48.959 

0 .000 

0.000 

60.496 

4 000 . 000 

56.044 

56.082 

47.711 

4 1 .796 

0 .000 

0.000 

59.455 

5000 .000 

55.399 

54 . 541 

44 .617 

32.428 

0.000 

0.000 

58 .208 

6300.000 

53.840 

52 . 942 

40 .974 

28.433 

0.000 

0 .000 

56 .553 

8000.000 

52.190 

51 .253 

36.227 

24.304 

0 .000 

0 .000 

54.821 

10000.000 

50.638 

49.614 

30.419 

20.447 

0 .000 

0.000 

53.192 

12500 .000 

49 .044 

47.890 

22.834 

16 . 590 

0.000 

0 .000 

51 . 523 

16000.000 

47.202 

45.851 

11.842 

12.323 

0.000 

0 .000 

49.591 

20000.000 

45.436 

43.863 

-0.924 

8.466 

0 .000 

0.000 

47.731 

25000.000 

43 .549 

41.710 

-16.833 

4 .609 

0 .000 

0.000 

4 5 . 737 

31 500 .000 

41.440 

39.279 

-37.092 

0.614 

0.000 

0 .000 

43.503 

40000.000 

39.065 

36.522 

-62.593 

-3.515 

0.000 

0.000 

40.987 


Table D3. Output File From Prediction Code for Test Case of Figure 91 


ONE-THIRD OCTAVE 
SOUND PRESSURE LEVELS 


PRESSURE SUCTION SEPARATION 


FREQUENCY { HZ ) 

SIDE TBL 

SIDE TBL 

SIDE TBL 

LAMINAR 

BLUNTNESS 

TIP 

TOTAL 

100 .000 

19.913 

43.883 

-19.803 

0 .000 

0.000 

-34.005 

43.900 

125.000 

23.788 

46 .159 

-0 .396 

0 .000 

0.000 

-24.312 

46 .184 

160 .000 

27.673 

48.459 

16.851 

0 .000 

0 .000 

-14.255 

48.498 

200 . 000 

30.853 

50 .372 

29.124 

0.000 

0.000 

-5.769 

50.452 

250 .000 

33 .746 

52.155 

38.723 

0 .000 

0 .000 

2.145 

52.407 

315 .000 

36.470 

53.894 

46.334 

0 .000 

0 .000 

9 .738 

54.662 

400 .000 

3 9 . 0 2 4 

55.609 

52 .245 

0.000 

0.000 

16.940 

57.320 

500 .000 

41.202 

57.165 

56.460 

0.000 

0 .000 

23.074 

59.897 

630 .000 

4 3 .274 

58.766 

59.996 

0 . 000 

0 .000 

28.824 

62.489 

800 .000 

4 5 .252 

60.360 

6 3 . 297 

0 .000 

0 .000 

34.121 

65.130 

1000 .000 

46.980 

60.940 

65.719 

0.000 

0 .000 

38.475 

67.016 

1 250 .000 

4 8.6 2 0 

60 .473 

65.697 

0 . 000 

0 .000 

42.257 

66.917 

1 600 .000 

50.364 

58.874 

62 .909 

0 .000 

0 .000 

45.774 

64.582 

2000 .000 

51.911 

57 .328 

59.818 

0 .000 

0.000 

48.349 

62.363 

2500 .000 

53.456 

55.775 

56.383 

0 .000 

0.000 

50.351 

60.580 

3150.000 

54 .709 

54 .122 

51 .975 

0 .000 

0 .000 

51.821 

59.364 

4000 .000 

54 .799 

52.336 

45 . 974 

0 .000 

0.000 

52 .694 

58.443 

5000.000 

53 . 761 

50.565 

38.550 

0 .000 

0.000 

52.917 

57.439 

6 300 .000 

52.162 

48.597 

28 . 510 

0.000 

0.000 

52.544 

56.204 

8000.000 

50.507 

46.387 

15.081 

0.000 

0 . 000 

51 . 512 

54.736 

10000.000 

4 8 . 9 3 6 

44 , 132 

-0.755 

0 . 000 

0 .000 

49.955 

53.078 

12500 .000 

47.311 

41.665 

-20 . 241 

0.000 

0.000 

47.826 

51 . 110 

16 000 .000 

45.415 

38.655 

-46 .603 

0.000 

0 .000 

44.802 

48.594 

20000.000 

43 . 583 

35-650 

-75,275 

a .ooo 

0 .000 

41.466 

46 .075 

25000 . 000 

41.611 

32.347 

-90 .000 

0 .000 

0.000 

37 . 557 

4 3 .405 

31500.000 

39.390 

28 . 582 

-90 .000 

0.000 

0 .000 

32 .904 

40 . 555 

40000 . 000 

36.873 

24 .291 

-90 .000 

0 . 000 

0.000 

27.449 

37.552 
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0001 


PROGRAM PREDICT 



0002 






0003 






0004 

C 





0005 

c 


* * * * * 

VARIABLE DEFINITIONS ***** 


0006 

c 





0007 






0008 

c 

VARIABLE 

NAME 

DEFINITION 

UNITS 

0009 

c 





0010 






0011 

c 

ALPHTIP 


TIP ANGLE OF ATTACK 

DEGREES 

0012 

c 

ALPSTAR 


SEGMENT ANGLE OF ATTACK 

DEGREES 

0013 

c 

ALPRAT 


TIP LIFT CURVE SLOPE 

— 

0014 

c 

C 


SEGMENT CHORDLENGTH 

METERS 

0015 

c 

CO 


SPEED OF SOUND 

METERS/SEC 

0016 

c 

FRCEN 


1/3 OCTAVE CENTERED FREQUENCIES 

HERTZ 

0017 

c 

H 


SEGMENT TRAILING EDGE THICKNESS 

METERS 

0018 

c 

IBLUNT 


FLAG TO COMPUTE BLUNTNESS NOISE 

— 

0019 

c 

I LAM 


FLAG TO COMPUTE LBL NOISE 

— 

0020 

c 

ITIP 


FLAG TO COMPUTE TIP NOISE 

— 

0021 

c 

ITRIP 


FLAG TO TRIP BOUNDARY LAYER 

— 

0022 

c 

ITURB 


FLAG TO COMPUTE TBLTE NOISE 

— 

0023 

c 

L 


SEGMENT SPAN LENGTH 

METERS 

0024 

c 

MAXFREQ 


MAXIMUM NUMBER OF FREQUENCIES 

— 

0025 

c 

MAX S EG 


MAXIMUM NUMBER OF SEGMENTS 

— 

0026 

c 

NFREQ 


NUMBER OF 1/3 OCTAVE FREQUENCIES 

— 

0027 

c 

NSEG 


NUMBER OF SEGMENTS 

— 

0028 

c 

PI 


PRESSURE ASSOCIATED WITH 


0029 

c 



TBLTE PREDICTION 

NT/M2 

0030 

c 

P 2 


PRESSURE ASSOCIATED WITH 


0031 

c 



TBLTE PREDICTION 

NT/M2 

0032 

c 

P3 


PRESSURE ASSOCIATED WITH 


0033 

c 



TBLTE PREDICTION 

NT/M2 

0034 

c 

P 4 


PRESSURE ASSOCIATED WITH 


0035 

c 



TOTAL PREDICTION 

NT/M2 

0036 

c 

P 5 


PRESSURE ASSOCIATED WITH 


0037 

c 



LBLVS PREDICTION 

NT/M2 

0038 

c 

P 6 


PRESSURE ASSOCIATED WITH 


0039 

c 



BLUNTNESS PREDICTION 

NT/M2 

0040 

c 

P 7 


PRESSURE ASSOCIATED WITH 


0041 

c 



TIP NOISE PREDICTION 

NT/M2 

0042 

c 

PHI 


DIRECTIVITY ANGLE 

DEGREES 

0043 

c 

PSI 


BLUNTNESS ANGLE 

DEGREES 

0044 

c 

R 


SEGMENT TO OBSERVER DISTANCE 

METERS 

0045 

c 

ROUND 


LOGICAL INDICATING ROUNDED TIP 

— 

0046 

c 

SPL 


TOTAL SOUND PRESSURE LEVEL 

DB 

0047 

c 

SPLALPH 


SOUND PRESSURE LEVEL ASSOCIATED 


0048 

c 



WITH TBLTE PREDICTION 

DB 

0049 

c 

SPLBLNT 


SOUND PRESSURE LEVEL ASSOCIATED 


0050 

c 



WITH BLUNTNESS PREDICTION 

DB 

0051 

c 

SPLLBL 


SOUND PRESSURE LEVEL ASSOCIATED 


0052 

c 



WITH LBL PREDICTION 

DB 

0053 

c 

SPLP 


SOUND PRESSURE LEVEL ASSOCIATED 


0054 

c 



WITH TBLTE PREDICTION 

DB 

0055 

c 

SPL S 


SOUND PRESSURE LEVEL ASSOCIATED 


0056 

c 



WITH TBLTE PREDICTION 

DB 

0057 

c 

SPLTBL 


TOTAL PRESSURE LEVEL ASSOCIATED 


0058 

c 



WITH TBLTE PREDICTION 

DB 

0059 

c 

SPLTIP 


SOUND PRESSURE LEVEL ASSOCIATED 


0060 

c 



WITH TIP NOISE PREDICTION 

DB 

0061 

c 

ST 


STROUHAL NUMBER 

— 

0062 

c 

THETA 


DIRECTIVITY ANGLE 

DEGREES 

0063 

c 

U 


SEGMENT FREESTREAM VELOCITY 

METERS/SEC 

0064 

c 

vise 


KINEMATIC VISCOSITY 

M2/SEC 

0065 






0066 






0067 


PARAMETER 

(MAXSEG 

= 20, MAXFREQ = 27) 


0068 






0069 






0070 


DIMENSION 

FRCEN ( MAXFREQ ) ,C(MAXSEG) , L ( MAXSEG ) 

0071 


1 

ST ( MAXFREQ ) , S P LLB L ( MAX F REQ ) , SPLTBL ( MAXFREQ ) , 

0072 


2 

U ( MAXSEG ) , SPLP ( MAXFREQ ) , SPLS ( MAXFREQ ) , 

0073 


3 

SPLALPH (MAXFREQ ) , S P L ( 7 , MAX FREQ ) , R ( MAX S EG ) 

0074 


5 

SPLBLNT (MAXFREQ ) ,PHI(MAXSEG) , SPLTI P ( MAX FREQ ) , 

0075 


7 

THETA( MAXSEG ) , ALPSTAR ( MAXSEG ) ,PSI(MAXSEG) r 


0076 

8 

H (MAXSEG ) 

, PI (MAXFREQ ) 

,P2 (MAXFREQ) 

0077 

9 

P3 ( MAXFREQ ) 

, P4 (MAXFREQ ) 

, P5 (MAXFREQ ) 

0078 

0079 

0080 
0081 
0082 
0083 

1 

REAL 

LOGICAL 

P6 (MAXFREQ ) 
L 

ROUND 

, P7 (MAXFREQ ) 



0084 C DEFINE DEFAULT VALUES FOR NAMELIST DATA 

0085 C 

0086 

0087 

0088 DATA C / MAXSEG*1.0 / 

0089 DATA L / MAXSEG*. 10 / 

0090 DATA R / MAXSEG * 1. / 

0091 DATA THETA / MAXSEG * 90. / 

0092 DATA PHI / MAXSEG * 90. / 

0093 DATA ALPSTAR / MAXSEG * 0.0 / 

0094 DATA H / MAXSEG * .0005/ 

0095 DATA PSI / MAXSEG * 14.0 / 

0096 DATA U / MAXSEG * 100. / 


0097 


DATA 

ITRIP 

/ 

0 

/ 





0098 


DATA 

I LAM 

/ 

0 

/ 





0099 


DATA 

ITURB 

/ 

0 

/ 





0100 


DATA 

I BLUNT 

/ 

0 

/ 





0101 


DATA 

IT I P 

/ 

0 

/ 





0102 


DATA 

ALPHTIP / 

o 

o 

/ 





0103 


DATA 

NS EG 

/ 

10 

/ 





0104 


DATA 

vise 

/ 1 

. 4 5 2 9 E- 5 / 





0105 


DATA 

CO 

/ 

340 . 

46 / 





0106 


DATA 

ALPRAT 

/ 

1 . 0 

/ 





0107 


DATA 

ROUND 

/ 

.FALSE. / 





0108 











0109 


DATA 

NFREQ 

/ 

27 

/ 





0110 











0111 

C 

SET UP VALUES OF 

1/3 

OCTAVE 

CENTERED 

FREQUENCIES 



0112 

C 

— 



— 






0113 











0114 


DATA 

FRCEN 

/ 100 


, 125 . 

, 160 . 

, 200 . 

, 250 . 


0115 


1 


315 

. 

, 400 . 

, 500 . 

, 630 . 

, 800 . 


0116 


1 


1000 

. 

, 1250. 

, 1600, 

, 2000. 

, 2500. 


0117 


3 


3150 


, 4000. 

, 5000. 

, 6300. 

, 8000. 


0118 


2 


10000 


, 12 500 . 

, 16 0 0 0 . 

, 20000 . 

,25000 . 


0119 


3 


31500 


, 40000 . 

/ 





0120 

0121 

0122 NAMELIST /INDATA / C ,L ,R 


0123 

1 

THETA 

, PHI 

, ALPSTAR 

0124 

2 

H 

, PSI 

,u 

0125 

1 

ITRIP 

, I LAM 

, ITURB 

0126 

2 

I BLUNT 

, IT I P 

, ROUND 

0127 

3 

ALPHTIP 

, NSEG 

i CO 

0128 

4 

vise 




0129 

0130 

0131 

0132 C READ IN NAMELIST DATA AND ECHO INPUT TO OUTPUT FILE 

0133 C 

0134 

0135 OPEN ( UNI T=4 , STATUS = 'OLD') 

0136 READ ( 4 , INDATA ) 

0137 

0138 OPEN ( UNI T=5 , STATUS * 'NEW') 

0139 WRITE ( 5 , INDATA) 

0140 

0141 

0142 C INITIALIZE ALL PREDICTED PRESSURES AND SOUND 

0143 C PRESSURE LEVELS TO ZERO 

0144 C 

0145 

0146 DO 6001 1=1 , NFREQ 

0147 P 1 ( I ) = 0.0 

0148 P 2 ( I ) = 0.0 

0149 P3 ( I ) = 0.0 

0150 P4(I) = 0.0 
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0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 
0161 
0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 
0181 
0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 
0 201 
0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 
0211 
0212 

0213 

0214 

0215 

0216 
0217 
0216 

0219 

0220 
0221 
0222 

0223 

0224 

0225 


P5 ( I ) = 0.0 
P 6 < I ) = 0.0 
P7 ( I ) = 0.0 

DO 6002 J=l,7 
SPL( J , I ) * 0.0 

6002 CONTINUE 
6001 CONTINUE 


C FOR EACH BLADE SEGMENT, MAKE A NOISE PREDICTION ACCORDING 

C TO THE MECHANISMS SELECTED. TIP NOISE IS PREDICTED FOR 

C THE LAST SEGMENT ONLY. 

C 

DO 6000 1 1 1 = 1 , NS EG 
IF (I LAM . EQ. 1) 

1 CALL LBLVS ( ALPSTAR ( III ) , C ( I II ) f U ( III ) , FRCEN, SPLLBL, 

2 THETA( III ) , PHI (III) r L(IIX) , R ( I I I ) , NFREQ , 

3 VISC, CO) 

IF ( ITURB .EQ. 1) 

1 CALL TBLTE ( ALPSTAR( I II ) , C ( III ) , U( III ) , FRCEN, ITRIP , SPLP f 

1 SPLS , SPLALPH , SPLTBL , THETA (III) , PHI ( III ) ,L(III) ,R( II I ) , 

2 NFREQ , VISC, CO ) 

IF ( I BLUNT .EQ. 1) 

1 CALL BLUNT ( ALPSTAR ( I II ), C( III ), U ( III ), FRCEN , ITRIP ,SPLBLNT , 

1 THETA(III) , PH I ( 1 1 1 ) , L ( 1 1 1 ) f R( III ) , H ( I I I ) , PS I ( 1 1 1 ) , 

2 NFREQ, VISC, CO ) 

IF ( ( IT I P .EQ. 1) .AND. (Ill . EQ . NS EG ) ) 

1 CALL TIPNOIS (ALPHTIP , ALPRAT , C ( III ) ,U< III ), FRCEN, SPLTIP , 

2 THETA, PHI ,R( III ) , NFREQ, VI SC ,C0 , ROUND) 


C ADD IN THIS SEGMENT'S CONTRIBUTION ON A MEAN-SQUARE 

C PRESSURE BASIS 

C 

DO 989 1*1, NFREQ 

IF ( I LAM -EQ. 1) THEN 

P5 ( I ) * P5 ( I ) + 10. **(SPLLBL(I )/10. ) 

ENDIF 

IF (ITURB .EQ. 1) THEN 

P1(I) = P1(I) + 10 . * * ( SPLP ( I )/10 . ) 

P2 ( I ) = P 2 ( I ) + 10 . *M SPLS ( I )/10 . ) 

P3 ( I ) = P 3 ( I ) + 10 .**( SPLALPH ( I )/10 . ) 

ENDIF 

IF ( I BLUNT .EQ. 1) THEN 

P6(I) = P 6 ( I ) + 10 . * * ( S PLBLNT ( I )/10 . ) 

ENDIF 

IF ( ( I TI P .EQ. 1) .AND. (Ill . EQ . NSEG ) ) THEN 
P 7 ( I ) = P 7 ( I ) + 10 . * * ( SPLTIP ( I )/10 . ) 

ENDIF 


C COMPUTE TOTAL PRESSURE FOR THE SEGMENT FOR ALL MECHANISMS 

C 

P4 ( I ) = Pl(I) + P2 ( I ) + P 3 ( I ) + P5 ( I ) + P6 ( I ) + P7 ( I ) 

989 CONTINUE 
6000 CONTINUE 

C CONTRIBUTIONS FROM ALL SEGMENTS ARE NOW ACCOUNTED FOR. 

C COMPUTE SOUND PRESSURE LEVELS FOR EACH MECHANISM AND 

C FOR THE TOTAL 

C 
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0226 DO 6003 1=1 , NFREQ 


0227 

IF 

(pim 

. NE . 

0 . 

SPL (1,1) 

= 

10 . * ALOG 1 0 ( PI ( I ) ) 

0228 

IF 

(P2 ( I ) 

.NE . 

0 . 

S PL ( 2,1) 

= 

10 . * ALOGl 0 ( P 2 ( I ) ) 

0229 

IF 

< P 3 ( I ) 

. NE . 

0 . 

SPL (3,1) 

= 

10 . * ALOG1 0 ( P 3 ( I ) ) 

0230 

IF 

< P4 < X ) 

.NE . 

0 . 

S P L ( 4,1) 

= 

10 . * ALOGl 0 ( P4 ( I ) ) 

0231 

IF 

( P5 ( I ) 

. NE . 

0 . 

SPL { 5,1) 

= 

10 . * ALOGl 0 ( P5 ( I ) ) 

0232 

IF 

(P6 ( I ) 

.NE . 

0 . 

SPL ( 6,1) 

= 

10 . * ALOG 1 0 ( P 6 ( I ) ) 

0233 

IF 

( P 7 ( I ) 

.NE . 

0 . 

S PL ( 7,1) 

= 

10 . * ALOGl 0 ( P7 ( I ) ) 


0234 6003 CONTINUE 

0235 

0236 

0237 C WRITE OUTPUT FILE 

0238 C 

0239 

0240 WRITE (5,7000) 

0241 

0242 DO 6005 1=1, NFREQ 

0243 

0244 WRITE ( 5 , 7100 ) FRCEN ( I ) , ( SPL ( J , I ) , J=1 , 3 ) , ( SPL { J , I ) , J=5 , 7 ) , 

0245 1 S PL ( 4 , I ) 

0246 IF ( MOD (1,5) . EQ . 0) WRITE ( 5 , 7 2 0 0 ) 

0247 6005 CONTINUE 

0248 

0249 7000 FORMAT ( 1H1 , 52X , r ONE-THIRD OCTAVE ' ,/,50X, 'SOUND PRESSURE LEVELS' 


0250 

1 

////,5X, 

r t 

9 

PRESSURE 

9 

9 

0251 

2 


' SUCTION ' , ' 

SEPARATION 

'/ 

0252 

3 

5X 

' FREQUENCY ( HZ ) ' , ' 

SIDE TBL 

9 

0253 

4 


' SIDE TBL ' , ' 

SIDE TBL 

9 

9 

0254 

5 


' LAMINAR ' , ' 

BLUNTNESS 

9 

0255 

6 


' TIP ' , ' 

TOTAL 

9 

0256 

7 


/ , 5X , 8 ( ' 

— ' ) ,/) 



0257 

0258 7100 FORMAT ( 8F15 . 3 ) 

0259 7200 FORMAT ( ' r ) 

0260 8000 FORMAT (13) 

0261 8002 FORMAT (4110 ) 

0262 

0263 

0264 STOP 

0265 END 
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I 



0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 
004 2 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 


SUBROUTINE 

1 

PARAMETER 


LBLVS ( ALPSTAR , C , U , FRCEN , S PLLAM , THETA , PH I 
NFREQ,VISC,C0 ) 

MAXFREQ = 27) 


L , R , 


C 

c ***** VARIABLE DEFINITIONS ***** 

C 


C 

C 

VARIABLE NAME 

DEFINITION 

UNITS 

C 

ALPSTAR 

ANGLE OF ATTACK 

DEGREES 

c 

C 

CHORD LENGTH 

METERS 

c 

CO 

SPEED OF SOUND 

METERS/SEC 

c 

D 

REYNOLDS NUMBER RATIO 

— 

c 

DBARH 

HIGH FREQUENCY DIRECTIVITY 

— 

c 

DELTAP 

PRESSURE SIDE BOUNDARY LAYER 


c 


THICKNESS 

METERS 

c 

DSTRP 

PRESSURE SIDE BOUNDARY LAYER 


c 


DISPLACEMENT THICKNESS 

METERS 

c 

DSTRS 

SUCTION SIDE BOUNDARY LAYER 


c 


DISPLACEMENT THICKNESS 

METERS 

c 

E 

STROUHAL NUMBER RATIO 

— 

c 

FRCEN 

1/3 OCTAVE FREQUENCIES 

HERTZ 

c 

G 1 

SOUND PRESSURE LEVEL FUNCTION 

DB 

c 

G 2 

OVERALL SOUND PRESSURE LEVEL 


c 


FUNCTION 

DB 

c 

G 3 

OVERALL SOUND PRESSURE LEVEL 


c 


FUNCTION 

DB 

c 

ITRIP 

FLAG TO TRIP BOUNDARY LAYER 

— 

c 

L 

SPAN 

METERS 

c 

M 

MACH NUMBER 

— 

c 

NFREQ 

NUMBER OF FREQUENCIES 

— 

c 

O AS P L 

OVERALL SOUND PRESSURE LEVEL 

DB 

c 

PHI 

DIRECTIVITY ANGLE 

DEGREES 

c 

R 

OBSERVER DISTANCE FROM SEGMENT 

METERS 

c 

RC 

REYNOLDS NUMBER BASED ON CHORD 

— 

c 

RC 0 

REFERENCE REYNOLDS NUMBER 

— 

c 

SCALE 

GEOMETRIC SCALING TERM 


c 

S PLLAM 

SOUND PRESSURE LEVEL DUE TO 


c 


LAMINAR MECHANISM 

DB 

c 

STPRIM 

STROUHAL NUMBER BASED ON PRESSURE 


c 


SIDE BOUNDARY LAYER THICKNESS 

— 

c 

ST 1 PR I M 

REFERENCE STROUHAL NUMBER 

— 

c 

STPKPRM 

PEAK STROUHAL NUMBER 

— 

c 

THETA 

DIRECTIVITY ANGLE 

DEGREES 

c 

U 

FREESTREAM VELOCITY 

METERS/SEC 

c 

vise 

KINEMATIC VISCOSITY 

M2/SEC 


DIMENSION STPRIM< MAXFREQ ) , S PLLAM ( MAXFREQ ) , FRCEN { MAXFREQ ) 

REAL L ,M 


C COMPUTE REYNOLDS NUMBER AND MACH NUMBER 

C 

M = U / CO 

RC = U * C/VISC 


C COMPUTE BOUNDARY LAYER THICKNESSES 

C 

CALL THICK(C,U , ALPSTAR , ITRIP , DELTAP T DSTRS , DSTRP ,C0 , VISC ) 


C 


COMPUTE DIRECTIVITY FUNCTION 


0076 C 

0077 

0078 CALL DIRECTH<M, THETA, PHI, DBARH ) 

0079 

0080 
0081 

0082 C COMPUTE REFERENCE STROUHAL NUMBER 

0083 C 

0084 

0085 IF (RC .LE. 1.3E+05) STlPRIM = .18 

0086 I F ( ( RC .GT. 1 . 3 E + 0 5 ) . AND . ( RC . LE . 4 . 0 E + 0 5 ) ) ST 1 PRIM= . 0 0 1 7 56 * RC * * . 3 9 3 1 

0087 IF (RC . GT . 4.0E+05) STlPRIM = .28 

0088 

0089 STPKPRM = 1 0 . * * ( - . 0 4 * ALP STAR ) * STlPRIM 

0090 

0091 

0092 

0093 C COMPUTE REFERENCE REYNOLDS NUMBER 

0094 C 

0095 

0096 IF (ALPSTAR . LE . 3.0) RC 0 = 1 0 . * * ( . 2 1 5 * ALPSTAR+ 4 . 9 7 8 ) 

0097 IF (ALPSTAR .GT. 3.0) RC0=10.**( .120 * ALPSTAR+5 .263) 

0098 

0099 

0100 
0101 

0102 C COMPUTE PEAK SCALED SPECTRUM LEVEL 

0103 C 

0104 

0105 D * RC / RC 0 

0106 

0107 IF (D .LE. .3 237) G 2 = 7 7 . 8 5 2 * ALOG 1 0 ( D ) + 1 5 . 3 2 8 

0108 IF < ( D .GT. . 3 2 37 ) . AND . ( D . LE . .5689)) 

0109 1 G 2 = 65 . 188 *ALOG10 ( D ) + 9.125 

0110 IF <(D .GT. . 5689 ). AND .( D . LE . 1.7579)) 

0111 1 G2 = -11 4.052 * ALOG 1 0 { D ) * * 2 . 

0112 IF ( ( D .GT. 1 . 7579 ) .AND .( D . LE . 3.0889)) 

0113 1 G 2 = -6 5. 18 8 * ALOG1 0 ( D ) + 9 . 1 2 5 

0114 IF (D .GT. 3.0889) G2 =- 7 7 . 8 5 2 * ALOG1 0 ( D ) +1 5 . 3 2 8 

0115 

0116 

0117 G 3 = 171.04 - 3.03 * ALPSTAR 

0118 

0119 SCALE = 10. * ALOG 1 0 ( DE LTAP *M* * 5 * DBARH * L/R * * 2 ) 

0120 
0121 
0122 

0123 C COMPUTE SCALED SOUND PRESSURE LEVELS FOR EACH STROUHAL NUMBER 

0124 C 

0125 

0126 DO 100 1 = 1 , NFREQ 

0127 

0128 STPRIM ( I ) = FRCEN(I) * DELTAP / U 

0129 

0130 E = STPRIM ( I ) / STPKPRM 

0131 

0132 

0133 

0134 1 

0135 

0136 1 

0137 

0138 1 

0139 

0140 

0141 SPLLAM(I) = G1 + G2 + G3 + SCALE 

0142 

0143 100 CONTINUE 

0144 

0145 RETURN 

0146 END 


IF (E . LT . .5974 ) Gl = 3 9 . 8 * ALOG 1 0 ( E ) -1 1 . 1 2 

IF ((E .GE. . 5974 ) .AND. ( E . LE . .8 54 5 )) 

Gl = 98.409 * ALOG 1 0(E) + 2.0 
IF ((E .GE. . 854 5 ) . AND , ( E . LT . 1.17)) 

Gl = -5.0 76+SQRT (2. 484-506. 25*( ALOG 1 0 ( E ) ) * * 2 . ) 
IF ((E .GE. 1.17). AND . ( E . LT . 1.674)) 

Gl = -98.409 * ALOG 1 0(E) + 2.0 
IF (E .GE. 1.674) G 1 =- 3 9 . 8 0 * ALOG 1 0 ( E ) - 1 1 . 1 2 
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0041 

C 

BMINB 

MINIMUM ' B * EVALUATED AT B 

DB 

0042 

C 

BMINBO 

MINIMUM 1 B ' EVALUATED AT BO 

DB 

0043 

C 

BRBO 

INTERPOLATION FACTOR 

DB 

0044 

C 

C 

CHORD LENGTH 

METERS 

0045 

C 

CO 

SPEED OF SOUND 

METERS/SEC 

0046 

C 

DBARH 

HIGH FREQUENCY DIRECTIVITY 

— 

0047 

C 

DBARL 

LOW FREQUENCY DIRECTIVITY 

— 

0048 

C 

DELK1 

CORRECTION TO AMPLITUDE FUNCTION 

DB 

0049 

C 

DELTAP 

PRESSURE SIDE BOUNDARY LAYER THICKNESS METERS 

0050 

C 

DSTRP 

PRESSURE SIDE DISPLACEMENT THICKNESS 

METERS 

0051 

c 

DSTRS 

SUCTION SIDE DISPLACEMENT THICKNESS 

METERS 

0052 

c 

FRCEN 

ARRAY OF CENTERED FREQUENCIES 

HERTZ 

0053 

c 

GAMMA 

USED IN r B ' COMPUTATION 

— 

0054 

c 

GAMMAO 

USED IN ' B ' COMPUTATION 

— 

0055 

c 

ITRIP 

TRIGGER TO TRIP BOUNDARY LAYER 

— 

0056 

c 

Kl 

AMPLITUDE FUNCTION 

DB 

0057 

c 

K2 

AMPLITUDE FUNCTION 

DB 

0058 

c 

L 

SPAN 

METERS 

0059 

c 

M 

MACH NUMBER 

— 

0060 

c 

NFREQ 

NUMBER OF CENTERED FREQUENCIES 

— 

0061 

c 

PHI 

DIRECTIVITY ANGLE 

DEGREES 

0062 

c 

PI 

PRESSURE SIDE PRESSURE 

NT/M2 

0063 

c 

P 2 

SUCTION SIDE PRESSURE 

NT/M2 

0064 

c 

P 4 

PRESSURE FROM ANGLE OF ATTACK 


0065 

c 


CONTRIBUTION 

NT/M2 

0066 

c 

R 

SOURCE TO OBSERVER DISTANCE 

METERS 

0067 

c 

RC 

REYNOLDS NUMBER BASED ON CHORD 

— 

0068 

c 

RDSTRP 

REYNOLDS NUMBER BASED ON PRESSURE 


0069 

c 


SIDE DISPLACEMENT THICKNESS 

— 

0070 

c 

RDSTRS 

REYNOLDS NUMBER BASED ON SUCTION 


0071 

c 


SIDE DISPLACEMENT THICKNESS 

— 

0072 

c 

SPLALPH 

SOUND PRESSURE LEVEL DUE TO ANGLE OF 


0073 

c 


ATTACK CONTRIBUTION 

DB 

0074 

c 

SPLP 

SOUND PRESSURE LEVEL DUE TO PRESSURE 


0075 

c 


SIDE OF AIRFOIL 

DB 
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0076 

0077 

0078 

0079 

0080 
0081 
0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 
0111 
0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 
0121 
0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 


0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 


SPLS 

SPLTBL 

ST P 
STS 
ST 1 

STlPRIM 
ST 2 

STPEAK 

SWITCH 

THETA 

U 

vise 

XCHECK 


PARAMETER (MAXFREQ 


SOUND PRESSURE LEVEL DUE TO SUCTION 
SIDE OF AIRFOIL 

TOTAL SOUND PRESSURE LEVEL DUE TO 
TBLTE MECHANISM 
PRESSURE SIDE STROUHAL NUMBER 
SUCTION SIDE STROUHAL NUMBER 
PEAK STROUHAL NUMBER 
PEAK STROUHAL NUMBER 
PEAK STROUHAL NUMBER 
PEAK STROUHAL NUMBER 
LOGICAL FOR COMPUTATION OF ANGLE 
OF ATTACK CONTRIBUTION 
DIRECTIVITY ANGLE 
VELOCITY 

KINEMATIC VISCOSITY 

USED TO CHECK FOR ANGLE OF ATTACK 
CONTRIBUTION 


27) 


DB 


DB 


DEGREES 

METERS/SEC 

M2/SEC 


DIMENSION SPLTBL ( MAX FREQ ) 

1 S P LALPH ( MAX FREQ ) 

1 STS ( MAXFREQ ) 

LOGICAL SWITCH 
REAL L,M,K1,K2 


, SPLP (MAXFREQ ) 
t STP ( MAXFREQ ) 

, FRCEN (MAXFREQ ) 


f SPLS ( MAXFREQ ) 


RC 

M 


= u * c / vise 

- U / CO 


COMPUTE BOUNDARY LAYER THICKNESSES 

CALL TH I CK ( C , U , ALP STAR , ITRIP , DELTAP , DSTRS , DSTRP , CO r VISC ) 
COMPUTE DIRECTIVITY FUNCTION 


CALL DIRECTL(M, THETA, PHI , DBARL ) 
CALL DIRECTH ( M , THETA, PHI , DBARH ) 


CALCULATE THE REYNOLDS NUMBERS BASED ON PRESSURE AND 
SUCTION DISPLACEMENT THICKNESS 


RDSTRS * DSTRS * U / VISC 
RDSTRP = DSTRP * U / VISC 

DETERMINE PEAK STROUHAL NUMBERS TO BE USED FOR 
' A ’ AND 'B' CURVE CALCULATIONS 


0132 

ST1 

= .02 * M * 

* ( “ . 6 ) 

0133 




0134 

IF 

(ALPSTAR .LE. 1 

.333 ) ST 2 = STl 

0135 

IF 

((ALPSTAR .GT. 

1.333) .AND. (ALPSTAR 

0136 

1 

ST 2 = STl * 10 . * * 

( . 0 054 * ( ALPSTAR— 1 . 3 3 

0137 

IF 

(ALPSTAR .GT. 1 

2.5) ST 2 = 4.72 * ST 

0138 




0139 




0140 

STl 

PRIM = (ST1+ST2 

)/2. 


CALL A0 COMP ( RC , A0 ) 

CALL AOCOMP ( 3 . *RC , A02 ) 

EVALUATE MINIMUM AND MAXIMUM 


CURVES AT A0 


CALL AMIN ( A0 , AMINA0 ) 
CALL AMAX ( A0 , AMAXA0 ) 
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0151 

0152 CALL AMIN ( AO 2 , AM IN AO 2 ) 

0153 CALL AMAX ( AO 2 f AMAXA0 2 ) 

0154 

0155 C COMPUTE 'A' MAX/MI N RATIO 

0156 C 

0157 

0158 ARAO = (20. + AMINAO ) / (AMINAO - AMAXAO ) 

0159 ARAO 2 = (20. + AMINAO 2 ) / ( AMINAO 2 - AMAXAO 2 ) 

0160 

0161 C COMPUTE BO TO BE USED IN 'B' CURVE CALCULATIONS 

0162 C 

0163 

0164 IF (RC . LT . 9.52E+04) BO = .30 

0165 IF ( ( RC .GE. 9 . 5 2 E + 0 4 ) . AND . ( RC . LT . 8.57E+05)) 

0166 1 BO = ( -4 . 48E-13 ) * < RC-8 . 57E+05 ) * *2 . + .56 

0167 IF (RC .GE. 8.57E+05) BO = .56 

0168 

0169 C EVALUATE MINIMUM AND MAXIMUM 'B' CURVES AT BO 

0170 C 

0171 

0172 CALL BMIN( BO , BMINBO ) 

0173 CALL BMAX ( BO , BMAXBO ) 

0174 

0175 C COMPUTE 'B' MAX/MI N RATIO 

0176 C 

0177 

0178 BRB 0 = (20. + BMINBO) / (BMINBO - BMAXBO) 

0179 

0180 C FOR EACH CENTER FREQUENCY , COMPUTE AN 

0181 C r A' PREDICTION FOR THE PRESSURE SIDE 

0182 C 

018 3 

0184 STPEAK = STl 

0185 

0186 DO 100 1=1 , NFREQ 

0187 STP(I) = FRCEN(I) * DSTRP / U 

0188 A = ALOG 1 0 ( STP(I) / STPEAK ) 

0189 CALL AMIN ( A , AMINA ) 

0190 CALL AMAX ( A, AMAXA ) 

0191 AA = AMINA + ARAO * (AMAXA - AMINA) 

0192 

0193 IF (RC .LT. 2.47E+05) K1 = -4.31 * ALOGIO(RC) + 156.3 

0194 I F ( ( RC .GE. 2 . 4 7 E + 0 5 ) . AND . ( RC . LT . 8.0E + 05)) 

0195 1 Kl = -9.0 * ALOG 1 0 ( RC ) + 181.6 

0196 IF (RC .GT. 8.0E+05) Kl = 128.5 

0197 

0198 IF ( RDSTRP . LE . 5000 .) DELK1 = - ALPSTAR * ( 5 . 2 9 - 1 . 4 3 * 

0199 1 ALOGIO ( RDSTRP ) ) 

0200 IF (RDSTRP . GT . 5000.) DELK1 = 0.0 

0201 

0202 SPLP(I) = AA + Kl — 3 . +10 . * ALOG 10 (DSTRP*M**5. * DBARH * L/R * * 2 . )+ DELHI 

0203 

0204 

0205 

0206 

0207 GAMMA = 27.094 * M + 3.31 

0208 BETA = 72.650 * M + 10.74 

0209 GAMMAO = 23.430 * M + 4.651 

0210 BETAO =-34.190 * M - 13.820 

0211 

0212 IF (ALPSTAR . LE . ( GAMMAO -GAMMA ) ) K2 = -1000.0 

0213 IF ( (ALPSTAR. GT. ( GAMMAO -GAMMA ) ) .AND. (ALPSTAR. LE. ( GAMMAO+GAMMA ) ) ) 

0 214 1 K 2 = SQRT ( BETA* *2 . - ( BETA/GAMMA ) * *2 . * ( ALPSTAR-GAMMAO ) * *2 . ) +BETA0 

0215 IF (ALPSTAR . GT . (GAMMAO+GAMMA)) K2 = -12.0 

0216 

0217 K2 = K2 + Kl 

0218 

0219 

0220 

0221 STS(I) = FRCEN(I) * DSTRS / U 

0222 

0223 C CHECK FOR 'A' COMPUTATION FOR SUCTION SIDE 

0224 C 

0225 
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0226 

0227 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0258 

0259 

0260 
0261 
0262 

0263 

0264 

0265 

0266 

0267 

0268 

0269 

0270 

0271 

0272 

0273 

0274 

0275 

0276 

0277 

0278 

0279 

0280 
0281 
0282 

0283 

0284 
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XCHECK = GAMMAO 
SWITCH = .FALSE. 

IF ((ALPSTAR .GE. XCHECK ). OR .( ALPSTAR . GT . 1 2 . 5 ) ) SWI TCH= . TRUE . 

IF (.NOT. SWITCH) THEN 

A = ALOG 1 0 ( STS ( I ) / ST1PRIM ) 

CALL AMIN ( A, AMINA) 

CALL AMAX ( A, AMAXA) 

AA = AMINA + ARAO * (AMAXA - AMINA) 

SPLS(I) - AA+Kl - 3 . + 1 0 . * ALOG 10(DSTRS*M**5. * DBARH * 

1 L/R* *2 . ) 

C 'B' CURVE COMPUTATION 

C 

B = ABS ( ALOG10 ( STS ( I ) / ST 2 ) ) 

CALL BMIN ( B , BMINB ) 

CALL BMAX ( B , BMAXB ) 

BB = BMINB + BRB0 * ( BMAXB-BMINB ) 

SPLALPH ( I ) =BB + K2 + 1 0 . *ALOG10 ( DSTRS *M * * 5 . * DBARH * L/R * * 2 . ) 

ELSE 


C THE 'A' COMPUTATION IS DROPPED IF 'SWITCH' IS TRUE 

C 


1 

1 


SPLS(I) * 0.0 + 10 . *ALOG10(DSTRS*M**5. * DBARL* 
L/R* *2 . ) 

SPLP(I) = 0.0 + 10.* ALOG 10(DSTRS*M**5. * DBARL * 
L/R* *2 . ) 

B - ABS { ALOG10 ( STS ( I ) / ST 2 ) ) 

CALL AM I N ( B , AM I N B ) 

CALL AMAX ( B , AMAXB ) 

BB = AMI NB + ARAO 2 * ( AMAXB— AMINB ) 

SPLALPH ( I ) =BB + K2 + 1 0 . * ALOG 10 ( DSTRS *M* *5 . * DBARL* 
L/R* *2 . ) 


ENDIF 


C SUM ALL CONTRIBUTIONS FROM r A' AND 'B f ON BOTH 

C PRESSURE AND SUCTION SIDE ON A MEAN-SQUARE PRESSURE 

C BASIS 

C 


IF 

( SPLP ( I ) 

. LT . 

-100 . 

SPLP ( I ) 

= -100 

IF 

( SPLS ( I ) 

. LT . 

-100 . 

SPLS ( I ) 

- -100 

IF 

(SPLALPH ( I ) 

. LT . 

-100 . 

) SPLALPH ( I ) 

= -100 


PI = 10 . * * { S PLP ( I ) / 10 . ) 

P2 = 10 . * * ( SPLS ( I ) / 10 . ) 

P4 = 10 .**( SPLALPH ( I ) / 10.) 

SPLTBL(I) a 10. * ALOG 1 0 ( P 1 + P2 + P4 ) 

100 CONTINUE 

RETURN 

END 
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0001 SUBROUTINE AMIN ( A , AMINA ) 

0002 

0003 C THIS SUBROUTINE DEFINES THE CURVE FIT CORRESPONDING 

0004 C TO THE A-CURVE FOR THE MINIMUM ALLOWED REYNOLDS NUMBER. 

0005 C 

0006 

0007 XI = ABS(A) 

0008 

0009 IF (XI .LE. .204 ) AMINA=SQRT ( 6 7 . 5 5 2 - 8 8 6 . 7 8 8 *X1 * * 2 . ) -8 . 2 1 9 

0010 I F ( ( X 1 .GT. . 204 ) . AND. ( XI . LE . . 2 4 4 ) ) AMI NA=- 3 2 . 6 6 5 * X 1 + 3 . 9 8 1 

0011 IF (XI. GT. .244 ) AMINA=-142 . 795*Xl**3.+103.656*Xl**2.-57.757*Xl+6.006 

0012 

0013 RETURN 

0014 END 


0001 SUBROUTINE AMAX ( A , AMAXA ) 

0002 

0003 C THIS SUBROUTINE DEFINES THE CURVE FIT CORRESPONDING 

0004 C TO THE A-CURVE FOR THE MAXIMUM ALLOWED REYNOLDS NUMBER. 

0005 

0006 XI = ABS(A) 

0007 

0008 IF (XI .LE. .13) AMAXA= SQRT (67. 552-886. 788*Xl**2. ) — 8 . 219 

0009 I F ( ( X 1 .GT. .13). AND. (XI . LE . . 3 2 1 ) ) AMAXA=- 1 5 . 9 0 1 * X 1 + 1 . 0 9 8 

0010 IF (XI. GT. .321 ) AMAX A=- 4 . 6 6 9 * X 1 * * 3 . + 3 . 4 9 1 *X 1 * * 2 . - 1 6 . 6 9 9 *X1 + 1 . 1 4 9 

0011 

0012 RETURN 

0013 END 


0001 SUBROUTINE BMI N ( B , BMI NB ) 

0002 

0003 C THIS SUBROUTINE DEFINES THE CURVE FIT CORRESPONDING 

0004 C TO THE B-CURVE FOR THE MINIMUM ALLOWED REYNOLDS NUMBER. 

0005 

0006 XI = ABS ( B ) 

0007 

0008 IF (XI .LE. . 13 )BMINB = SQRT( 16 . 888-886 . 788 *Xl**2 . )-4 . 109 

0009 I F ( (XI .GT. .13). AND. (XI . LE . . 1 4 5 ) ) BMI NB-- 8 3 . 6 0 7 * X 1 + 8 . 1 3 8 

0010 IF (Xl.GT. .145) BMINB=-8 1 7.81*X1**3.+355.21*X1**2.-135.024*X1+10.619 

0011 

0012 RETURN 

0013 END 


0001 SUBROUTINE BMAX ( B , BMAXB ) 

0002 

0003 C THIS SUBROUTINE DEFINES THE CURVE FIT CORRESPONDING 

0004 C TO THE B-CURVE FOR THE MAXIMUM ALLOWED REYNOLDS NUMBER. 

0005 

0006 XI - ABS(B) 

0007 

0008 IF (XI .LE. .1) BMAXB=SQRT (16. 888-886. 788*X1**2. ) - 4 . 1 0 9 

0009 I F ( (XI .GT. .1). AND. (XI . LE . .187) ) BMAXB = - 3 1 . 3 1 3 * Xl + 1 . 8 5 4 

0010 IF (Xl.GT. .187) BMAXB=-80 . 541*Xl**3.+44.174*Xl**2.-39.381*Xl+2.344 

0011 

0012 RETURN 

0013 END 
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0001 

0002 

0003 

C 

SUBROUTINE A0 COMP ( RC , A0 ) 

THIS SUBROUTINE DETERMINES WHERE 

THE A-CURVE 

0004 

C 

TAKES ON A 

VALUE OF -20 dB . 


0005 

0006 
0007 


IF (RC . LT . 
IF ( ( RC . GE 

9 . 5 2 E + 0 4 ) A0 = . 57 

. 9 . 52E + 0 4 ) . AND . ( RC . 

LT. 8 . 57E + 05 ) ) 

0008 


1 A0 = (-9 

. 57E-13 ) * (RC-8 . 57E+05 

) **2 . + 1.13 

0009 

0010 
0011 


IF (RC .GE. 

RETURN 

END 

8 . 57E + 05 ) A0 = 1.13 



0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 


SUBROUTINE DIRECTH (M, THETA, PHI , DBAR ) 

C THIS SUBROUTINE COMPUTES THE HIGH FREQUENCY 

C DIRECTIVITY FUNCTION FOR THE INPUT OBSERVER LOCATION 

REAL M,MC 

DEGRAD = .017453 

MC ~ . 8 * M 

THETAR = THETA * DEGRAD 
PHIR = PHI * DEGRAD 

DBAR= 2 . *SIN ( THETAR/ 2 . )**2.*SIN(PHIR)**2./( (1. +M*COS ( THETAR ) ) * 
1 (l. + {M-MC)*COS( THETAR ) ) * * 2 . ) 

RETURN 

END 


0001 


SUBROUTINE DIRE CTL(M, THETA, PHI, DBAR ) 

0002 




0003 

C 

THIS SUBROUTINE COMPUTES THE 

LOW FREQUENCY 

0004 

C 

DIRECTIVITY FUNCTION FOR THE 

INPUT OBSERVER LOCATION 

0005 




0006 


REAL M , MC 


0007 




0008 


DEGRAD = .017453 


0009 




0010 


MC = . 8 * M 


0011 


THETAR = THETA * DEGRAD 


0012 


PHIR = PHI * DEGRAD 


0013 




0014 


DBAR = ( SIN ( THETAR ) *SIN ( PHIR ] 

1 ) ** 2 /( 1 . +M* COS (THETAR ) ) *M 

0015 




0016 


RETURN 


0017 


END 
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0001 

0002 

0003 

0004 


SUBROUTINE BLUNT ( ALPSTAR , C r U f FRCEN , I TRI P f SPLBLNT , THETA , PHI , 
1 L , R , H , PS I ,NFREQ f VISC,C0 ) 

0005 

C 




0006 

C 

***** 

VARIABLE DEFINITIONS ***** 


0007 

0008 

C 




0009 

C 

VARIABLE NAME 

DEFINITION 

UNITS 

0010 

0011 

c 




0012 

c 

ALPSTAR 

ANGLE OF ATTACK 

DEGREES 

0013 

c 

ATERM 

USED TO COMPUTE PEAK STROUHAL NO. 

— 

0014 

c 

C 

CHORD LENGTH 

METERS 

0015 

c 

CO 

SPEED OF SOUND 

METERS/SEC 

0016 

c 

DBARH 

HIGH FREQUENCY DIRECTIVITY 

— 

0017 

c 

DELTAP 

PRESSURE SIDE BOUNDARY LAYER 


0018 

c 


THICKNESS 

METERS 

0019 

c 

DSTARH 

AVERAGE DISPLACEMENT THICKNESS 


0020 

c 


OVER TRAILING EDGE BLUNTNESS 

— 

0021 

c 

DSTRAVG 

AVERAGE DISPLACEMENT THICKNESS 

METERS 

0022 

c 

DSTRP 

PRESSURE SIDE DISPLACEMENT THICKNESS 

METERS 

0023 

c 

DSTRS 

SUCTION SIDE DISPLACEMENT THICKNESS 

METERS 

0024 

c 

ETA 

RATIO OF STROUHAL NUMBERS 

— 

0025 

c 

FRCEN 

ARRAY OF 1/3 OCTAVE CENTERED FREQ. 

HERTZ 

0026 

c 

F4TEMP 

G 5 EVALUATED AT MINIMUM HDSTARP 

DB 

0027 

c 

G 4 

SCALED SPECTRUM LEVEL 

DB 

0028 

c 

G 5 

SPECTRUM SHAPE FUNCTION 

DB 

0029 

c 

G 5 0 

G 5 EVALUATED AT PSI = 0.0 

DB 

0030 

c 

G 5 1 4 

G 5 EVALUATED AT PSI = 14.0 

DB 

0031 

c 

H 

TRAILING EDGE BLUNTNESS 

METERS 

0032 

c 

HDSTAR 

BLUNTNESS OVER AVERAGE DISPLACEMENT 


0033 

c 


THICKNESS 

— 

0034 

c 

HDSTARL 

MINIMUM ALLOWED VALUE OF HDSTAR 

— 

0035 

c 

HDSTARP 

MODIFIED VALUE OF HDSTAR 

— 

0036 

c 

ITRIP 

TRIGGER FOR BOUNDARY LAYER TRIPPING 

— 

0037 

c 

L 

SPAN 

METERS 

0038 

c 

M 

MACH NUMBER 

— 

0039 

c 

NFREQ 

NUMBER OF CENTERED FREQUENCIES 

— 

0040 

c 

PHI 

DIRECTIVITY ANGLE 

DEGREES 

0041 

c 

PSI 

TRAILING EDGE ANGLE 

DEGREES 

0042 

c 

R 

SOURCE TO OBSERVER DISTANCE 

METERS 

0043 

c 

RC 

REYNOLDS NUMBER BASED ON CHORD 

— 

0044 

c 

SCALE 

SCALING FACTOR 

— 

0045 

c 

SPLBLNT 

SOUND PRESSURE LEVELS DUE TO 


0046 

c 


BLUNTNESS 

DB 

0047 

c 

STPEAK 

PEAK STROUHAL NUMBER 

— 

0048 

c 

STPPP 

STROUHAL NUMBER 

— 

0049 

c 

THETA 

DIRECTIVITY ANGLE 

— 

0050 

c 

U 

FREESTREAM VELOCITY 

METERS/SEC 

0051 

0052 

0053 

0054 

0055 

c 

vise 

PARAMETER (MAXFREQ 

KINEMATIC VISCOSITY 
= 27) 

M2/SEC 

0056 

0057 

0058 

0059 


DIMENSION SPLBLNT ( MAXFREQ ) , FRCEN ( MAXFREQ ) , STPPP ( MAXFREQ ) 

REAL M , L 

0060 

c 

COMPUTE NECESSARY 

QUANTITIES 


0061 

0062 

0063 

0064 

0065 

0066 

c 

M = U /CO 

RC = U * C / VISC 



0067 

c 

COMPUTE BOUNDARY LAYER THICKNESSES 


0066 

0069 

0070 

0071 

c 

CALL TH I CK ( C , U , ALPSTAR , ITRIP , DELTAP , DSTRS , DSTRP , CO , VISC ) 


0072 

c 

COMPUTE AVERAGE DISPLACEMENT THICKNESS 


0073 

0074 

0075 

c 

DSTRAVG = (DSTRS + 

DSTRP) / 2. 




0076 HDSTAR = H / DSTRAVG 

0077 

0078 DSTARH a 1. /HDSTAR 

0079 

0080 C COMPUTE DIRECTIVITY FUNCTION 

0081 C 

0082 

0083 CALL DIRECTH <M f THETA , PHI , DBARH ) 

0084 

0085 

0086 C COMPUTE PEAK STROUHAL NUMBER 

0087 C 

0088 

0089 ATERM = .212 - .0045 * PSI 

0090 

0091 IF (HDSTAR .GE. .2) 

0092 1 STPEAK = ATERM / < 1 . + . 2 3 5 * DSTARH- . 0 1 3 2 * DSTARH * * 2 . ) 

0093 IF (HDSTAR . LT . .2) 

0094 1 STPEAK = .1 * HDSTAR + .095 - .00243 * PSI 

0095 

0096 C COMPUTE SCALED SPECTRUM LEVEL 

0097 C 

0098 

0099 IF (HDSTAR . LE . 5.) G4 = l 7 . 5 * ALOG1 0 ( HDSTAR ) +1 5 7 . 5-1 . 1 1 4 * PS I 

0100 IF (HDSTAR . GT . 5.) G4=169.7 - 1.114 * PSI 

0101 
0102 

0103 C FOR EACH FREQUENCY, COMPUTE SPECTRUM SHAPE REFERENCED TO 0 DB 

0104 C 

0105 

0106 DO 1000 1=1, NFREQ 

0107 

0108 STPPP(I) = FRCEN(I) * H / U 

0109 ETA = ALOG10 ( STPPP ( I )/STPEAK) 

0110 

0111 HDSTARL = HDSTAR 

0112 

0113 CALL G5COMP ( HDSTARL , ETA, G514 ) 

0114 

0115 HDSTARP = 6.724 * HDSTAR * * 2 . - 4 . 0 1 9 * H DSTAR+ 1 . 1 0 7 

0116 

0117 CALL G5COMP ( HDSTARP , ETA, G50) 

0118 

0119 

0120 G 5 = G50 + .0714 * PSI * (G514-G50) 

0121 IF { G5 . GT . 0 . ) G5 = 0 . 

0122 CALL G5COMP ( . 25 , ETA, F4TEMP ) 

0123 IF ( G 5 . GT . F 4 TEMP ) G5 = F4 TEMP 

0124 

0125 

0126 SCALE s* 10. * ALOG 1 0 ( M * * 5 . 5 *H * DBARH * L/R * * 2 . ) 

0127 

0128 S PLBLNT ( I ) = G4 + G5 + SCALE 

0129 

0130 

0131 1000 CONTINUE 

0132 

0133 RETURN 

0134 END 
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0001 SUBROUTINE G 5COMP { HDSTAR , ETA , G 5 ) 

0002 

0003 

0004 REAL M,K,MU 

0005 

0006 

0007 IF (HDSTAR . LT . .25) MU = .1211 

0008 IF ((HDSTAR . GT . . 2 5 ) . AND . ( HDST AR 

0009 1 MU=- . 2175 *HDSTAR + .1755 

0010 IF ((HDSTAR .GT. . 6 2 ) . AND . ( HDSTAR 

0011 1 MU = -.0308 * HDSTAR + .0596 

0012 IF (HDSTAR .GE. 1 . 1 5 ) MU = .0242 

0013 

0014 IF (HDSTAR . LE . .02) M = 0.0 

0015 IF ((HDSTAR . GE . . 0 2 ). AND .( HDSTAR . LT . .5)) 

0016 1 M=68 . 724 *HDSTAR - 1.35 

0017 IF ((HDSTAR . GT . . 5 ) . AND . ( H DSTAR . LE . .62)) 

0018 1 M = 308.475 * HDSTAR - 121.23 

0019 IF ((HDSTAR .GT. . 6 2 ) . AND . ( HDSTAR . LE . 1.15)) 

0020 1 M = 224.811 * HDSTAR - 69.354 

0021 IF ((HDSTAR . GT . 1.15) .AND. (HDSTAR . LT . 1.2)) 

0022 1 M = 1583.28 * HDSTAR - 1631.592 

0023 IF (HDSTAR . GT . 1.2) M = 268.344 

0024 IF { M . LT . 0.0) M = 0.0 

0025 

0026 ETAO = -SQRT ( ( M*M*MU* *4 )/( 6 . 2 5+M*M*MU*MU ) ) 

0027 

0028 K = 2 . 5 * SQRT ( 1 .-( ETAO /MU ) * * 2 . ) - 2 . 5 -M* ETAO 

0029 

0030 IF (ETA .LE. ETAO) G5 = M * ETA + K 

0031 IF ((ETA .GT. ETAO ). AND .( ETA . LE . 0 . ) ) G 5= 2 . 5 * SQRT ( 1 . - ( ETA/MU ) * * 2 . ) - 2 . 5 

00 3 2 IF ( ( ETA.GT. 0 . ) . AND. ( ETA. LE . . 03616 ) )G5 = SQRT( 1 . 5625-1194. 99*ETA**2. )-1.25 

0033 IF (ETA . GT . .03616) G5=-155.543 * ETA + 4.375 

0034 

0035 RETURN 

0036 END 


. LE . . 62 ) ) 

. LT . 1.15)) 
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0001 


SUBROUTINE 

TIPNOIS ( ALPHTIP , ALPRAT , C , U , FRCEN , 

SPLTIP ,THETA r PHI , 

0002 

0003 

0004 

C 

1 

R, NFREQ, VISC, CO , ROUND) 



0005 

C 


***** VARIABLE DEFINITIONS ***** 



0006 

0007 

C 





0008 

C 

VARIABLE 

NAME DEFINITION 


UNITS 

0009 

c 





0010 

0011 

c 

ALPHTIP 

TIP ANGLE OF ATTACK 


DEGREES 

0012 

c 

ALPRAT 

TIP LIFT CURVE SLOPE 


— 

0013 

c 

ALPTIPP 

CORRECTED TIP ANGLE OF ATTACK 

DEGREES 

0014 

c 

C 

CHORD LENGTH 


METERS 

0015 

c 

CO 

SPEED OF SOUND 


METERS/SEC 

0016 

c 

DBARH 

DIRECTIVITY 


— 

0017 

c 

FRCEN 

CENTERED FREQUENCIES 


HERTZ 

0018 

c 

L 

CHARACTERISTIC LENGTH FOR TIP 

METERS 

0019 

c 

M 

MACH NUMBER 


— 

0020 

c 

MM 

MAXIMUM MACH NUMBER 


— 

0021 

c 

NFREQ 

NUMBER OF CENTERED FREQUENCIES 

— 

0022 

c 

PHI 

DIRECTIVITY ANGLE 


DEGREES 

0023 

c 

R 

SOURCE TO OBSERVER DISTANCE 

METERS 

0024 

c 

ROUND 

LOGICAL SET TRUE IF TIP 

IS ROUNDED 

0025 

c 

SCALE 

SCALING TERM 


— 

0026 

c 

SPLTIP 

SOUND PRESSURE LEVEL DUE 

TO TIP 


0027 

c 


MECHANISM 


DB 

0028 

c 

STPP 

STROUHAL NUMBER 


— 

0029 

c 

TERM 

SCALING TERM 


— 

0030 

c 

THETA 

DIRECTIVITY ANGLE 


DEGREES 

0031 

c 

U 

FREESTREAM VELOCITY 


METERS/SEC 

0032 

c 

UM 

MAXIMUM VELOCITY 


METERS/SEC 

0033 

0034 

c 

vise 

KINEMATIC VISCOSITY 


M2/SEC 

003 5 
0036 


PARAMETER 

( MAXFREQ =27) 



0037 


DIMENSION 

SPLTIP (MAXFREQ ) , FRCEN ( MAXFREQ ) 



003 8 


REAL L , M , MM 



0039 


LOGICAL ROUND 



0040 

0041 

0042 


ALPTIPP ■ 

ALPHTIP * ALPRAT 



0043 

0044 


M 

U / CO 



0045 


CALL DIRECTH { M , THETA, PHI , DBARH ) 



0046 

0047 


IF (ROUND) 

THEN 



0048 


L = .008 

* ALPTIPP * C 



0049 


ELSE 




0050 


IF { ABS ( 

ALPTIPP ) . LE . 2 . ) THEN 



0051 


L = ( . 

023 + . 0169 *ALPTIPP ) * C 



0052 


ELSE 




0053 


L = ( . 

0378 + . 0095*ALPTIPP ) * C 



0054 


ENDIF 




0055 


ENDIF 




0056 

0057 

0058 

0059 


MM = ( 

1. + . 036 *ALPTIPP ) * M 



0060 


UM = MM * CO 



0061 

0062 


TERM = M* 

M * MM **3 . *L**2 . * DBARH/R * * 2 . 



0063 


IF (TERM . 

NE . 0.0) THEN 



0064 


SCALE = 

1 0 . * ALOG 1 0 ( TERM ) 



0065 


ELSE 




0066 


SCALE = 

o 

o 



0067 


ENDIF 




0068 

0069 


DO 100 1=1 

, NFREQ 



0070 


STPP 

= FRCEN ( I ) * L / UM 



0071 


SPLTIP ( I 

) = 1 26 . -30 . 5* < ALOG1 0 ( STPP ) + . 3 ) * * 2 . 

+ SCALE 


0072 


100 CONTINUE 




0073 


RETURN 




0074 


END 
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0001 

0002 


SUBROUTINE THICK(C,U , ALP STAR , ITRI P , DELTAP , DSTRS , DSTRP , CO , VI SC ) 


0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 
00 51 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 


UOUj 

c 


— 

— 



0004 

C 


* * 

* * *. 

VARIABLE DEFINITIONS ***** 


0005 

c 


-- 

— 



0006 







0007 

c 

VARIABLE 

NAME 


DEFINITION 

UNITS 

0008 

c 






0009 







0010 

c 

ALPSTAR 



ANGLE OF ATTACK 

DEGREES 

0011 

c 

C 



CHORD LENGTH 

METERS 

0012 

c 

CO 



SPEED OF SOUND 

METERS/SEC 

0013 

c 

DELTAO 



BOUNDARY LAYER THICKNESS AT 


0014 

c 




ZERO ANGLE OF ATTACK 

METERS 

0015 

c 

DELTAP 



PRESSURE SIDE BOUNDARY LAYER 


0016 

c 




THICKNESS 

METERS 

0017 

c 

DSTRO 



DISPLACEMENT THICKNESS AT ZERO 


0018 

c 




ANGLE OF ATTACK 

METERS 

0019 

c 

DSTRP 



PRESSURE SIDE DISPLACEMENT 


0020 

c 




THICKNESS 

METERS 

0021 

c 

DSTRS 



SUCTION SIDE DISPLACEMENT 


0022 

c 




THICKNESS 

METERS 

0023 

c 

ITRIP 



TRIGGER FOR BOUNDARY LAYER TRIPPING 

0024 

c 

M 



MACH NUMBER 

— 

0025 

c 

RC 



REYNOLDS NUMBER BASED ON CHORD 

— 

0026 

c 

U 



FREESTREAM VELOCITY 

METERS/SEC 

0027 

c 

vise 



KINEMATIC VISCOSITY 

M2/SEC 

0028 







0029 







0030 

c 

COMPUTE ZERO ANGLE 

OF ATTACK BOUNDARY LAYER 


0031 

c 

THICKNESS 

( METERS ) 

AND REYNOLDS NUMBER 


0032 

c 



— 



0033 







0034 


M = 

u / 

CO 



0035 







0036 


RC 

U * 

C/VISC 


0037 







0038 


DELTAO 

10 . * 

Ml. 

6 56 9-. 9 0 4 5 * ALOG 1 0 (RC) + 


0039 


1 

. 0 596 * ALOG 10 ( RC ) * * 2 . ) *C 


0040 


IF (ITRIP 

• EQ. 

2) DELTAO = .6 * DELTAO 



COMPUTE PRESSURE SIDE BOUNDARY LAYER THICKNESS 


DELTAP ■ 10. **( — .04175 * ALPSTAR+ .00106 * ALP STAR * * 2 . ) * DELTAO 


COMPUTE ZERO ANGLE OF ATTACK DISPLACEMENT THICKNESS 


IF ((ITRIP .EQ. 1) .OR. (ITRIP . EQ . 2)) THEN 

IF (RC .LE. . 3 E+ 0 6 ) DSTR0 = .0601 * RC **(-.114)*C 

IF (RC .GT. . 3 E + 0 6 ) 

1 DSTR0=10. **(3.411-1. 5 3 9 7 * ALOG 1 0 ( RC ) + . 10 59 *ALOG10 ( RC > * * 2 . ) *C 

IF (ITRIP .EQ. 2) DSTRO = DSTRO * .6 

ELSE 

DSTR 0 = 10. **(3.018 7-1. 5397 * ALOG 10 ( RC ) + . 10 59 *ALOG10 ( RC ) * *2 . ) *C 
ENDIF 

PRESSURE SIDE DISPLACEMENT THICKNESS 


DSTRP = 10 . * * ( - . 04 32 *ALPSTAR+ . 001 1 3 *ALPSTAR* * 2 . ) *DSTR0 
IF (ITRIP .EQ. 3) DSTRP = DSTRP * 1.48 

SUCTION SIDE DISPLACEMENT THICKNESS 


IF (ITRIP .EQ. 1) THEN 

IF (ALPSTAR .LE. 5.) DSTRS= 1 0 . * * { . 0 6 7 9 * ALPSTAR ) * DSTRO 
I F ( ( ALPSTAR .GT. 5 .). AND .( ALPSTAR . LE . 12.5)) 

1 DSTRS = . 3 8 1 * 10 . * * ( . 1 51 6 *ALPSTAR ) *DSTR0 

IF (ALPSTAR .GT. 1 2 . 5 ) DSTRS= 1 4 . 2 9 6 * 1 0 . * * ( . 0 2 5 8 * ALPSTAR ) * DSTRO 
ELSE 
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0076 

IF (ALPSTAR 

.LE. 7.5) DSTRS =10 - * * ( .0679 * ALPSTAR ) 

* DSTRO 

0077 

IF ( (ALPSTAR 

.GT. 7 . 5 ) .AND .( ALPSTAR . LE . 12.5)) 


0078 

1 DSTRS = .0162*10. **( .3066 * ALPSTAR ) * DSTRO 


0079 

IF (ALPSTAR 

.GT. 12.5) DSTRS = 5 2 . 4 2 * 1 0 . * * ( . 0 2 5 8 

* ALPSTAR ) * DSTRO 

0080 

ENDIF 



0081 




0082 

RETURN 



0083 

END 
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Table 1. Test and Data Processing Conditions for Tripped Boundary-Layer Airfoil Cases of Figures 11 to 43 
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Table 2. Test and Data Processing Conditions for Untripped Boundary- Layer Airfoil Cases of Figures 44 to 74 
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Table 2. Concluded 
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